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ABSTRACT 

At high-energy 7 -rays (> 100 MeV) the Large Area Telescope (LAT) on the Fermi satellite 
already detected more than 145 rotation-powered pulsars (RPPs), while the number of pul¬ 
sars seen at soft 7 -rays (20 keV - 30 MeV) remained small. We present a catalogue of 18 
non-recycled RPPs from which presently non-thermal pulsed emission has been securely de¬ 
tected at soft 7 -rays above 20 keV, and characterize their pulse profiles and energy spectra. 
For 14 of them we report new results, (re)analysing mainly data from RXTE, INTEGRAL, 
XMM-Newton and Chandra. The soft 7 -pulsars are all fast rotators and on average ^ 9.3x 
younger and ~ 43x more energetic than the Fermi LAT sample. The majority (11 members) 
exhibits broad, structured single pulse profiles, and only 6 have double (or even multiple. 
Vela) pulses. Fifteen soft 7 -ray pulsar show hard power-law spectra in the hard X-ray band 
and reach maximum luminosities typically in the MeV range. For only 7 of the 18 soft 7 -ray 
pulsars pulsed emission has also been detected by the LAT, but 12 have a pulsar wind nebula 
(PWN) detected at TeV energies. For six pulsars with PWNe, we present also the spectra of the 
total emissions at hard X-rays, and for IGRJ18490-0000, associated with HESS J1849-000 
and PSRJ1849-0001, we used our Chandra data to resolve and characterize the contribu¬ 
tions from the point-source and PWN. Finally, we also discuss a sample of 15 pulsars which 
are candidates for future detection of pulsed soft 7 -rays, given their characteristics at other 
wavelengths. 

Key words: radiation mechanisms: non-thermal - stars: neutron - X-rays: general - 
gamma-rays: general - pulsars: individual: PSRJ0205H-6449, PSRB0531H-21, PSRJ0537- 
6910, PSRB0540-69, PSRB0833-45, IGR J14003-6326, PSRB1509-58, PSRJ1617-5055, 
PSRJ 1640-4631, PSRJ1811-1925, PSRJ1813-1246, PSRJ1813-1749, AXJ1838.0-0655, 
PSRJ1846-0258, IGR J18490-0000, PSR J1930 h-1852, PSR J2022h-3842, PSRJ2229h- 6114, 
PSR Ill 19-6127, PSR J1124-5916, PSR J1357-6429, PSR J1833-1034, PSR J1420-6048, 
PSR J1418-6058, PSR J1023-5746, PSR J1838-0537, PSR J1826-1256, PSR J1301-6305, 
PSR J1341-6220, PSR J1614-5048, PSR J1803-2137, PSR J1856 h-0245, IGR J11014-6103 


1 INTRODUCTION 

Rotation-powered pulsars (RPPs) form a well-established galactic 
7 -ray source population known to emit steadily from radio frequen¬ 
cies up to high-energy 7 -rays. Despite decades of theoretical mod¬ 
eling the origin and nature of the non-thermal high-energy (> 10 
keV) emission are still being debated. The most developed mod¬ 
els assume as production site of the high-energy emission either 
a location in the pulsar magnetosp here near (or starting near) the 
magneti c pole (polar cap mod els IHarding et al.ir200 ^: s lot gap 
models dMuslimov & H^din3 12004 : Dvks & Rudakl 2003h ~) or in 
vacuum gaps along the last cl osed field lines in the outer mag¬ 
netosphere (outer gap models dCheng et al.lll98(d : [Romanll l l996l : 
ICheng et al.ll 2000 h l. or e ven in the stripe of the pulsar wind outside 
the light cylinder radius dPetri et alj|2012al) . Modeling its emission 
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characteristics (e.g. phase-resolved spectra, pulse morphology and 
polarization ) drives ever-more sophi s ticated electrodynamic calcu¬ 
lations (e.g. IWang & Hirot anill201 ll: i KalTOOth p-akos et ^l2012al : 


iLi et alj[201 J P 6 tri et aLlbOOhl !^ HirotanilbOlS 1. Simulations have 

been performed for the different geometrical models to compute 
beaming patterns and light curves as a function of magnetic inclina¬ 
tion angle and viewing dire ction, resulting in atlasses of gamma-ray 
profil e s and radio lags (e. g. Watters et alj2009l: Romanj_& Wattosi 
I2OIOI : iBai & Spitkwskvl I2010al f9: Kalapotharakos et"^ l 2012 bh . 
These recent theoretical developents were driven by the major ob¬ 
servational progress achieved by the Large Area Telescope (LAT; 
20 MeV - 300 GeV ) aboard the Fermi Gamma-ray Space Tele¬ 
scope, launched June 11, 2008. After more than four years of oper¬ 
ations Fermi LAT had detected and characterized the pulsed emis- 
sion (see e.g. lAbdo et alj20103 : lRav et alj201 lal : lAbdo et alj2013l) 
of about 145 pulsars (status May 8 , 2013), a significant fraction 
of these is 7 -ray only i.e. Geminga-like with no obvious counter- 
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parts at other (less energetic) wavelengths regimes. This number 
can be compared with only seven pulsars securely detected above 
100 MeV with EGRET on the preceding Compton Gamma-Ray 
Observatory (CGRO; April 5, 1991 - June 4, 2000). 

The observational progress at softer 7 -rays proceeded much 
slower. After the CGRO mission (till ~ 2005), at soft 7 -rays above 
20 keV up to ~ 30 MeV the pulsed signals from only four pulsars 
had been detected, hampered mainly by the intrinsic source weak¬ 
nesses in this energy band with the relatively poorer sensitivities of 
the operational instruments. These were the Crab (PSRB0531-1-21) 
and the Vela (PSRB0833-45) pulsars, detected up to GeV ener¬ 
gies, PSRB1509-58 in MSH 15-52, detected up to ~ 10 MeV, and 
the twin of the Crab pulsar in the Large Magellanic Cloud (LMC) 
PSRB0540-69 detected up to ~ 50 keV (see for references section 
5). 

Significant progress could be achieved in the hard X-ray/soft 
7 -ray band when the high-resolution imaging {Chandra) and large 
sensitive area X-ray observatories {RXTE, XMM-Newton, Suzaku 
and recently NuSTAR) became available, allowing the detections 
of weak energetic point-sources at soft/medium X-rays (0.1-10 
keV) in young supernova remnants (SNRs), often discovered at 
radio-frequencies, or in location-error boxes of unidentified CGRO 
EGRET/Terrai LAT (> 100 MeV), INTEGRAL ISGRI (20-300 
keV) or H.E.S.S.A/ERITAS/MAGIC (> 30/100 GeV) sources. 
In the X-ray band below 10 keV, the instruments aboard 
these observatories were sufficiently sensitive to allow for the 
measurements of pulsar ephemerides, when there were no radio 
ephemerides available (many are radio-weak or even radio-quiet). 
For energies above 20 keV, the operational instruments were rel¬ 
atively insensitive, and could not detect independently the pulsed 
signals. But, once these timing solutions were determined at lower 
energies, phase folding of the event arrival times measured with 
RXTE PCA & HEXTE and INTEGRAL ISGRI in the hard X-ray 
band could be performed to search for and reveal the pulsed signals 
above 20 keV. 

Currently, there are 18 pulsars for which pulsed emission has 
been detected in the soft 7 -ray band above 20 keV up to ~ 50—150 
keV. In this paper we present the catalogue of soft 7 -ray pulsars, 
and for each entry an overview with the latest observational results. 
For 14 of them we present in this work new results, mainly from 
analyses of archival RXTE PCA & HEXTE and INTEGRAL ISGRI 
data for energies above 20 keV, but also from analyses of XMM- 
Newton, Chandra and AS'CA data below 10 keV, or CGRO COMP- 
TEL and BATSE data at MeV energies. For one pulsar and its neb¬ 
ula, the INTEGRAL source IGR118490-0000/PSR11849-0001 we 
present the results from our follow-up Chandra observations of this 
source together with new results from analyses of archival XMM- 
Newton, RXTE and INTEGRAL data. 

The structure of the paper is as follows. Section 2 describes the 
FermiLAT high-energy 7 -ray pulsar population of which the hard 
non-thermal power-law spectra are expected to extrapolate into the 
soft gamma-ray band. We will therefore compare the average char¬ 
acteristics of the pulsars in the two catalogues in the discussion. 
Section 3 presents briefly the properties of the instruments used in 
this work, and section 4 briefly the timing-analysis methods. In sec¬ 
tion 5, a status review is given for each of the 18 pulsars with new 
results for 14 of them. Section 6 presents the observational status 
for very young (< 15 kyr) and energetic {Lsd ^ 10^® erg/cm^ s) 
pulsars so far not detected at hard X-rays, but that are good candi¬ 
dates for future detections of pulsed soft 7 -rays, given their char¬ 
acteristics at other wavelengths. Finally, in section 7 the results are 


summarized and discussed, including the above-mentioned com¬ 
parison between the Eermi-LlXI and soft 7 -ray pulsar populations. 


2 THE Fermi-LAT HIGH-ENERGY (> 100 MEV) 
GAMMA-RAY PULSAR POPULATION 


The launch of Fermi heralded a new era in high-energy gamma- 
ray pulsar researc h. Analyzing the first three years of observations, 
lAbdo et al.l ( 1201 3h present the second Fermi LAT pulsar catalogue, 
listing 117 gamma-ray pulsars (an additional 28 were discovered 
by the LAT collaboration and other groups during the preparation 
of the catalogue). The 117 pulsars are categorized in three groups: 
radio-loud pulsars (number 42, obtained through pulse-phase fold¬ 
ing using radio-ephemerides), radio-quiet pulsars (number 35, dis¬ 
covered at 7 -energies through blind searches) and millisecond (i.e. 
recycled) pulsars (number 40). All have spin-down powers Lad 
above 3 x 10®® erg/s, which seems to be a lower-limit for detectabil¬ 
ity as 7 -ray (> 100 MeV) pulsar. 

The recycled millisecond 7 -ray pulsars are not addressed in 
this work, although some of the most energetic (Lad > 1 x 10 ®® 
erg/s) ones like PSR J0218-I-4232, PSR B1937-I-21 and PSR B1821- 


24, have been detected earlier as h ard X-ray pulsars with pu lsed 
emission detected up to 20-30 keV dKuiner et al.ll2003'3 . l2004h . 

Of the 77 { = 42 - 1 - 35) young/middle-aged (tc < 4 x 10® 
y) non-recycled pulsars, 58 (~ 75%) show at high-energy 7 -rays 
pulse profiles with two strong, well separated, caustic peaks, with 
often significant bridge emission. The radio phase lag S between the 
7 -ray leading pulse and the radio main pulse (often fiducial point) 
and the 7 -ray peak separation A show a strong anti-correlation, 
being a predicted general property of out er-magnetosphere mod¬ 
els w ith caustic pulses (see e.g. Fig. 3 of iRomani & Yadigaroglul 
Il995h . On average the 7 -ray peak separation is in the 0.4-0.5 range 
(about half a rotation), while the radio lag distribution reaches its 
maximum in the 0.15 - 0.25 range. 

In the spectral domain the pulsed emission of the 
young/middle-aged 7 -ray pulsars can well be represented by a 
power law model in combination with a simple exponential cutoff 
, i.e. F-f = k ■ {E-f/Eo)^ ■ exp( — {E-f/Ec)^) with /3 = 1 and E-/, 
the photon energy, Eq, the pivot energy minimizing the co-variance 
between the parameters, L, the photon index, Ec, the cutoff energy 
and k the normalization. 

Outer-gap and slot-gap models predict the production of 7 - 
rays at high altitudes in the magnetosphere with such a simple ex¬ 
ponential cut-off (P ~ 1 ), while in polar cap models the production 
occurs near the magnetic poles with a sharp hyper-exponential turn 
over (P ~ 2) of the spectrum at high energies. The Fermi spec¬ 
tral results thus disfavours the polar-cap models (see e.g. the de¬ 
tailed d iscussions of the fits to the Vel a pulsar spectrum 1 Abdo et al.l 
l2009ch and the Crab pulsar spectrum dAbdo et al.ll 2010 dy i 

The cutoff energies E^ are on average in the 1-3 GeV range, 
whil e the photon index distribution peaks in the -1.8 - -1.3 range. 

lAbdo et al.l ( l2013h also considered the properties of the LAT 
detected pulsars in the soft X-ray band (0.3-10 keV). Though the 
X-ray coverage is rather uneven, with the newly discovered LAT 
pulsars having often only snapshots observations with Swift XRT, 
while deep observations with Chandra and/or XMM-Newton exist 
for the well-know gamma-ray pulsars, all LAT pulsars do have X- 
ray coverage. Of the 77 young/middle-aged LAT pulsars 49 (30 
radio-loud; 19 radio-quiet) have an X-ray counterpart with a non- 
thermal (i.e. flagged with a power law) spectral component ob¬ 
tained through model fitting in terms of a power-law in combination 
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with one (or even two) blackbody(ies). The X-ray counterpart as¬ 
sociation was established either through the detection of X-ray pul¬ 
sations or positional coincidence with accurately determined radio- 
or LAT timing positions. 

For only 2 10 of the young/middle-aged LAT pulsars the pulsed 
fingerprints are detected in the canonical X-ray band (< 10 keV). 


3 INSTRUMENTS 

In this section we discuss briefly the general characteristics of the 
instruments aboard RXTE, INTEGRAL, XMM-Newton and Chan¬ 
dra, that we extensively used in this work to obtain new results 
on the (hard) X-ray/soft 7 -ray properties of the pulsar sample 
presented in this paper. For the Vela pulsar (PSRB0833-45) and 
PSRB1509-58 we obtained new pulse profiles analysing archival 
CGRO COMPTEL 0.75-30 MeV and CGRO BATSE 317-1102 
keV data, respectively, from observations covering the full CGRO 
mission lifetime. Eor PSR J2229-I-6114 archival ASCA GIS data 
were analysed to obtain new pulse profiles and a new spectrum for 
the pulsed emission. Eor the latter instruments references are given 
in those sections where data from these instruments are analysed 
by us. 


3.1 Rossi X-ray Timing Explorer {RXTE) 

In this work extensive use has been made of data from observations 
of rotation powered pulsars with the two non-imaging X-ray in¬ 
struments aboard RXTE, the Proportional Counter Array (PCA; 2- 
60 keV) and the High Energy X-ray Timing Experiment (HEXTE; 
15-250 keV). 


3.1.1 RXTE PCA 

The PCA jjahoda et al.lll99^ consisted of five collimated Xenon 
proportional counter units (PCUs) with a total effective area of 
- 6500 cm^ over a ~ 1° (FWHM) field of view. Each PCU 
had a front Propane anti-coincidence layer and three Xenon lay¬ 
ers which provided the basic scientific data, and was sensitive to 
photons with energies in the range 2-60 keV. The energy resolu¬ 
tion was about 18% at 6 keV. All PCA data used in this study have 
been collected from observations in GoodXenon mode allowing 
high-time resolution (0.9/is) analyses in 256 spectral bins. Since 
the launch of RXTE on Dec. 30, 1995 the PCA has experienced high 
voltage breakdowns for all constituting PCUs at irregular times. To 
reduce further breakdowns during its mission, which ended on Jan. 
5, 2012, not all PCUs were simultaneously operating. The most 
stable PCU was PCU-2, which was on for almost all of the time. 
On average one (50%) or two (40%) PCUs was/were operational 
during a typical observation. 


^ Table 15 of lAbdo et alj ItoTjI) labeled 13 LAT pulsars with “P” (Pulsed 
X-ray emission detect ed), omitting by mistake th is label for PSR B1509- 
58 (PSR J1513-590 8. ISeward & Harndeij h982h i and PSR J2021+3651 
iHessels et alj|2004h . while pulse d X-ray emission ha s been detected for 
6 LAT pu l sars; P SR J1459-6053 jPancrazi et ai]|2012ll. Pp J2021+4026 
iLinetalJ l2013h. PSR J 1741 -2054 iMarelli et all l20lil. PSR J1813- 
1246 iMmelli et alj2014lj) . PSR J1 836+5925 jArumugasamv et alj2014bl : 
iLin et alj2014h and PSR B0540-69 iMartin et alj2014ll after the publication 
of the catalogue. 


3.1.2 RXTE HEXTE 

The HEXTE instrument teothschild et alJ[T998h consisted of two 
independent detector clusters A& B, each containing four Na(Tl)/ 
CsI(Na) scintillation detectors. The HEXTE detectors were me¬ 
chanically collimated to a 1° (FWHM) field of view and cov¬ 
ered the 15-250 keV energy range with an energy resolution of ~ 
15% at 60 keV. The collecting area was 1400 cm^ taking into ac¬ 
count the loss of the spectral capabilities of one of the detectors. 
The maximum time resolution of the tagged events was 7.6ps. In 
its default operation mode the field of view of each cluster was 
switched on and off source to provide instantaneous background 
measurements. However, also HEXTE suffered from aging, and 
since July 13, 2006 HEXTE cluster-A operated in staring mode at 
an on-source position, while on March 29, 2010 cluster-B was com¬ 
manded to stare at an off-source position. Due to the co-alignment 
of HEXTE and the PCA, celestial targets have been observed si¬ 
multaneously. 


3.2 INTEGRAL 


The INTEGRAL spacecraft jWinkler et alJl20()^ . launched 17 Oc¬ 
tober 2002 , carries two main 7 -ray instrum ents: a high-angular- 
resolution imager IBIS (|Ubgrtini_etaLl 2003 ) and a high-energy- 
resolution spectrometer SPI jVedrenne et al.ll2003l) . The payload is 
further equiped with an X-ray monitor ing instrument, th e Joint Eu¬ 
ropean Monitor for X-rays (JEM-X; iLund et alj|2003h providing 
complementary observations in the X-ray band. These three high- 
energy instruments make use of coded aperture masks enabling im¬ 
age reconstruction in the hard X-ray/soft 7 -ray band. 

In this work, guided by sensitivity considerations, we only 
used data recorded by the INTEGRAL Soft Gamma-Ray Imager IS- 
GRI ( ILebrun et al.ll200^ . the upper detector system of IBIS, sen¬ 
sitive to photons with energies in the range ~15 keV - 1 MeV 
(effectively about 300 keV). 

With an angular resolution of about 12' and a source location 
accuracy of better than 1' (for a > 10a source) ISGRI is able to 
locate and separate high-energy sources in crowded fields within 
its 19° X 19° field of view (50% partially coded) with an unprece¬ 
dented sensitivity (~ 960 cm^ at 50 keV). Its energy resolution of 
about 7% at 100 keV is amply sufficient to determine the (contin¬ 
uum) spectral properties of hard X-ray sources in the ~ 20 - 300 
keV energy band. 

The timing accuracy of the ISGRI time stamps recorded on 
board is about 61^s. The time alignment between INTEGRAL 
and RXTE is better than 50/rs, verified using data from si¬ 
multaneous RXTE and INTEGRAL observations of the accretion- 
powe red millisecond puls ars IGR J00291-I-5934 and IGR J17511- 
3057 (lEalangaetalj200ll201ll. respectively). 


3.3 XMM-Newton 

XMM-Newton launched Dec. 10, 1999, carries three different 
Wolter type-1 multi-shell grazing incidence mirrors focusing X- 
rays to different X-ray instruments. One of the telescopes focusses 
the X-rays directly on an EPIC pn type CCD camera, while the 
other two telescopes have reflection grating arrays (RGA) in their 
light paths dispersing by reflection about 40% of the incoming 
X-ray radiation to a linear strip of CCDs (RGS detector) at their 
secondary foci, while about 44% of the incoming X-ray radia¬ 
tion passes to the prime foci onto MOS type CCD cameras. With 
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the EPIC cameras very sensitive imaging observations can be per¬ 
formed over a 30' (diameter) field of view across the 0.15-12 keV 
band with a spectral resolution (E/ AE) of 20-50 and angular res¬ 
olution of ~ 6 " 6 (FWHM). 

In th is work we only use d data collected by the EPIC pn in¬ 
strument dStriider et al.ll200lh . consisting of 12 back-illuminated 
CCDs, because the operation modes of this CCD camera are well 
suited for accurate studies of timing phenomena occuring at milli 
second scales. In particular, the Small window mode offers a fast 
read-out/frame time of about 5.67 ms, amply sufficient for most of 
the soft 7 -pulsars studied in this work, while imaging information 
is kept only for the central CCD number 4 in a frame of 63 x 64 
pixels of 4'.'l X 4'.'l size each offering thus a 4(3 x 4(37 field of 
view. 

For one pulsar in our sample, PSR10537-6910, we ana¬ 
lyzed the EPIC pn data collected in timing mode, in which the 
time stamps are registered at a ~ 30/rs resolution, while one¬ 
dimensional imaging information is conserved. The other pul¬ 
sars for which we analyzed XMM-Newton EPIC pn data were 
PSRJ0205-H6449, PSRB0833-45, PSRJ1813-1246, PSRJ1813- 
1749 and IGR118490-0000, and in these cases the data were taken 
in Small window mode. 

3.4 Chandra 

The Chandra X-ray Observatory CXO dWeisskopf et al.l2000ll was 
launched on luly 23, 1999, and combines a sub-arcsecond reso¬ 
lution X-ray telescope with advanced imaging and spectroscopic 
instruments in the focal plane of the High Resolution Mirror As¬ 
sembly (HRMA), which is composed of a nested set of four Wolter 
type-1 grazing-incidence X-ray mirror pairs with an outer diameter 
of 1.2 m and focal length of 10 m. Within the optical path (between 
the HRMA and focal plane instruments) two different types of grat¬ 
ing assemblies, the High-Energy Transmission Grating (HETG) - 
comprised of a High Energy Grating (HEG) and Medium Energy 
Grating (MEG) on a single structure - and the Low Energy Trans¬ 
mission Grating (LETG), can be placed mechanically on command. 

The Science Instrument Module (SIM) at the focal plane 
houses two instruments, the Advanced CCD Imaging Spectrome¬ 
ter (ACIS) and High Resol ution Camera (HRC ). 

The ACIS instrument dCarmire et al.l20()^ contains 10 planar, 
1024 X 1024 pixels (pixel scale is 0'.'4920 ± O'.'OOOl) CCDs of 
which 4 are arranged in a 2x2 array (16(9 x 16(9 field of view) for 
imaging purposes (ACIS-I) and the other 6 in a 1 x 6 array (8(3 x 
50(6 field of view) serving either as grating readout or as imaging 
detector (ACIS-S; containing two back-illuminated CCDs). 

The energy resolving power (E/AE) is moderate and de¬ 
pends on the CCD type (front or back-illuminated), is energy and 
time (Charge Transfer Inefficiency) dependent, and is in the range 
~ 10-60 for the FI CCDs and ~4-45 for the BI CCDs. The BI type 
CCDs (~ 0.18-10 keV) have higher effective area at lower energies 
than the FI type CCDs (^0.48-10 keV). The effective area reaches 
a maximum of about 600 cm^ at 1.4 keV for the FI CCDs. The 
nominal frame time in Timed Exposure (TE) mode is 3.2 s with a 
40.96 ms read-out time. This relatively long integration time can 
cause pile-up (more than one event registered in a pixel during a 
frame time) effects even for moderate count rates. 

Although in TE mode other frame times are possible using 
subarray read-outs, the default TE mode is not suitable for pulsar 
timing studies at millisecond time scales, but it is very useful to 
obtain sub-arcsecond (energy resolved) images of the pulsar en¬ 
vironment. The other operation mode of the ACIS is Continuous 


Clocking (CC). In this mode the integration time is 2.85 ms, and 
this is reached at the expense of spatial information in one direc¬ 
tion. _ 

_ The other focal plane instrument is the HRC jMurrav et al.l 

l 200 (]h . which is comprised of two multichannel plate imaging de¬ 
tectors: the HRC-I offering wide-field imaging (30' x 30' field of 
view with O'.'13175/pixel), and the HRC-S serving mainly as read¬ 
out for the LETG. The HRC detectors do not/hardly have energy 
resolving power. Due to a wiring problem the time resolution of 
the HRC-I/S detector is governed by the actual count rate and is 
much worse than the anticipated resolution of about 16 ps. How¬ 
ever, the central segment of the HRC-S ( 6 ' x 30'field of view) can 
be operated in a special mode, designed for accurate timing studies, 
in which the original 16 ps time resolution can be recovered. 

In this work we present (Pl-obtained) Chandra data for 
IGR 118490-0000 from two 25 ks observations, one with the HRC- 
S in timing mode and one with the ACIS-S. Moreover, we have esti¬ 
mated the genuine pulsed fraction of PSR 11617-5055 using Chan¬ 
dra HRC-S timing mode data. 


4 ANALYSIS METHODS 

In the analysis of the high-energy data from the non-imaging RXTE 
instruments PCA and HEXTE, the coded-mask imager ISGRI 
aboard INTEGRAL and the X-ray telescopes aboard Chandra and 
XMM-Newton some common analysis procedures exist. To obtain 
pulsed flux measurements accurate timing models (ephemerides) 
are required describing every revolution of the neutron star ac¬ 
curately. These so-called phase coherent ephemerides are mainly 
based, in this work, on Time-of-Arrival (ToA) determinations of a 
pulse template using RXTE PCA data, and are listed in Table[T]for a 
limited set of soft 7 -ray pulsars, for which RXTE PCA monitoring 
observations exist. The method to o btain accurate timing mode ls is 
extensively described in Sect. 4.1 of iKuiper & HermserJ ( l2009h . 

Irrespective of the high-energy instrument with which the 
data have been collected a typical timing analysis consists of 
1) data selection using e.g. the HE AS ARC browse interface 
for RXTE, XMM-Newton and Chandra observations fulfilling the 
source search criteria or the browse facility at the ISDC for INTE¬ 
GRAL data, 2) screening the data by creating good-time intervals 
(GTI), 3) barycentering the event time tags using the instantaneous 
spacecraft ephemeris (position and velocity) information, the JPL 
solar system ephemeris information (DE200/DE405) and an ac¬ 
curate source position to convert the times from Terrestial Time 
scale (TT or TDT, which differs from Coordinated Universal Time 
(UTC) by a number of leap seconds and a fixed offset of 32.184 s) 
into Barycentric Dynamical Time (TDB) scale, a time standard for 
Solar system ephemerides. 

Subsequently, the TDB ti me tags are converte d to pulse phases 
(see formula ( 1 ) in Sect. 3.4 of iKuiper et al.l20T3) adopting appro¬ 
priate timing models (phase folding). Next, the events are sorted 
on pulse phase and energy to obtain 2 d event distributions, storing 
thus pulse profiles for any energy band covering the passband of the 
involved instrument. The detection of a pulsed sig nal is expressed 
in a nu mber of ct’s applying the Z'^ test described in iBuccheri et al.l 
(Il983h . Typical timing analyses of RXTE PCA, HEXTE, XMM- 
Newton and INTEGRAL ISGRI data are described in Sect. 3.1, 
Sect. 3.7 and Sect. 3.8 of I Kuiper et al.l ( 1201 Oh and Sect. 4.3 of 
[Kuiper & Hermsetil ( l2009tl . respectively. A ty pical timing analysis 
for Chandra HRC data is given in Sect. 4 of iKuiper et alj ( l2002h 
and Sect. l5.15] of this work. 
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From the pulse-phase distributions per energy band we derived 
the number of pulsed excess counts by fitting with (a) pulse profile 
template(s) or a truncated Fourier series. The pulsed fluxes were 
obtained by applying forward folding spectral fitting procedures 
taking into account the instrument spectral response, effective area, 
exposure time and absorption in t he interstellar medium. Th is is 
explained for example in Sect. 5 of iKuiper & Flermsenl ( l2009l) for 
P CA, HEXTE and IS GRI data. For XMM-Newton (see e.g. Sect. 4 
of iKuiper et alj201^ and Chandra equivalent procedures exist. 


5 THE SOFT y-RAY PULSAR POPULATION 
5.1 PSRJ0205-I-6449 

PSR J0205-I-6449 is a young (r ~ 5.4 kyr) 65 ms pulsar located 
near the center of super nova remnant/pulsar wind nebula 3C58 
with a spin-down luminosity of 2.7 x 10^^ erg/s. Its pulsation was 
discovered in the X-ray band using Chand ra data and confirmed 
using RXTE PCA data jMurrav et alj|2002h . and its weak pulsed 
signal was later detected at radio frequencies. 

At TeV-energies (> 300 GeV) VERITAS did not (yet) de¬ 
tect emission from the sky region including 3C58 , yielding a 99% 
flux u pper-limit of about 2% of the Crab nebula dHumenskv et aU 
I 2 OO 9 I) . However, the upgraded stereoscopic MAGIC telescope de¬ 
tected 3C58 as unresolved source at a 5.7cr level with an integra l 
flux above 1 TeV of about 0.65% of the Crab jAleksic et al.l20149) . 

Detailed timing st udies in the X-ray and radi o-bands have 
been p erformed by both lUivingstone et al.l (l2009h and iKuiper et al.l 
( 120101) . The X-ray pulse profile shows two sharp peaks separated 
in phase by 0.488 ± 0.002 with the main peak lagging the sin¬ 
gle radio pulse by 0.089 in phase. The pulsed X-ray/soft 7 -ray 
spectrum over the energy band 0.56-267.5 keV (based on XMM- 
Newton EPIC-pn, RXTE PCA and HEXTE measurements) for the 
sum of the two puls es is very hard with photon index —1.03 ±0.02 
dKuiper et al.l201^ . No difference in spectral shape is found for the 
two pulses in the X-ray/soft 7 -ray band. High-energy 7 -ray pul¬ 
sations (> 100 MeV) from PSR J0205-I-6449 have been detected 
using Eermi LAT data dAbdo et alj|2009ah . The 7 -ray lightcurve 
(> 100 MeV) shows also two pulses, aligned with the X-ray pulses, 
but now the strongest 7 -ray pulse coincides with the weakest X- 
ray pulse, indicating different spectral behaviours of the two pulses 
when considering the full 0.5 keV to 10 GeV band. Eor more 
detailed information about t he high-energy char acteri stics of the 
pulsed emission we refer to iKuiper et all ) l 2010 l) and lAbdo et ^ 
(l2009td) . 


An intriguing feature of its pulsed emission is that the double- 
peaked pulse-profile with a peak separation of about 0.4 (in phase) 
keeps its identity and preserves its alignment approximately (the 
X-ray/ 7 -ray main pulse preceeds the radio main pulse by about 
300/rs) across the electro-magnetic spectrum (see e.g. Fig. 6 of 
IKuiper et alj2003l) . though the relative intensities of the two pulses 
vary. 


The total emission of the Crab (nebula and pulsar) has long 
time served as flux/spectral calibration target for various instru¬ 
ments because of its detectability across the ‘full’ electro-magnetic 
spectrum and its stability over time. The latter property has been 
abandoned recently by the dis covery of considerable var iations in 
the hard X-ray/soft 7 -ray band dWilson-Hodee et alj 2011 [) and flar¬ 
ing a ctivity at high-energy -v-ravs l lTavani et ^ 2011 ; Abdo et ^ 
l201lh . Contrarily, the pulsed emission of PSR B0531-1-21 as seen by 
RXTEVCN (3.2-35 keV) seems rathe r stable and decreases stead ily 
consistent with the pulsar spin-down dWilson-Hodee et^l201 ll) . 

The high-energy spectrum of the pulsed emission of 
PSRB0531-I-21 across 8 decades in energy from soft X-rays up 
to high-energy 7 -rays is shown in the multi-source spectral com¬ 
pilation depicted in Fig.as red-coloured datapoints. Data from 


CGRO OSSE, COMPTEL and EGRET, BeppoSAX LEGS, MECS 
and PDS and the ball oon-borne GRIS instruments are included 
(see IKuiper et alifTOOlL for more details), supplemented by RXTE 
PCA (5-50 keV, E_2 50us_12 8M_14_ls event mode; this work) 
and HEXTE (~ 13-235 keV; this work) pulsed flux measurements. 
The HEXTE data have been collected during observation run 40805 
covering the time periods March 17-31, 1999 and December 18- 
19, 1999, for an effective deadtime corrected ON source time of 
22.7 ks and 23.8 ks for cluster A and B, respectively. Also included 
in Fig.l^is the model s pectrum of the pulse d emission (100 MeV - 
10 GeV) as obtained bv lAbdo et al.l ) l2010cl) using Fermi LAT data. 
Maximum pulsed flux is reached near 100 keV, while beyond 30 
MeV a flattening occurs, followed beyond ~ 10 GeV by a power- 
law like d ecay, detected up to about 400 GeV by the VERITAS TeV 
telescope dAliu et al.ll201 ih . 

The (unabsorbed; using a Nh of 3.61 x 10^*^ cm~^) 2-10 
keV, 20-100 keV 1-10 MeV and 0.1-1 GeV pulsed fluxes for the 
Crab pulsar are (1.948 ± 0.005) • 10"®, (2.766 ± 0.008) ■ 10“®, 
(1.36 ± 0.05) • 10"® and (7.87 ± 0.28) • 10"^° erg/cm®s, respec¬ 
tively, derived from RXTE PCA/HEXTE data (2-10 and 20-100 
keV, this work), CGRO COMPTEL (see e.g. IKuiper et al.ll200lL 
from the pulsed flux measurements reported in that work we de¬ 
rived the 1-10 MeV pulsed flux value) and Fermi LAT (flux value 
derived from the spect ral parameters and their la statistical errors 
as given in Sect. 4.3 of lAbdo et al.ll 2010 ch . 


5.3 PSRJ0537-6910 


5.2 The Crab pulsar, PSR B0531+21 / PSR J0534-(-2200 

The 33 ms pulsar, PSRB0531-H21 or PSR J0534-H2200, in the 
Crab ne bula (3C 144; NGC 1952; CTB 18) was discovered in 
1968 bv IStaelin & ReifensteinI l ll968l) . Since its detection at radio¬ 
frequencies the pulsed fingerprint of PSR B0531-1-21 has been 
detected across the full electro-magnetic spectrum from radio, 
optical/infra-red, soft X-rays (0.1-2 keV), hard X-rays (2-20 keV), 
soft 7 -rays (2 0-500 keV), mediu m-energy 7 -rays (0.5 MeV - 30 
MeV, see e.g. IKuiper et alj|200ll for a high-energy picture from 
soft X-rays up to hig h-energy 7 -rays), h igh-energy 7 -rays (30 MeV 
- 100 GeV, see e.g. lAbdo et ^ 2010cl) and recently even at TeV 
energies (> 0.1 T eV, see e.g. Aliu et alJl2008l ; lAiekkc et al.ll 201 iL 
|2014a| ; |Aliu et al.ll201 ll) . 


iMarshall et Zl ) ll998all9) reported the discovery of a young (r ~ 
4.9 kyr) fast 16 ms pulsar in LMC supernova remnant N157B 
(NGC 2060) using RXTE PCA data taken from 2 observation runs 
in 1996 centered on the nearby SN 1987A and archival ASCA 
GIS observations to pin down its location and derive its spin- 
down. So far, PSR J0537-6910 is the fastest and most energetic 
(E ~ 4.9 X 10^® erg/s) non-recycled spin-down powered pulsar. 
The pul sed signal was also found in archival BeppoSAX MECS (2- 
10 k eV; ICusumano et alJll998all9) and ROSAT HRI data (0.1-2.4 

and in follow-up Chandra HRC-S 


keV; 


Wang & Gotthell 


data I Wang et al.ll200lh . 


On January 19, 1999 an intensive RXTE monitoring (observa¬ 
tion id. 40139) of PSR J0537-6910 commenced, which ended with 
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Table 1. Phase-coherent ephemerides, derived in this work, for the subset of soft 7 -ray pulsars for which RXTE PCA monitoring data exist. 


Pulsar 

# 

Start 

[MJD] 

End 

[MJD] 

to, Epoch 
[MJD,TDB] 

V 

[Hz] 

V 

xl0~“ Hz/s 

V 

X 10-22 Hz/s2 

Validity range 
(days) 

PSRJ0537-6910 

51197 

51263 

51237.0 

62.0402460583(64) 

-19.9221(2) 

99(40) 

67 

PSRJ0537-6910 

51294 

51423 

51366.0 

62.0380675984(19) 

-19.92497(3) 

84(5) 

130 

PSRJ0537-6910 

51422 

51547 

51480.0 

62.0361053238(9) 

-19.92024(2) 

46(6) 

126 

PSRJ0537-6910 

51576 

51706 

51628.0 

62.0335856000(14) 

-19.92746(3) 

72(3) 

131 

PSRJ0537-6910 

51715 

51786 

51761.0 

62.0313152172(24) 

-19.9279(1) 

248(16) 

72 

PSRJ0537-6910 

51795 

51818 

51795.0 

62.03072978(2) 

-19.920(2) 

0.0 

24 

PSRJ0537-6910 

51833 

51875 

51833.0 

62.0300843300(32) 

-19.9294(2) 

0.0 

43 

PSRJ0537-6910 

51886 

51955 

51886.0 

62.02918045(1) 

-19.9357(8) 

276(30) 

70 

PSRJ0537-6910 

51964 

51996 

51964.0 

62.0278657375(87) 

-19.9354(7) 

0.0 

33 

PSRJ0537-6910 

51997 

52145 

52049.0 

62.0264019706(14) 

-19.92818(3) 

63(2) 

149 

PSRJ0537-6910 

52132 

52166 

52132.0 

62.0249729876(73) 

-19.9153(6) 

0.0 

35 

PSRJ0537-6910 

52175 

52230 

52175.0 

62.0242443562(98) 

-19.9364(9) 

167(38) 

56 

PSRJ0537-6910t 

52252 

52368 

52252.0 

62.0229447197(59) 

-19.9342(2) 

68(5) 

117 

PSRJ0537-6910 

52389 

52446 

52389.0 

62.02059598(1) 

-19.936(1) 

134(56) 

58 

PSRJ0537-6910 

52459 

52540 

52459.0 

62.019403934(12) 

-19.9359(7) 

103(22) 

82 

PSRJ0537-6910 

52551 

52648 

52551.0 

62.0178456346(56) 

-19.9374(3) 

94(7) 

98 

PSRJ0537-6910 

52647 

52718 

52647.0 

62.016192268(14) 

-19.9327(9) 

148(29) 

72 

PSRJ0537-6910 

52745 

52792 

52745.0 

62.0145138255(43) 

-19.9360(3) 

0.0 

48 

PSRJ0537-6910 

52822 

52884 

52822.0 

62.013203356(14) 

-19.9401(9) 

141(36) 

63 

PSRJ0537-6910 

52889 

53008 

52889.0 

62.0120638182(43) 

-19.9409(2) 

132(3) 

120 

PSRJ0537-6910 

53019 

53122 

53019.0 

62.0098457965(56) 

-19.9405(3) 

95(7) 

104 

PSRJ0537-6910 

53128 

53143 

53128.0 

62.007969319(34) 

-19.939(6) 

0.0 

16 

PSRJ0537-6910 

53147 

53285 

53147.0 

62.0076662548(52) 

-19.9426(2) 

95(3) 

139 

PSRJ0537-6910 

53290 

53366 

53290.0 

62.005227527(14) 

-19.9488(8) 

240(29) 

77 

PSRJ0537-6910 

53365 

53444 

53365.0 

62.003935314(18) 

-19.9404(9) 

192(30) 

80 

PSRJ0537-6910 

53446 

53549 

53446.0 

62.002556305(10) 

-19.9471(4) 

165(10) 

104 

PSRJ0537-6910 

53551 

53682 

53609.0 

61.9997679030(28) 

-19.94013(3) 

106(5) 

132 

PSRJ0537-6910 

53710 

53860 

53792.0 

61.9966405925(17) 

-19.93914(2) 

83(3) 

151 

PSRJ0537-6910 

53862 

53947 

53862.0 

61.995449381(11) 

-19.9510(5) 

230(17) 

86 

PSRJ0537-6910 

53952 

53996 

53952.0 

61.993899826(27) 

-19.945(3) 

531(154) 

45 

PSRJ1811-1925 

52341 

52541 

52341.0 

15.4612089999(14) 

-1.056106(19) 

0.0 

201 

PSRJ1811-1925 

52595 

52904 

52595.0 

15.4609772117(10) 

-1.056411(9) 

0.0 

310 

PSRJ1811-1925 

52904 

53333 

52904.0 

15.4606951476(22) 

-1.056670(25) 

-2.03(15) 

430 

AXJ1838.0-0655 

54513 

55002 

54513.0 

14.1847581888(11) 

-0.992954(11) 

1.95(6) 

488 

AX J 1838.0-0655 

55038 

55537 

55136.0 

14.1842449934(1) 

-0.999054(6) 

3.8(1.2) 

497 

IGRJ18490-0000 

55525 

55545 

55526.0 

25.961259794(16) 

-0.961(2) 

0.0 

19 

PSRJI930-tl852 

52529 

52633 

52529.0 

7.30613850(9) 

-4.008(2) 

0.0 

104 

PSRJ2022-r3842 

55223 

55232 

55223.0 

20.58651709(2) 

-3.649(6) 

0.0 

10 


fThis entry has been used for XMM-Newton EPIC-pn observation 0113020201 in timing mode 


cessation of the RXTE mission in early January 2012, aimed at fol¬ 
lowing its timing behaviour with time. Analyzing the first 2.6 years 
of this monitoring campaign re sulted in a glitch rate of 2.3 per year, 
the highest seen for any pu lsar 1 Marshall_etaLy00J). An update on 
this work was provided bv iMiddleditch et alj ( I2OO6I1 analyzing 7.6 
years of RXTE monitoring data, and finding 23 glitches at a rate 
of 3.3 per year. The latter authors found a correlation between the 
time interval between two glitches and the amplitude of the first 
glitch (of the interval), yielding a prediction of these intervals with 
an accuracy of a few days. 


While its timing characteristics have been studied intensively, 
the high-energy spectral characteristics of PSR J0537-6910 are 
poorly determined. Therefore, we invoked the huge RXTE data 
archive on PSR J0537-6910 to establish its high-energy spectral 
char acteristics. For this purpo se we used a (slightly larger) dataset 
than iMiddleditch et al.l fcoOfil) and revisited their analyses to obtain 
the alignment for each of the 24 timing segments enabling us to co¬ 


add the pulse profiles for each of the 24 timing segments directly 
(Note : this alignment information is not given in lMiddleditch et al.l 

l2006lfl 

We started with the derivation of ToA’s for all observations 
between Jan. 19, 1999 and Oct. 5, 2006 (time span MJD 51197 - 
53996; 7.7 years). From these ToA’s we determined in a consistent 
way phase coherent ti ming models, which are in good agreement 
with those obtained by iMiddleditch et al.l l l2006h and are listed in 
Table[T] Across this period we found 24 glitc hes of which 23 have 
been reported on bv iMiddleditch et^ ( l2006l ). resulting in 24 tim¬ 
ing segments for which phase coherent solutions (all with proper 
alignment) have been derived. 


^ Analyzing the full RXTE database on PSR J0537-6910 is outside the 
scope of this paper, but work on this, in particular to derive the charac¬ 
teristics for all the glitches across the full RXTE data period, is in progress. 
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Figure 1. PSR J0537-6910; pulse profiles for XMM-Newton EPIC pn (0.7- 
2 keV; panel a), RXTE PCA (1.7-27.3 keV; panel b) and RXTE HEXTE 
(15.6-50.3 keV; panel c). Pulsed emission has been significantly detected 
up to ~ 50 keV. 


Subsequent pulse-phase folding of PCA events (from all de¬ 
tector layers), screened according to standard criteria and removing 
break-down periods of the detector units, yielded high-statistics 2 
dimensional event distributions sorted on pulse-phase and energy 
(PHA). In Fig.[T]panel b the 120 bin 1.7-27.3 keV profile is shown, 
analyzing PCA data from time segments 1-18, which provides us 
already with very high pulsed-signal statistics. Similar procedures 
for HEXTE yielded, now analyzing all 24 segments given the poor 
pulsed-signal statistics contrary to the PCA case, a lOcr pulsed 
signal for the 15.6-50.3 keV band (see Fig.[T]panel c), and even a 
5.2 (t signal for the 31.0-50.3 keV band, establishing PSRJ0537- 
6910 as a soft 7 -ray pulsar. Indications for pulsed emission were 
found for energies above 50 keV. 

In order to determine the timing characteristics of PSR J0537- 
6910 below 2 keV we also analyzed EPIC pn timing mode data 
from a 35.9 ks XMM-Newton observation (obs. id. 0113020201) 
performed on November 19, 2001 (data period MJD 52232.949 
- 52233.366). Default screening and optimizing on the RAWX 
event parameter by allowing only events with RAWX in the range 
[26, 42] yielded highly significant pulse profiles after phase fold¬ 
ing using the PSR10537-6910 ephemeris entry with epoch MJD 
52252 (see Table[TE. The resulting 0.7-2 keV XMM-Newton EPIC 


® iMiddleditch et alj i200d) reported a glitch (#8 in their table 4) occuring at 



10 "^ 10 "^ 10 “' 
E [MeV] 


Figure 2. PSR J0S37-6910; high-energy (0.7-250 keV) unabsorbed pulsed 
spectrum as derived from measurements by XMM-Newton EPIC-pn (0.7-10 
keV; red), RXTE PCA (2.5-30 keV; aqua) and RXTE HEXTE (15-250 keV; 
blue). The best fit model (“curved” power-law) to all these pulsed flux mea¬ 
surements is superposed as solid black line. The red dotted line represents 
the best power-law fit model to solely EPIC-pn flux measurements. 


pn pulse profile is shown in Fig[T] panel a, and the deviation from 
uniformity is about 16.4 (t applying a Zftest. 

In order to derive the pulsed spectrum of PSR 10537-6910 we 
fitted a model composed of a (flat) background component and a 
component describing the shape of the single pulse, as derived from 
the high-statistics RXTE PCA profile, to the pulse profiles of any 
of the energy bands in the band width of XMM-Newton EPIC pn, 
RXTE PCA and HEXTE. The resulting pulsed excess counts are 
converted to pulsed fluxes in a forward folding procedure assuming 
an underlying (photon) power-law model and adopting interstellar 
abso rption under a column of Nh of 7 x 10^^ cm“^ (see e.g. Table 
4 of iTownslev et al.ll2006L for X-ray spectroscopic information of 
some sources in the 30 Doradus region, including PSR 10537-6910 
and SNR N157B, using Chandra ACIS data). 

In Fig. 12 the pulsed flux measurements across the 0.7-250 
keV band are shown for XMM-Newton EPIC pn (red data points), 
RXTE PCA (aqua) and RXTE HEXTE (blue) along with fits to 
solely EPIC-pn data (power-law model; red dotted line; 0.7-10 
keV) and to all data (“curved” power-law; black solid line; 0.7- 
250 keV). The later model provides a 4.3cr improvement over a 
simple power-law model fit considering the full dataset, and thus 
proves the necessity of a gradually breaking spectrum. The unab¬ 
sorbed pulsed flux in the 2-10 keV band for the “curved” PL-model 
is (7.42±0.26) x 10“^® erg/cm^s, about 10% higher than the 2-10 
keV flux of (6.72 ± 0.38) x 10“^® erg/cm^s derived from the PL- 
fit using solely EPIC-pn data (the photon index was —1.57 ± 0.03 
in this case). Both values are, however, consistent with the value of 


MJD 52241.6 it 7.8, however, the XMM-Newton data can still be correctly 
described by the post-glitch ephemeris entry and not with the pre-glitch 
ephemeris, indicating that the glitch must have occured somewhere between 
MJD 52230 (the end time of the pre-glitch ephemeris) and MJD 52232.949, 
the start time of the XMM-Newton observation. 
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Figure 3. PSRB0540-69; pulse profile (broad single-pulse with a dip in the 
middle) for the 2-20 keV band as measured by RXT E PCA (upper panel ) 
combining about 700 ks of effective exposure time (see lde Plaa et alj2003l) . 
In the lower panel the RXTE HEXTE lightcurve of PSR B0540-69 is shown 
for the ~ 13-48 keV band (Zj ~ 8.Sc) using about 250 ks of dead-time 
and off-axis corrected exposure time for clusters A and B. 


(6.7 ± 0.6) X 10“^® erg/cm^s given by [Marshall et alj h998bh . 
Comparing the unabsorbed pulsed flux in the 2-8 keV band of 
(6.04 ± 0.24) X 10“^® erg/cm^s (this work) with the unabsorbed 
flux v alue of 3.09 x 10~^^ erg/cm^s as given by iTownslev et al.l 
(l2006h for the total emission of PSR J0537-6910 indicates that 
about 18% of the emission is pulsed in this band, similar to what 
has been found for the Crab pulsar/nebula system. 

Noteworthy is the recent detection of PSR J0537-6910 by IN¬ 
TEGRAL ISGRI in the 20-60 keV band, at a flux of 0.44 ± 0.08 
mCrab (total emission), resulting from a very deep survey of the 
LMC region dOrebenev et al.ll20O . 

On the other hand no pulsed radio and high-energy 7 -ray 
emissions have been detected so far fr om PSR J0537-6910 (see 
ICrawford et ^ 1 19981: lAbdo et al.l 1201 3l, for the radio and high- 
energy 7 -ray band, respectively). In the TeV domain, however, 
the H.E.S.S. collaboration reported recently on the detection, 
at a 14(7 level, of a point source positionally consistent with 
N157B/PS R J0537-6910 making this the first extragalactic PWN 
association dH.E.S.S. Collaboration et ^l2012h . 


5.4 PSR B0540-69 / PSR J0540-6919 


PSRB0540-69 was discovered at X-rays by ISeward et all dl984h 
using Einstein X-ray obs ervatory data, soon follo wed by a detec¬ 
tion in the op tical band i Middleditch et al.lll985l) and a weak ra¬ 
dio detection dManchester et al.lll993l) ab^t a decade later. It is a 
young Crab-like pulsar located in the Large Magellanic cloud, and 
it is one of the few known spin-down powered pulsars resident out¬ 
side our own galaxy (see iRidlev et al.ll2013l . for an update on the 
spin-down powered pulsar content of the LMC, now 23 members, 
and SMC, 5 members). 

Since its discovery PSRB0540-69 has been a primary 
target for many d ifferent X-r ay _ instruments e.g. Ginga 


dNagase et alj Il990l : Deeter et ^ 1999ll . ROSAT dfinle^jt^ 
1993 |l . ASCA dHiravama et al.l l2002l) . Bemo SAX dMineo et al. 


Cusumano et al 


2003h . Chandr a 1 Gotthelf&Wan 


Kaarete^ iJ l200lh and RXTE dde Plaaet alJ 120031: 


Livingstone et akluOPSl: Campana et akluOPa). 

Ide Plaa et al.l d2003h reported for the first time the detection 
of pulsed emission from PSR B0540-69 in the hard X-ray band up 
to ~ 50 keV using RXTE PCA and HEXTE data. Figure [3] shows 
the pulse profiles of PSRB0540-69 for the ~ 2-20 and 13^8 keV 
bands, as measured by RXTE PCA (pa nel a) and HEXTE (p anel b), 
respectively, using the same data as in ide Plaa et al.l d2003l) . More¬ 
over, these authors found that the pulsed spec trum softens towards 
higher energies. An update of the findin gs by Ide Plaa et alj d2003h 
was presented bv ICampana et al.| d2008ll using about 3 times more 
RXTE PCA and 5 times more HEXTE exposure time. In general, 
these authors reached similar conclusions, namely that the pulsed 
spectrum softens towards higher energies, but obtained a differ¬ 
ence in flux normalization of about 30% and a somewhat differ¬ 
ent spectral shape. The first discrepancy can be traced back al¬ 
most completely to improper re sponse functions for the PCA units 
(i.e the sensitive areas) used by Ide Plaa et al.l (l2003l) at the time of 
the production of their paper, thereby overshooting the normaliza¬ 
tion by more than 20 %, while the shape discrepancy can be traced 
back by the different Hydrogen column densities used by th e dif- 
fere nt authors, i.e. 4.(3 x 10 ^^ cm~^ by Ide Plaa et alJ l l2003l) (see 
e.g. iKaaret et al.ll 200 lh . while IcWpana et alT i 2008 1 derived/used 


a slightly lower value of 3.7 x 10^^ cm ^. We p refer the higher val¬ 
ued nu mber, which finds support in the work of ISerafimovich et akl 

(I2OO4I) . 

The softening of the pulsed spectrum explains also the non¬ 
detections by CGRO BAT p > 20 keV and CGRO OSSE (in 
the 80-210 keV band) (see lWilson eraDll993bl:|Hertz et alJll9^ 
respectively). The unabsorbed pulsed fluxes of PSRB0540-69 
for the 2-10 and 20-100 keV bands are ( 6.7t.n't) ■ 10~^^ an d 
( 6 . 1^2 f) erg/cm^s, respectively (see lCampana et alj2008l) . 

The former flux value was derived fr om the absorbed flux v alue of 
(6.51 q 7 ) • erg/cm^s given bv ICampana et al.l ( I2OO8I) adopt¬ 
ing their estimated Hydrogen column density Nh of 3.7 x 10^^ 
cm“^. 


A comparison of the spectral energy distributions of 
PSRB0540-69 and its PWN for the ~ 1-100 keV band using Swift 
XRT (WT mode), RXTE PCA and HEXTE for the pulsed emis¬ 
sion, and Swift XRT (PC mode) and INTE GRAL ISGRI for the to - 
tal (=pulsar plus PWN) emission is given in ICampana et al.l (l2008l) . 
and shows that about 25% of the t otal emission is pulsed . Notewor¬ 
thy is also the work presented bv iGrebenev et alj (12013l) on a deep 
hard X-ray survey of the LMC using INTEGRAL ISGRI (4.8 Ms 
effective exposure; PSRB0540-69 (plus PWN) was detected at a 
22 .Oct significance level in the 20-60 keV band with a 20-60 keV 
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flux of 1.68 ± 0.08 mCrab) and INTEGRAL JEM-X (1.8 Ms effec¬ 
tive exposure; 3-20 keV flux of 1.28 ± 0.07 mCrab 18 ct)). 

At GeV energies recently the pulsed fingerprint of 
PSRB0540-69 has been detected at a 6.8 (t level u sing a RXTE 
based ephemeris covering a time span of -^3.5 years jMartin et al.l 
l2014h . The GeV pulse profile exhibits a single broad pulse and 
the GeV spectrum of the pulsed emission is soft and properly 
described by a cutoff power-law model connecting smoothly to the 
soft 7 -ray part of the spectrum. 

So far, (pulsed) emissio n has not been detected from 
PSRB0540-69 at TeV energies dKomin et al.l2012h . 


5.5 The Vela pulsar, PSR B0833-45 / PSR J0835-4510 


Targeting t he V ela Supernova remnant at radio-frequencies 
iLarge et al.l (Il968h discovered short period pulsations with P ~ 89 
ms from a young (11 kyr) and energetic {Lsd ~ 6.9 x 10^® 
erg/s) neutron star, known since as PSRB0833-45. Since then the 
Vela pulsar has been monitored intensively at radio-frequencies, 
and it turns out to be a frequent glitcher, e.g. nine major glitches 
occured in the time period 1969 - 1994. Detailed single-pulse 
studies showed that the radio pulse (single; see e.g. Fig. |4] panel 
a) cou ld be described as a composite of four (Gaussian) compo¬ 
nents dKrishnamohan & Downs! Il 9831) covering a phase extent of 
~ 0.055. Morphology changes of the radio pulse as a function 
of frequency (1.4-24 GHz) can be described by merely intensity 
variations of the constituting fou r Gaussian pulses, while maintain¬ 
ing their wid ths and positions dKeith et alj|201 ih . Further single 
pulse studies jjohnston et al.l200l[) reveal the presence of giant mi¬ 
cropulses in the leading edge of the radio pulse and on the trailing 
edge a large Gaussian component. 

In the optical domain puls ations from PSRB0833-45 were 
detected by [Wallace et aljj 19771) fro m a sky region including the 
candidate proposed by Laskeil ( 19761) . The statistics of the (double 
pulsed) optical profile has increased sig nificantly since the first de¬ 
tection (see e.g [Manchester et al.ll 198011 with a most ac curate char¬ 
acterisation of the profile given by Gouiffes. C. I ( [ 19981) using fast- 
photometry observations performed at the ESO 3.6 m telescope 
in 1993 and 1994. The latter observations revealed a third narrow 
pulse coinciding with the radio pulse (see e.g. Fig.|4]panel b), the 
presence of which was confirmed at smaller wavelengths by HST 
STIS NUV and FU V observations (see e.g. Fig.|4]panels c and d, 
[Romani et al.l200^ . The latter observations revealed even a narrow 

fourth pulse preceding the radio-pulse. _ 

_ Multi-band optical photometry with E SO NTT l|Nasij(i_et_^ 

[l997h . ESO VLT (ISAAC, J , and H bands [Shibanov et al.|[2003[) 
H5T STIS (NUV and FUV, [Romani et al.|[2005l; [K^galtsev et al.[ 
[2007h and Gemini GeMS (X^-band. [Zvuzin et al.[ [ 2 OI 3 I) along 
with spectroscopy observations (VLT FORS2-, 4000- 11000 A, 
[Mignani et al.[[200% yielded a detailed spectral picture from near- 
IR to FUV of the Vela puls ar revealing a very flat spectrum with 
index a = —0.01 ± 0.01 ([Zvuzin et al.[[201^ across this range 
confirming, inline with its high pulsed fraction of about 80%, the 
non-thermal magnetospheric nature of the near-IR to FUV emis¬ 
sion. 

Optical astrom etry with (initially) th e ESO 3.5 m NTT 
([Nasuti et alj [ 19971). HST WFCAVFPC2 jde Luca et all [20001; 


[Caraveo et al.[ 200 1[) made accurate proper motion and even paral¬ 
lax measurements of the Vela pulsar possible, yielding a distance of 
2941^® pc, consistent with the more accurate distance estimate of 
2871,® pc derived s omewhat later from VLBI (radio) observations 
([Podson et al.l2()o3) . 


In the soft X-ray band (0.1-4 keV) the first detailed imaging 
studies of the Vela pulsar environment have been performed with 
the IPC (-^ spatial resolution) and HRI (~ 2" resolution) instru¬ 
ments aboard the Einstein satellite (Nov. 12, 1978 - April 1981). 
Structures a t four different spatia l scales could be discerned from 
the images ( [Hamden et al.[ [ l98^ . in particular a) a pointlike ob¬ 
ject coincident with the pulsar, b) a relatively bright ~ 4' diffuse 
nebula about the pulsar, c) diffuse hard emission at ~ 1° scale be¬ 
tween the pulsar and the radio object Vela X and finally d) thermal 
emission from the entire ~ 5° Vela SNR. No evidence for X-ray 
pulsations could be found from the pointlike object associated with 
PSRB08 33-45, indicating a rather low p ulsed fraction of < 9% 
(see also. [6gelman & Zimmermann[[ 19891 using EXOSAT 0.03-2.4 
keV imaging observations). 

The first unambiguous detectio n of pulsed X- r ay em ission 
from PSRB0833-45 was reported bv [6gelman et alj ([l993[) using 
ROSAT PSPC (0.1-2.4 keV) data. This was confirmed by ROSAT 
HRI monitoring obser vations of PSRB08 33-45 yielding a pulsed 
fraction of about 12% iSeward et alj2000h . The soft X-ray (0.1-2.4 
keV) pulse shape turned out to be composed of one broad and two 
narrower pulses. 

In the meantime Chandra and XMM-Newton observations 
considerably improved the overall soft X-ray picture of the Vela 
pulsar and its PWN. The subarcsecond angular resolution of Chan¬ 
dra made it possible to separate the pulsar from its bright PWN, 
and to determine the genuine underlying pulsa r spectrum. The fea- 
tureless phase-averaged pulsar spectmm (see [Pavlov et al.[[2001aL 
using a combination of HRC-S/LETG and ACIS-S/HETG in CC- 
mode data) could best be described in terms of a two component 
model composed of a Hydrogen NS atmosphere model (to describe 
the soft part) and a power-law model with best fit photon-index of 
— 1.5 ± 0.3, nicely connecting to both the optical part and high- 
energy part of the Vela pulsar spectrum. 

Monitoring observations of the Vela PWN with the ACIS-S 
show a highly dynamical system with a complex variability and 
morphological changes of structures (e. g. inner/outer and counter 
jets; inner/outer ar cs; shell) in the PWN 1 Pavlov et alj2001b[ . [2003l ; 
fourant et aljTOljh . 

In the timing domain Chandra observations with the HRC-I 
instrument showed the tripled peaked (energy integrated) soft X- 
ray profile in d etail and aligned acc ording to a contemporaneous ra¬ 
dio ephemeris ([Helfand et al.l200l]) . These authors derived a pulsed 
fraction of 7.1 ± 1.1 across the 0.1-10 keV band (gravity point near 
1 keV) not diluted by PWN emission. 

The Vela pulsar/PWN region was also observed (early in 
the mission on Dec. 1-2, 200(0) by the (large collecting area) 
XMM-Newton satellite with the EPIC pn operating in Small Win¬ 
dow mode at a temporal resolution of 5.7 ms, amply sufficient to 
study the energy d ependency of the pulse morphology (see e.g. 
[Manzali et al.[[200^ . for phase-resolved spectroscopy; however, the 
moderate angular resolution of about 676 FWHM and 10" extrac¬ 
tion radius used, imply a severe (not accounted for) PWN nebula 
contribution at every phase slice making their published (phase- 
resolved) spectral results unreliable). 

We also analysed the XMM-Newton Vela observation of Dec. 
2, 2000 (obs. id. 0111080201; 58.6 ks pn time) to study the energy 
dependency of the pulse profile. The radio-aligned profiles for the 
0.2-0.8 and 0.8-2 keV ranges are shown in Fig.|4]panels e^(the 2- 
8 keV band pn profile is statistically equivalent to the RXTE PGA 


Since April 27, 2006 there is a yearly monitoring of the Vela pulsar region 
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Pulse Phase 

Figure 4. PSRB0833-45; multi-wavelength radio-phase aligned pulse-profile col lage from radio-freq uencies up to high-energy 7 -rays. The left 4 panels (a-d) 
show th e Tow’ energy profi les: a) radio profile at 13.8 GHz, b) optical profile iGouiffes, C. Ill998li and c—d) NUV (1400-3270 A) & FUV (1140-1730 A) 
profiles jRomani et al.ll200^ . Panels e—i show the soft-hard X-ray pulse-profiles as obtained from XMM EPIC pn (panels e—f) and RXTE P CA (panels g-/) 
data. In panels j-l the soft 7 -ray profiles are displayed obtained from multi-year INTEGRAL ISGRI (panels j-k) and CGRO OSSE (panel /. IStrickman et al J 
1 1996ll data. Finally, the right hand panels m-p show the medium- and high-energy 7 -ray pulse profiles as obtained by analyzing multi-year CGRO COMPTEL 
(panels m-n, this work), CGRO EGRET (panel o) and Eermi EAT (panel p) data. Notice the dramatic morphology changes as function of energy from single 
peaked (radio) to double peaked with bridge emission at high-energy 7 -rays, with at intermediate energies very complex shapes with a plethora of pulses at 
optical. X-rays and soft/medium energy 7 -rays. 


2-8 profile and therefore is not shown here). Drastic changes can 
be discerned between the two - soft and intermediate (transition to 
the non-thermal regime) - bands. 

At medium energy X-rays (RXTE PCA; ~ 2 — 30 keV) the 
morphology of the lightcurve drastically changed with two (nar- 
row) peaks roughly alig ned in phase with the high-energy 7 -pulses 
JStrickman et alj[l999t 93 ks obs. performed during Jan 12-23, 
1997). Much deeper RXTE PCA observations (274 ks exposure col¬ 
lected during April-May and July-August 1998), yielding much- 
improved statistic s, showed the presen ce of several (narrow) emis¬ 
sion components jHarding et alj20()3) . In this work we reanalyzed 
all available RXTE PCA data on PSRB0833-45, totaling 360 ks 
screened exposure time (for PCA unit-2). The resulting lightcurves 
for the 2-4, 4-8 and 8-24 keV bands are shown in Fig.|4]panels g-i. 
Notice the presence of at least 5 (narrow) pulse components, one 
coincident with the radio and narrow optical/NUV/FUV pulse at 
phase zero, rendering the Vela pulsar as the most complex pulsar at 
X-ray energies. 

Preceded by non-detections jKnight et alJI 19821 HEAD 7; 15 
keV - 11 M eV) and controversial (short duration) balloon flight 
results (e.g. llTarnden et al.lll972L iTiimer et^ 1 1984 ISacco et S] 


Il990h . the first unambiguous detection (4.6cr in the 70-600 keV 
band) of the Vela pulsar at soft 7 -rays was mad e by the OSSE in¬ 
strument aboard CGRO dStrickman et al.|[l996h . The pulse profile 
(see Fig.|4]panel 1; 70-600 keV) showed two peaks - a narrow and 
a more structured broad one - coincident with the pulses seen at 
high-energy 7 -rays. 

The ISGRI instrument (~ 20-300 keV) aboard the INTEGRAL 
satellite offered a new look at the Vela pulsar at soft 7 -rays. In this 
work a timing analysis of 2244 INTEGRAL ISGRI observations 
of the Vela region, spread over the period June 12, 2003 - Oct. 
13, 2010 (INTEGRAL Revolution period 81-976), has been per¬ 
formed to shed light on the pulse morphology in the so far poorly 
explored soft 7 -ray band. The combined ISGRI dataset represents 
an effective on-axis exposure time of 6.12 Ms. The pulse profiles 
of PSRB0833-45, generated through folding upon proper radio- or 
Fermi LAT based timing models, are shown in Fig. |4] panels j-k 
for the 20-50 and 50-120 keV bands, yielding non-uniformity sig¬ 
nificances of 5a and 4.9(j, respectively, adopting a Zn test. The 
20-300 keV pulse profile represents even a 7.9a significance, with 
one prominent narrow first pulse and a structured second broader 
pulse. 



























































































The soft 'j-ray pulsar population 11 


At medium 7 -ray energies 1-30 MeV), using data from 
the COMPTEL telescope aboard CGRO dSchbnfelder et al.lll993h . 
secure detections of the Vela pulsar have been reported by 
iHermsen et alj dl993L[l9^ and lKuiner et alj dl998h . In the latter 
work phase-resolved spectra were shown combining COMPTEL 
MeV data collected from the first four observation cycli of CGRO. 
In this work we present some new, unpublished so far, COMPTEL 
Vela pulsar results based on observations performed over the full 
CGRO mission lifetime, spanning nine observation cycli and cover¬ 
ing the period May 10, 1991-March 21, 2000 (23 viewing periods, 
typically lasting 2 weeks, with PSRB0833-45 within 30° from the 
pointing axis). 

Applying event selections using optimized parameter win¬ 
dows on Time of Elight (TOE), Puls e Shape Discrimina tion (PSD) 
and (jiarm (see e.g. sections 2 & 3 of iKuiper et alJl 200 lL and refer¬ 
ences therein) and sorting the barycentered timetags on pulse phase 
and energy resulted in pulse-phase distributions as shown in Fig.|4] 
panels m-n, for the 0.75-10 and 10-30 MeV bands, respectively. 
The Zg-test significances represent a 6.2 (t and 12.7a deviation 
from uniformity for both ranges, respectively. Clearly visible are 
emission peaks coincident with the two high-energy 7 -ray pulses, 
PI and P2, with in between bridge emission (most notably in the 
0.75-10 MeV band). Indications for emission structures in the off- 
pulse region exist for energies above 3 MeV. 

PSRB0833-45 was the first pul sar recognized in high - energ y 


7 -ray data from the SAS-2 satellite l lThompson et al 197^ 1977), 


followed by more detailed studies with COS- B 1 Bennett et al.l 
I 1977 I : iKanbach et alj 119801: ICrenier et al.l Il988h and later with 
CGRO EGRET jKanbach et alJll994 Ipierro et alJll^ . The 7 - 


ray pulse profile shows two pulses with clear bridge emission with 
the main pulse trailing the radio-pulse ~ 0.12 in phase (see e.g. 
Fig. |4]panel p for a > 100 MeV Fermi LAT profile). Fermi LAT 
made high-statistics (phase-resolved) studies at high-energy 7 -rays 
possible, revealing a third component in the bridge region, which 
moves towards higher phase with increasing en ergy. The two 7 -ray 
pulses become narrower the higher the energy l lAbdo et ^l2009a . 
l2010dh with an even vanishing first (main) pulse component. This 
has been confirmed now by the upgraded TeV telescope systeem 
H.E.S.S. (phase 11, including a fifth 28m telescope at the center lo¬ 
cation; energy threshold about 30 GeV), yie lding a ~ 11.Oa pulse d 
signal (H-test) at a mean energy of 40 GeV dStegmann et al.l2014h . 

Considering the full pulse-profile collage shown in Fig.|4]it is 
clear that PSRB0833-45 shows extremely complex emission pat¬ 
terns, which change drastically as a function of energy. Therefore, it 
will be a major theoretical challenge to understand the multi-beam 
emission structures of PSRB0833-45 and so to unveil the topology 
and plasma filling of its magnetosphere. 

In Fig. the total pulsed flux of PSRB0833-45 across the 
~ 2 keV - 10 GeV band (thermal part excluded) is shown in blue 
color as solid curves and data s ymbols. Published data and spec¬ 
tral fit results from CGRO OSSE dStrickman et al.ll996h and Fermi 
LAT dAbdo et alj2009cl) are shown as well as newly derived pulsed 
fluxes for the ~ 2- 30 keV (RXTEPCA), 20-175 keV {INTEGRAL 
ISGRI) and 0.75-30 MeV {CGRO COMPTEL) bands. 

Given the complexity of the Vela lightcurve at X-rays (see e.g. 
Fig. |4] panel h) it is rather difficult to define the unpulsed or DC 
level. For the RXTE data we have used the phase stretch 0.68-0.78 
to determine this level. A spectral analysis adopting a power-law 
model with Nh fixed to 3.3 x 10 ^° cm“^ (see e.g. IPavlov et'H] 
l2001al . Table 1) yielded a photon index of —1.06 ± 0.05 across 
the 2.2-28.1 keV band. Pulsed fluxes (unabsorbed) are: 2-10 keV 
(8.0 ± 1.0) • 10”^® erg/cm^s {RXTE PCA); 20-100 keV (1.1 ± 


0.4) • 10"“ erg/cm^s {INTEGRAL ISGR0); 1-10 MeV (6.2 ± 
0.8) ■ 10“^'^ erg/cm^s {CGRO COMPTEL; systematic errors are of 
the order of 30% because of uncertainties in the unpulsed level); 
0.1-1 GeV (3.92 ± 0.49) • 10“® erg/cm^s {Eermi LAT). 

The Vela pulsar is often considered as the canonical high- 
energy 7 -ray pulsar reaching maximum luminosity at GeV en¬ 
ergies. Note, however, that the detection at hard X-rays/soft 7 - 
rays is just possible because the pulsar is relatively nearby. Would 
PSRB0833-45 be placed at 2 kpc, as the Crab pulsar, then the 
pulsed hard X-ray/soft 7 -ray emission would reach flux levels be¬ 
low the detection limits of current instrumentation. 


5.6 PSR J1400-6325 / IGR J14003-6326 

IGR J14003-6326 was discovered bv iKeek et alj ( l2Q06h in a deep 
1.9 Ms INTEGRAL ISGRI mosaic image of the Circinus region 
for the 20-60 keV band using INTEGRAL observations performed 
between Febr. 28, 2003 and Mar. 3, 2005. A 20-60 keV flux of 
0.97 ± 0.15 mCrab and a positional error of 3(6 at 90% confidence 
were derived. 

A ~ 5 ks follow up observation with Chandra on June 
29, 2008 revealed the presence of an extended X-ray source, 
nearly circular with a 1(5 radius, within the ISGRI error region 
jTomsick et alj[200^ . Its morphology with at the center region a 
PWN surroundin g a so-far unde t ected pulsar made a SNR identi¬ 
fication plausible. iTomsick et al.l d2009l) derived spectral character¬ 
istics for the composite emission region within 1(5 from the cen¬ 
tral source by fitting an absorbed power-law model over the 0.3- 
10 keV range, and obtained an absorbing Hydrogen column Nh of 
(3.1 ± 0.3) X 10^^ cm“^, a photon index of —1.83 ± 0.13 and an 
unabsorbed 2-10 keV flux of 1.9 x 10“^^ erg/cm^s (~ 1 mCrab). 

iRenaud et alj feOlOh reanalyzed the Chandra data and decom¬ 
posed the emission into 3 components: 1) a point-source, extracted 
for a circular region of 1'.'5 centered on source with background 
taken from an annulus with radii l('5 and 2”, 2) a PWN-region 
extracted from an annulus with radii T.'5 and 39” and 3) the SNR- 
region extracted from an annulus with radii 0(98 and 1(27. Assum¬ 
ing absorbed power-law models they found the following photon- 
indices and unabsorbed 2-10 keV fluxes for an absorbing column 
of (2.09 ± 0.12) X 10^^ cm-^: _1.22 ± 0.15, -1.83 ± 0.09 and 
-2.56 ± 0.18, and (1.95 ± 0.5) x 10■^^ (1.51 ± 0.2) x 10"“ 
and (1.1 ± 0.2) x 10“^^ erg/cm^s for the 3 components, respec¬ 
tively. The PWN is thus by far the most dominant source, and about 
8 times brighter than the hard point-source. 

The detection of pulsed X-ray emission of IGR 114003-6326 
was also reported by iRenaud et alj ( l2010h using 2 sets of RXTE 
PCA observations of 44 ks each, taken one year apart on Sept. 29, 
2008 and Sept. 30, 2009. They found an energetic (75 = 5.1 x 10®^ 
erg/s) pulsar spinning with a 31.2 ms pulse period and an esti¬ 
mated age of 12.7 kyr. This period was later confirmed at radio¬ 
frequencies with the Parkes telescope, resulting furthermore in a 
(very) high DM value of 563 ± 4 cm~® pc and flux densities of 250 
and 110 /rJy at 1.3 7 4 and 3.1 GHz, respectively. 

IRenaud et alj (1201 Ol) also derived the pulsed spectrum of 
IGR J14003-6326, through an on-off method, for the 2-10 keV 
band using only the top detection layer of each involved PCU. They 


® Because of the complexity of the lightcurve we fitted the 20-100 keV 
profile in terms of a constant and 7 harmonics to estimate the unpulsed 
level, and so the number of pulsed excess counts. Varying the number of 
harmonics, 5, 7, 9, 11, yielded consistent results. 
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Pulse Phase 

Figure 5. IGR J14003-6326; pulse profiles for RXTE PCA in different en¬ 
ergy bands. The significance in the 12.2-32.2 keV band is S.Tct (panel c). 


found an absorbed 2-10 keV pulsed flux of 3.0 x 10~^® erg/cm^s 
and a (soft) photon index of — 2 . 0 lg’ 3 . 

We reanalyzed the PCA data including now all PCU detector 
layers to be more sensitive to the hard X-ray photons with ener¬ 
gies in excess of 10 keV. Pulse profiles for three PCA bands are 
shown in Fig. |5] Significant pulsed emission at ~ 3.7cr is found 
for the 12.2- 32.2 keV band. Pulsed excess counts have been de¬ 
rived through pulse template fitting using as extractor the 3 har¬ 
monics fit to the lightcurve of the 4-87 PHA band. These excess 
counts have been converted to pulsed fluxes adopting an absorb- 
ing Hydrogen column of 2.09 x 10^^ cm“^ jRenaud et al.ll2O10|) 
and assuming an underlying (absorbed) power-law model for the 
photon fl ux. We obtained a p hoton index of —1.95 ± 0.04, consis¬ 
tent with lRenaud et alJOOlOh . and an unabsorbed 2-10 keV flux of 
(1.3±0.1) X 10~^^ erg/cm^s, about ~ 4 times larger than found by 
[Renaud et al.l ( l2010h . Our pulsed flux repre sents ~ 67% of the total 
2-10 keV flux, while [^naud et af] ||2010|) found a pulsed fraction 
of only 18%. The pulsed flux measurements of 1GRJ14003- 
6326 for 2 — 32 keV band are shown in Fig. (left panel) as 
yellow data points. 

We also analyzed HEXTE data. As expected, we could not 
find evidence for pulsed emission for energies in excess of ~ 15 
keV. 

In the TeV domain no point-like source has been detected 
abov e 5a in the latest map o f the H.E.S.S. Galactic Plane Sur¬ 
vey dChaves et alJl2008L l2009h . Also at the Ge V energies neither 
a point-like source is found in Fermi EAT data jNolan et alj|2012ti 
at the loc ation of IGR J140 03-6326 nor pulsed emission has been 
detected dAbdo et^l2013h . 


5.7 PSR B1509-58 / PSR J1513-5908 


PSR B1509-58 was discovered at X-ra ys as a 150 ms pulsar in Ein - 
stein HRl and IPC (0.2-4 keV) data bv ISeward & HarndenI (Il982h 
and la ter confirmed at radio-wavelengths by Manchester et al.l 
( Il982h . Its characteristic age of only 1570 years and spin-down lu¬ 
minosity of 1.8 X 10®^ erg/s, as inferred from its timing parameters, 
qualifies PSR B1509-58 as a very young and energetic spin-down 
powered pulsar. Its derived surface polar magnetic field strength 
of 3.1 X 10^® G is close to the quantum critical field strength 
of 4.413 X 10^® G, above which exotic quantum electrodynamic 
effects play an important role in the high-energy pulsar emission 
mechanisms. 

PSR B1509-58 is associated with supernova remnant MSH 
15-52 (G320.4-1.2), and various detailed studies using X-ray data 
have been performed since the early eighties to disentangle the 
complex morphology of MSH 15-52 with at the north west¬ 
ern rim an excess near Ha nebula RCW 89 and near the cen¬ 
tre a clump containing a diffuse synchrotron nebula surrounding 
the pulsar ( its PWN). Using the superb X-ray imaging quality 
of Chandra iGaensler et al.l ( l2002li showed that, while the overall 
PWN shows a clear symmetry axis, the PWN also reveals sev¬ 
eral new components. Extended TeV emission (> 280 GeV) from 
MSH 15-52 was detected at a 25a confidence level by H.E.S.S. 
jAharonian et alj|2005bll . and later at hi gher TeV 7 -ray s (> 810 
GeV) by CANGAROO-111 at a 7a level iNakamori et'^l2008h . In 
the GeV ba nd extended 7 -rav emission has been detected by the 
Fermi EAT dAbdo et^bOlOa) . 

Pulsed emi s sion i n the medium/hard X-ray band was detected 
by iKawai et alj dl99lh using Ginga (2-60 keV) data. These au¬ 
thors also found that, using contemporaneous X-ray (Ginga) and 
radio (Parkes) data, the radio pulse leads the asymmetric X-ray 
pulse by 0.25 ± 0.02 in phase. Soon after the launch of the Comp¬ 
ton Gamma-Ray Observatory (CGRO) on April 16, 1991, pulsed 
hard X-ray/soft 7 -ray emission above ~ 60 keV was detected 
by BATSE using dat a collected in the fol ded-on-board data m ode 
I Wilson et^ll993al f9). confirmed later bv lUlmer et al.l (Il993h and 
iMatz et alJ ( Il994h using CGRO OS SE (0.05-10 MeV) data. The 
radio- soft 7 -ray pulse lag found by lUlmer et alJ ( Il993h combin¬ 
ing BATSE and OSSE data was 0.32 ± 0.02, somewhat larger than 
the lag found at lower energies. In the spectral domain iMatz et al.l 
(Il994 ) derived for the pulsed spectrum from 50 keV to 5 MeV, 
fitting a power-law model, a photon index of —1.68 ± 0.09, signif¬ 
icantly softer than the valu e of — 1.30 ± 0.05 obtained for the 2-60 
keV band by iKawai et al.l ( Il99 ih analyzing Ginga data, strongly 
indicating a softening towards higher energies. 

Pulsed emission was also detected in the medium energy 7 - 
ray band bv lKuiper et Zl ( 1 19991) using CGRO COMPTEL (0.75-30 
MeV) data. They found a 5.4 ct pulsed signal, a broad asymmetric 
pulse, for the 0.75-30 MeV band analyzing COMPTEL data col¬ 
lected over 6 CGRO observation cycli. The pulse reaches its max¬ 
imum at phase 0.38 ± 0.03 with respect to the radio pulse. The 
apparent shift of the pulse maxima from soft X-rays to medium 
energy 7 -rays can be explained by different spectral behaviours 
of two Gaussian shaped pulses comprising the broad asymmetric 
pulse with one pulse peaking at 0.250 ± 0.008 and the other at 
0.386 ± 0.012 as derived from high-statistics RXTE PCA data. 
The first (narrower) pulse seems to vanish towards higher energies 
dKuiper et al.lll999l) . The spectrum of the pulsed emission clearly 
breaks above ~ 10 MeV, though there were weak indications for 
pulsed emission beyond 10 MeV in both COMPTEL and CGRO 
EGRET data (also in the spatial domain). 
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Figure 6 . PSRB1509-58 (the prototype soft 7 -pulsar); pulse profiles 
from 2 keV up to 1 GeV using multi-year RXTE PCA & HEXTE data (pan¬ 
els a-c), CGRO BATSE (Cycli I-IX; panel d) & COMPTEL (Cycli I-VI; 
panel e) and Fermi EAT (covering 3.4 years; panel/). 


At high-energy 7 -rays (> 100 MeV ) significant progres s has 
been made a nalyzing Fermi LAT d ata dAbdo e t al.| 1201 Obi) and 
AGILE data jPellizzoni et alj l2009l ; IPilia et alT 201(1) by detect- 
in g pulsed emission from PSRB1509-58 well above 100 MeV. 
In Iden Hartog et all ll2014l) the Fermi LAT results are presented 
analyzing ~ 5 years of data detecting pulsed emission from 
PSRB1509-58 at a 18.2cr level across the ~ 30 - 1000 MeV band, 
enabling detailed studies of the morphology of the lightcurve and 
the breaking part of the spectrum at these energies. 

To achieve pulse profile morphology comparisons across the 
~ 2 keV - 1 GeV range at a high statistics level we analyzed and 
processed multi-year data from RXTE PCA (2-30 keV; April 25, 
2002 - Oct. 9, 2003), RXTE HEXTE (15-250 keV; Eeb. 10, 1999 - 
Jan. 1, 2012), CGRO BATSE (22 - 4000 keV; Oct. 31, 1991 - May 
9, 2000; MID 48560-51673; CGRO Cycli I-IX), CGRO COMP¬ 
TEL (0.75-30 MeV; Oct. 17, 1991 - Nov. 3, 1997; MJD 48546.6- 
50755.6; CGRO Cycli I-VI) and Fermi LAT (30-1000 MeV; Aug. 
4 18:00:00 - Jan. 1, 2012; MJD 54682.75-55927). 

RXTE PCA data (in goodXenon mode) were collected dur¬ 
ing observation run 70701 and part of 80803. The time intervals 
are fully covered by the validity interval of an ephemeris from 
the ATNE archive ensuring proper timing information. Subsequent 
screening adopting default criteria rendered for PCU-2 an effective 
exposure of 74.53 ks. This yielded ample statistics to study its tim¬ 
ing characteristics in great detail, because PSRB1509-58 is a rela¬ 


tively strong pulsar at medium/hard X-ray energies. To determine 
the pulsed spectrum across the 15-250 keV band we used RXTE 
HEXTE data collected over a much longer time period than used for 
the PCA. We combined all available HEXTE data on PSRB1509- 
58 taken during observation cycli P40704 - 96803 (Eeb. 10, 1999 
- Jan. 1, 2012) to obtain the highest possible statistics. Eor more 
information on these RXTE PCA and HEXTE data we refer to the 
section on RXTE in Iden Hartog et al.l ( |2014|) . The resulting pulse 
profiles for the PCA (2-20 keV) and HEXTE (20-100 and 100-250 
keV) are shown in panels a-c of Eig.|^ The Z|-test significance of 
the 100-250 keV profile is about IT.lcr. 

To study the pulse profile of PSRB 150 9-58 at soff/medium 
energy 7 -rays we analyzed BATSE (see e.g. lEishman et al.|[l9^ 
and references therein) LAD folded-on-board data types collected 
during the full duration of the CGRO mission yielding a net 
screened exposure of about 10.6 Ms (integrated sum for all 8 
LADs) for the time period Oct. 31, 1991 - M ay 9, 2000. The pro- 
cessing of these data is analoguous to that in lKuiper et al.l ( 1200 ll) 
for the Crab pulsar (see Section 5.3 of that paper). The resulting 
317-1102 keV BATSE profile of PSRB1509-58 is shown in panel 
d of Eig.|^ The deviation from a statistically uniform distribution is 
about 10.4cr adopting a Z|-test. Panel e o f Pig.[ 6 ]shows the C OMP- 
TEL 0.75-30 MeV profile as obtained bv lKuipCT et al.l (Il999l) . 

Einally, panel f of Eig.| 6 ]shows the Fermi LAT 30-1000 MeV 
lightcurve using data from a 5° aperture around PSR B1509-58 col¬ 
lected during Aug. 4, 2008 18:00:00 and Jan. 1, 2012 (about 3.4 
years) obtained after folding the SSB corrected arrival times with 
RXTE PCA based timing models. The Z^-test significance of the 
30-1000 MeV profile is about 10.2cr, and a single Gaussian can de¬ 
scribe the measured distribution accurately. Pulse maximum occurs 
at phase 0.342±0.007 and the width (Gaussian a) is 0.064±0.007. 
This indicates that the first pulse indeed vanishes, while the second 
pulse becomes somewhat narrower, when the energy increases. 

Equiped with the multi-year RXTE PCA/HEXTE data cover¬ 
ing at high statistics the ~ 2 — 250 keV band, the m ulti-year Fermi 
LAT (30-1000 MeV) data dden Hartog et al.ll2014h . s upplemented 
by th e (archival) COMPTEL (0.75-30 MeV) data dKuiper et al.l 
I 1999 I) and INTEGRAL ISGRI (20-175 keV; 550 ks of dedicated 
PSRB 1509-58 observation dmmg INTEGRAL AOl', this work), we 
attempted a broad-band spectral fit of the pulsed emission across 
the ~ 2 keV - 1 GeV range adopting a photon flux model (a kind 
of power-law wifh a (modified) exponential cutoff) with the fol¬ 
lowing E.y dependency: 

F^ = k- {E-,/Eof ■ eM-{E-,/E,f) 

Interstellar absorption effects, only applicable to RXTE PCA data, 
are modeled out a dopting an absorbing Hydrogen column Nh of 
9.5 X 10^^ cm“^ dGaensler et al.ll200a) . Eor our spectral dataset 
the normalization energy Eq, which minimizes the correlation be¬ 
tween the 4 fit parameters, was 0.100518 MeV, while the best fit 
ixl, 59_4 = 1.35 i.e. acceptable) yielded the following fit parame¬ 
ters, k = (1.574 ± 0.012) • 10"^ ph/cm^s MeV; T = -1.233 ± 
0.005; Ec = 0.078 ± 0.003 MeV and /3 = 0.286 ± 0.005. The 
broad-band ~ 2 keV - 1 GeV pulsed spectrum of PSRB1509-58 
is shown in Eig.|7] along with the best fit model. Maximum lumi¬ 
nosity per energy decade is reached at ~ 2.5 MeV and therefore 
PSRB1509-58 can be considered as the “canonical” soft 7 -ray pul¬ 
sar. 

Pulsed fluxes (unabsorbed), obtained from the best fit model, 
in the 2-10 keV, 20-100 keV, 1-10 MeV and 0.1-1 GeV are, 
(2.52 ± 0.05) ■ 10”“, (9.59 ± 0.13) • 10"“, (4.4 ± 0.3) • 10"^° 
and (1.2 ± 0.4) • 10“^^ erg/cm^s, respectively. 
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Figure 7. PSRB1509-58; unabsorbed high-energy spectrum of the pulsed 
emission across the 2 keV - 1 GeV range, combining RXTE PCA, 
HEXTE, INTEGRAL ISGRI, CGRO COMPTEL and Fermi EAT measure¬ 
ments (see text). Maximum luminosity is reached at about 2.5 MeV making 
PSRB1509-58 an archetypal soft 7 -ray pulsar. 


5.8 PSRJ1617-5055 

iTorii et all il998l) discovered an energetic 69 ms pulsar, 
PSRJ1617-5055, - 7' from SN-remnant RCW 103, analysing X- 
ray data obtained with the ASCA GIS instrument. It is noteworthy 
that lAoki et ^ lll992h already detected pulsed emission from the 
RCW 103 region with the Ginga large area counters in data taken 
in March 1989. 

Weak 0.5 mj y at 1.4 GHz) radio pulsations were detected 
bv lKasni et al.l l ll998h using the Parkes radio telescope. Combining 
more pulse period measurements from X-ray observations a char¬ 
acteristic age of about 8 ky r could be derive d as well as strong in¬ 
dications for a giant glitch iTorii et alJl 2000 ll that occured between 
August 1993 and September 1997. 

RXTE had the sky region containing PSR11617-5055 several 
times in its field of view: on Jan. 2, 1998 and during March 5-6, 
1999 in observation run 30210 with prime target IE 161348-5055 
(the point source in the middle of RCW 103), on Jan. 22, 2001 dur¬ 
ing run 50428 with prime target RCW 103, and finally in a dedi¬ 
cated observation run, 80090, performed during September 20-21, 
2003 for about 80 ks. We combined RXTE PCA (total screened 
PCU-2 exposure 146.53 ks) and HEXTE (dead-time and off-axis 
corrected cluster-0/1 exposures of 47.7 and 49.3 ks, respectively) 
data from these 3 observation runs and significantly detected pulsed 
X-ray emission up to 64 keV. Pulse profiles in differential energy 
bands, all single peaked with a sharper rise than fall, are shown in 
Fig.H] The HEXTE profile for energies in excess of 35.3 keV (Eig. 
m has a Zi significance of 3.8cr. 

We estimated the pulsed excess counts through pulse-profile 
model fitting with a truncated Fourier series (fundamental plus 2 
harmonics) and converted these to photon flux values, adopting 
an absorbed power-la w model with a fixed colum n density Nh 
of 3.2 X 10^^ cm“^ teecker & Aschenbachll20()2h . This yielded 
for the pulsed emission in the 2.5-30 keV PCA band a hard pho¬ 
ton index of —1.42 ± 0.02 and a 2-10 keV unabsorbed flux of 
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Figure 8. PSR J1617-5055; pulse profiles in various energy bands for 
RXTE PCA (panels a-d) and RXTE HEXTE (panels e-f). Pulsed emission 
has been detected up to ~ 64 keV. 


(3.30 ± 0.16) X 10“^^ erg/cm^s. We used all PCA detection lay¬ 
ers to be more sensitive to hard X-ray photons. The reconstructed 
RXTE PCA pulsed flux values are shown in Fig. as aqua colored 
points along with the best fit. HEXTE pulsed flux measurements 
are also shown (blue squares). 

Unfortunately, neither X-ray monitoring- nor regular radio ob¬ 
servations (too weak, requiring long exposure times) have been per¬ 
formed on this pulsar, that would have enabled us to construct phase 
coherent timing models over long time stretches. Therefore, we 
could not do a timing analysis for the soft 7 - band using t h e low- 
countrate INTEGRAL ISGRI data. However, lUandi et ^ ( l2007h 
demonstrated in a spatial analysis in ISGRI maps the presence of 
a hard X-ray source coincident with PSR 11617-5055 with a sig¬ 
nificance in the 18-60 keV band of ~ 7.4 (t, combining data col¬ 
lected since the beginning of the mission (November 2002) up to 
April 2006. Our own spectral analysis of ISGRI data (Revs. 46-411; 
March 2, 2003 - Feb. 24, 2006) yielded for the total emission in the 
20-300 keV band a photon index of —1.93 ± 0.28 and a 20-100 
keV flux of (1.25 ± 0.18) x 1 0~^^ erg/cm^s, consistent with the 
values obtained bv iLandi et alj ( l2007h . The ISGRI total flux mea¬ 
surements in the 20-300 keV band are indicated in Fig.[9]as purple 
datapoints along with the best power-law fit. 

The total emission spectrum of PSR 1 1 617-5055 below lOkeV 
was a lre ady derived by Tpriiet^ (l200(]h: [Becker & Aschenbachl 
( |2002|) : iLandi et^ l l2007h : Kargaltsev et al. ( l2009h using X- 
ray instruments on different spacecrafts - ASCA, XMM-Newton, 
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to determine the local background, we derived a genuine pulsed 
fraction (0.2-10 keV) of 78±2%± 10% taking into account the 10- 
15 % contribution of the inner PWN to the total emission from the 
8 ” circular (source) extraction region, which could not be resolved 
in this observation. The first error is related to the uncertainty in the 
inner PWN contribution and the second one to the uncertainty in the 
number of pulsed coun ts, derived through the method outlined by 
ISwanepoel et al.l l ll996ll . 

Also, above 20 keV there is very little room for pulsar DC 
and/or PWN emission. At GeV energies no pulsed high-energy 
7 -rav emission h as been detected so far for PSR J1617-5055 
dAbdo et al.ll2013h . 

At TeV energies (> 200 GeV) a bright source, HESS J1616- 
508, has been detected in the neighbourhood of PSR J1617-5055 
jAharonian et alj [20061) . However, due to (severe) misalignment 
a convincing identification at lower energies is lacking (see e.g. 
iMatsumoto et alJl2007l : lKargaltsev et alj2009ll . 

5.9 PSR J1640-4631 in SNR G338.3-0.0/HESS J1640-465 


E [MeV] 

Figure 9. PSRJ1617-S055; high-energy (0.8-300 keV) total and pulsed 
spectra as derive d from measurem e nts by Chandra ACIS (total pulsar 0.8- 
8 keV; solid red; iKargaltsev et all d2009r) i. XMM-Newton MOS (total, in¬ 
cluding compact PWN, 0.8-10 keV; solid orange; [Becker & Aschenbachl 
i2002t) i. joint fit XMM-Newton/B eppoSAX/lSGRl (to tal, including compact 
PWN, 0.8-10 keV; solid brown; iLandi et alj d2007h i. RXTE PCA (pulsed 
2.5-30 keV; aqua; this work), RXTE HEXTE (pulsed 15-150 keV; blue; this 
work) and INTEGRAL ISGRI (total 20-300 keV; purple; this work). 


BeppoSAX and Chan dra. The most recent determination by 
IKargaltsev et alJ ll2009l) . using Chandra ACIS-I data, resolves the 
emission around PSRJ1617-5055 in a point-like component, asso¬ 
ciated with the pulsar, an inner compact pulsar wind nebula compo¬ 
nent of ~ 1" size within (T.'S and 1'.'5 from the pulsar and an outer 
nebula extending up to l' from the pulsar. Fitting a power-law 
to the spectral data extra cted from a circular ape rture of 0”5 cen¬ 
tered on PSR J1617-5055. [Kargaltsev et al.U20(i9l) derived a power- 
law index of —1.14 ± 0.06 and a 0.5-8 keV (unabsorbed) flux of 
(3.6 ± 0.1) X 10“^^ erg/cm^s for the total (=pulsed-l-pulsar DC) 
emission of the pulsar. This converts to a 2-10 keV (unabsorbed) 
flux of 3.58 X 10“^^ erg/cm^s. The flux of the inner compact 
PWN is about 10-15 % of the total pulsar flux, and its combination 
is consistent with the total emission spectra deri ved using the X-ray 
instruments with much less spatial resolution; iLandi et ^ l l2007r) 
find a 2-10 keV flux of 4.2 x 10“^'^ erg/cm^s using BeppoSAX 
MEGS, XMM-Newton MOS 1+2 and INTEGRAL ISGRI in com¬ 
bination, wh ile the 2-10 keV flux, estimat ed from the best fit val¬ 
ues given bv iBecker & Aschenbachl | |2002|) , amounts 4.0 x 10~^^ 
erg/cm ^ s. The best fit mo d els derived by [Becker & Aschenbachl 
( l 2002 h ; ILandi et al.l j2007h ; IKargaltsev et alj (l2009l) are shown as 
solid orange, brown and red, respectively, lines in Fig.j^ 

Comparing the pulsed (RXTE PCA; this work) and total 2-10 
keV flux values (Chandra', MOS 1 - 1 - 2 ; BeppoSAX) yields a (high) 
pulsed fraction in the range 79-92 %. This value is consistent with 
the value, derived in this work, using data from merely Chandra 
HRC-S in timing mode. During a 78.2 ks exposure of IE 161348- 
5055 on July 2, 2007 (obs.id. 7619) with the HRC-S PSRJ1617- 


5055 was observed 7f3 off-axis. Using a 8 ” extraction radius, be¬ 
cause of the strongly degraded PSF, and using an annulus centered 
on the pulsar with inner and outer radii of 12 ” and 16”, respectively, 


HESS J1640-465 was discovered during the Gala ctic plane survey 
with H.E.S. S. performed during May-July 2004 jAharonian et al.l 
1 20051 l2006h . The source is marginally extended at TeV energies 
and its location is consistent with the S'diameter broken-shell 
SNR G338.3-0.0, which lies near the boundary of a bright H ii 
region. At the center of SNR G338.3-0.0 an X-ray source, AX 
J 16 40.7-4632, was detec ted during the ASCA Galactic plane sur¬ 
vey jSugizaki et al.ll 200 TI) . 

HESS 11640-465 was observed by XMM-Newton for 21.8 ks 
on August 20, 2005 wit h the EPIC-pn and MOS camera’s operat¬ 
ing in full-frame mode jpunk et alj|2007bl) . This observation was 
strongly affected by soft proton flares reducing the effective ex¬ 
posure to only 7.3 ks. Three sources were detected in this obser¬ 
vation of which XMM J164045.4-463131 is coincident with AX 
J 1640.7-4632, and is extended in nature with a compact core and 
a faint tail, resembling morphologically a typical PWN. A spectral 
analysis of XMM J164045.4-463131, fitting an absorbed power- 
law model, yielded a rather strong absorbing Hydrogen column Nh 
of (A.l^Qg) X cm“'^ and a photon index of —1.74 ±0.12. No 
shell-like X-ray emission was apparent in the XMM-observation. 

The extended nature of XMM J164045.4-463131 was con¬ 
firmed by lUemiere et alj j2009l) . who analysed a Chandra ACIS 
observation of HESS J1640-465 taken in May 2007 with an effec¬ 
tive exposure time of 26.4 ks. Employing the sub-ar csecond scale 
spatial resolution of Chandra lUemiere et al.l j2009l) performed a 
spatially resolved spectral analysis of the near field of XMM 
J164045.4-463131 and found that the putative pulsar spectrum is 
very hard with photon index —1.1 ± 0.4 heavily absorbed through 
a Hydrogen column Nh of 1.4 x 10^® cm“^. The unabsorbed 2-10 
keV flux of the putative pulsar was 1.5 x 10“^® erg/cm'^s. The com¬ 
pact PWN has a softer photon index of —2.5 ± 0.3 and unabsorbed 
2-10 keV flux of 4.2 x 10“^® erg/cm'^s, while the extended PWN 
is somewhat steeper with index —2.7 ± 0.5 and has an unabsorbed 
2-10 keV flux of 4.6 x 10“^'^ erg/cm'^s. 

Observations with NuSTAR (3-79 keV) of HESS 11640-465 
on June 22, 2013 and Sept. 29, 2013 for 48.6 and 89.9 ks, respec¬ 
tively, r evealed eventually th e pulsed nature of XMM J164045.4- 
463131 jCotthelf et alj2()14l) . An energetic, E ~ 4.4 x 10^® erg/s, 
206 ms pulsar was found with a characteristic age of 3.35 kyr. The 
pulsar shows a single pulse - relatively s harp compared to a sinu- 
soid - profile in the 3-25 keV band (see lGotthelf et al.]|2014L and 
also Fig. Impanel q), and thus can be considered as a soft 7 -ray 
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Figure 10. PSR J1811-1925; pulse profiles in various energy bands for RXTE PCA (panels a-c), RXTE HEXTE (panels d-f) and INTEGRAL ISGRI (panels 
g-i). Pulsed emission has been detected up to ~ 135 keV. 


pulsar. The pulsed fraction can be as high as ~ 82% taking into 
account both the ’normal’ background and PWN contamination. 
The (total) pulsar spectrum combining Chandra and NuSTAR data, 
adopting an absorbed power-law model, has a hard photon index of 
— l-3lo 5 absorbed 2-10 keV flux of (1.8 ± 0.4) x 10“^® 

erg/cm^s, which translates to an unabsorbed 2-10 keV flux of 
(3.63 ± 0.81) X 10“^® erg/cm^s taking into account the very high 
Hydrogen column density Nh of (1.8 ± 0.6) x 10^® cm“^. 

ICotthelf et al] j2014h also searched for pulsed high-energy 
7 -ray emission (> 100 MeV) using Fermi LAT data collected 
over a more than 5 years time period, but found no significant 
pulsed signal. In the sp atial domain, how ever, the Fermi LAT 
data (first reported on by ISlane et al.ll 20 i^ do show two distinct 
sour ces above 100 MeV in the reg ion containing HESS J1640- 
465 dLemoine-Goumard et alj|2014h . In light of the new/updated 
findings of H.E.S.S. TeV measurements on HESS J1640-465 
using all data collected between May 2004 and Sept. 2011 
iIh.E.S.S. Collaboration et aU l2014h and on the recently discov¬ 
ered (hard) nearby HESS J1641-463 (coincident with radio SNR 
G338.5-1-0.1), unnoticed using default detection techniques, but de¬ 
tected using much higher energy thre sholds jAbramowski et all 
l2014h . iLemoine-Goumiard et alJ l2014h revisited the Fermi LAT 
data and reported a power-law spectrum with index —1.99 ± 
0.04 ± 0.07 for HESS J1640-465, which smoothly connects to 
the TeV data points, and a soft power-law spectrum with index 
-2.47 ± 0.05 ± 0.06 for HESS J1641-463, requiring a signifi¬ 


cant hardening towards TeV energies. The new GeV/TeV data sug¬ 
gest that a significant part of the high-energy 7 -ray emission from 
HESS 11640-465 o riginates from proton acceleration i n the super¬ 
nova remnant shell 1 H^RS^. Collaboration et aLll2014ll . 

ICastelletti et"^ ^201117 presented a high-resolution imaging 
study of the G338.3-0.0 supernova remnant at radio frequencies ac¬ 
quired by the Giant Metrewave Radio Telescope GMRT (235, 610 
and 1280 MHz) and ATCA (1290 and 2300 MHz), and searched 
for a radio pulsar at 610 and 1280 MHz in a beam covering 
the XMM 1164045.4-463131 location (coincident with PSR 11640- 
4631). This search yielded no radio pulsar counterparts and resulted 
in flux density upper limits of 2 and 1 mly for 610 and 1280 MHz, 
respectively, at 8 (t confidence level. 


5.10 PSRJ1811-1925 

PSR 11811 -1925 was discovered bv lTorii et alJ (Il997t) during ASCA 
(0.5-10 keV) observations of G11.2-0.3, known to be a composite 
supernova remnant, pro posed to be associa ted with the historical 
supernova of A.D. 386 l lVasisht et alJll996ll . They found a 65 ms 
pulsar exhibiting a hard X-ray spectrum. Subsequent ASCA and 
BeppoSAX observations iTorii et ^Il999h allowed the measure- 
ment of the spin-down and yielded a characteristic age of 2.4 x 10^ 
yr. No pulsed radio emission has so far been detecte d from the sky 
region containing the pulsar iCrawford et al.lll998ll. A high nre- 
cision Chandra X-ray imaging observation ~l lKasDi et al.ll2()0lh re- 
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Figure 11. PSR J1811-1925; high-energy (1-300 keV; unabsorbed) total 
and pulsed spectra as derived from measurements by Chandra ACIS (1-7 
keV; pulsar Pulsed plus DC (=total; orange solid); PWN (orange dotted); 
PWN plus pulsar (orange long dashed)), RXTE PCA (aqua) and HEXTE 
(blue; both pulsed) and INTEGRAL ISGRI (total (brown) / pulsed (pur¬ 
ple)). The ISGRI total flux measurements still contain contributions from 
the PWN. The spectral measurements are generally consistent with reach¬ 
ing a pulsed fraction of 100% near 50 keV. 


vealed at the very center of the SN-remnant the pulsar counterpart 
at R.A. 18'’11"’29^22, Decl. -19°25'27''6 (Epoch J2000), embed¬ 
ded in a diffuse nebula (PWN). This provided strong evidence that 
the system is much younger than the characteristic age suggests. 
Assuming conventional spin-down with a constant magnetic field 
and breaking index n (from the assumed spin-down law: u oc u") 
both ages can only be reconciled if the pulsar was initially bom 
with a spin period of ~ 62 ms unless the breaking index is unusally 
large. 

Subsequent RXTE PCA X-ray monitoring over 992 days 
(March 8 , 2002 - November 24, 2004; MJD 52341 - 53333) al¬ 
lowed us to generate phase coherent timing models (see Table [T) 
and using these to detect th e pulsed signal up to ~ 100 keV in 
RXTE HEXTE data (see also lRoberts etaDl2004 . A second RXTE 
PCA monitoring campaign was initiated on August 25, 2007 and 
ended on October 31, 2011 near the cessation of the RXTE mis¬ 
sion. The data from this latter campaign are not included in this 
study. 

We performed a timing analysis of INTEGRAL ISGRI data 
for the time span covering INTEGRAL revolutions 46-249 (Eebru- 
ary 28, 2003 - October 28, 2004; MJD 52698 - 53306), for which 
we derived accurate ephemerides (see Table [T). We found a 6.5a 
signal for the 15-135 keV energy band showing a morphology con¬ 
sistent with that derived from RXTE PCA and HEXTE data. Pulse 
profiles in differential energy bands for RXTE PCA & HEXTE and 
INTEGRAL ISGRI are shown in Fig.[T0] The significances (adopt¬ 
ing Zj-statistics) for a deviation from a uniform distribution are: 
35.6a, 35.9a and 15.9a for the RXTE PCA ranges 2.8-8.2, 8.2- 
16.5 and 16.5-32.3 keV (panels a-c), respectively; 11.2a, 8 . 6 a and 
4.4a for the RXTE HEXTE ranges 15.6-29.0, 35.2-64.1 and 64.1- 
132.6 keV (panels d-f), respectively, and finally, 3.7a, 3.6a and 


3.4a for the INTEGRAL ISGRI ranges 15.0-30.0, 30.0-60.0 and 
60.0-135.0 keV (panels g-i), respectively. We derived the pulsed 
excess counts in differential energy bands with pulse-profile fit¬ 
ting, which subsequently have been converted to photon fluxes tak¬ 
ing into account the energy responses of the involved instruments. 
These (unabsorbed) pulsed flux measurements are shown in Fig. 
II Has aqua (PCA), blue (HEXTE) and purple (ISGRI) data points, 
along with the best-fit “curved” power-law model (black) resulting 
from a combined fit. Pulsed fluxes derived from this fit for the 1- 
10 and 20-100 keV energy bands are (2.21 ± 0.47) x 10“^^ and 
(1.2 2 ± 0.21) X 1 0~^^ erg/cm^s, respectively. 

[Roberts et al.l ( l2003l) presented detailed results from spatially 
resolved spectral analysis for the point-source (=pulsar) in G11.2- 
0.3, its PWN and other structures in its near environment using 
Chandra ACIS data. The power-law model fits to the spectral data 
of the pulsar (=total emission), PWN and their combination are 
shown in Fig.[TT]for the 1-8 keV band. The listed absorbed (private 
communication M. Roberts) energy flux of the pulsar in the 1-10 
keV band, (2.82 ± 0.12) x 10“'^^ erg/cm^s has been converted 
to its unabsorbed value of (3.48 ± 0.15) x lO”*^^ erg/cm^s using 
the estimated column density Nn of 2.22 x 10^^ cm“^ . Using our 
value for the pulsed flux, the pulsed fraction in the 1-10 keV band 
becomes 0.64 ± 0.14. 

In the imaging domain we detected soft 7 -ray emission up to 
~ 150 keV using archival data {INTEGRAL revolutions 46-495; 
February 28, 2003 - October 5, 2006, 11.35 Ms GTI exposure; 
this work; see also e.g. lSird et al.ll200'7l ; lDean et al.ll2008h . The to¬ 
tal (= pulsed plus unpulsed from PSR J1811-1925 and DC from 
the PWN and SN-remnant) emission spectrum of this source for 
energies above 20 keV can be described by a power-law model 
with a (hard) index of —1.61 ± 0.15 and a 20-100 keV flux of 
(1.54 ± 0.12) X 10“'^^ erg/cm^s (see Fig.[TTJ brown d ata points). 
This i s consistent with the index and flux derived by IPean et Tl 
(l2008h . who used slightly less INTEGRAL exposure on the source. 
The pulsed fraction in the 20-100 keV band is 0.79 ± 0.15, and 
above ~ 50 keV the pulsed PCA/HEXTE/ISGRI spectrum is con¬ 
sistent with the total ISGRI spectrum, thus the pulsed fraction be¬ 
coming consistent with 100% for energies beyond. Therefore, the 
spectrum of the PWN (and SN-remnant) above ~ 50 keV has to 
bend down. 

At GeV energies and in the TeV domain no (pulsed) source 
with a plausible assoc iation has been dete c ted in the vicin¬ 
ity of PSR J1811 -1925 dAharonian et al.ll2007l; iNolan et al.ll20lT 
lAbdo et al.l20l3l) . 

5.11 PSRJ1813-1246 

PSRJ1813-1246 wa s detected i n blind frequency searches using 
Fermi LAT data by lAbdo et Tl ( l2009dh . With a pulse period of 
48.1 ms and characteristic age of 43 kyr it is the second most en¬ 
ergetic {Lsd = 6.3 X 10^® erg/s) (radio-quiet) blind-search pulsar. 
It showed a double peaked 7 -ray (> 100 MeV) pulse profile with 
a peak separation of 0.485(3) and with bridge emission between 
the peaks. Two timing glitches have be en detected in the 7 -ray data 
dRav et al.l201 lbl;lMMelli et al.l2014bh . A plausible X-ray counter- 
part to PSR J 1813-1246 was found in a short Swift XKT observation 
dAbdo et al.ll2009dl;lRav et alj2011bl) . 

In a deep 108.9 ks XMM-Newton observation performed on 
March 10, 2013 with the EPIC-pn operated in small window mode 
highly significant pulsations were found with a lightcurve consist¬ 
ing of two peaks se parated ~ 0.5 in pha se lagging the 7 -ray one 
by ~ 0.25 in phase dMarelli et al.ll2014^ . The pulsed emission is 
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Figure 12. PSR J1813-1246; RXTE PCA 120 bins pulse profile for the 1.9- 
29.3 keV band. A double Gaussian plus background fit is superposed, yield¬ 
ing a 10.5(7 improvement with respect to a flat background description. 

very hard with a photon index of -0.85(3) and the pulsed fraction 
is very high, 96 ± 3%. We re-ana lysed the XMM-Newt on EPIC-pn 
data and confirmed the findings of iMarelli et al.l(l2014bll . The unab¬ 
sorbed 2-10 keV pulsed flux is (0.96 ± 0.015) x 10“^^ erg/cm^s. 
The 8-12 keV EPIC-pn pulse profile of PSR J1813-1246 is shown 
in Fig. Impanel r. We verified that very significant (12.1 ct) pulsed 
emission is detected above 10 keV. 

In the on-line version of the INTEGRAL IBIS 9-year 
Galactic hard X-ray survey iKrivonos et al.l 120121 see also 
http://hea. iki. rssi. ru/integral/nine-years-galactic-survey I the 35-80 
keV map of the Galactic bulge region (6.78 Ms exposure time) 
shows a 3.1a excess at the location of PSR J1813-1246 represent¬ 
ing a (total) flux of 0.43 ± 0.14 mCrab, which corresponds to 
(1.05 ± 0.34) X 10“® ph/cm^s keV. 

The hardness of the pulsed spectrum in the XMM-Newton 
bandpass and the detection of a point-source in the 35-80 keV IN¬ 
TEGRAL ISGRI band makes it very plausible that PSR J1813-1246 
shows pulsed emission at detectable levels in the hard X-ray/soft 7- 
ray band. In the X-ray archives we found a 26.8 ks RXTE PCA ob¬ 
servation (obs. id. 20090; taken in EN25usJ>4MJ)Ns event mode 
with a ~ 122.07/rs time resolution, and all PCU’s operational) 
performed on Nov. 7, 1997 targeting at GRO J1814-12, which has 
PSR 11813-1246 in the field of view at an off-axis angle of 23(7. 
Also, the much stronger LMXB 4U 1812-12 was in the field of 
view at an off-axis angle of 24(3 yielding increased (non-Gaussian) 
background levels in timing searches of PSR 11813-1246 in PCA 
data. In spite of the more than 11 years time lapse between the de¬ 
tection as (blind search) 7-ray pulsar and the RXTE observation a 
restricted period search {Z^ test, which turned out to be optimal 
for the EPIC-pn data, for each trial frequency) around the predicted 
pulse period adopting the ti ming parame t ers as given in supple¬ 
mentary material provided bv lAbdo et alj (1201 3h yielded a highly 
significant pulsed signal, 8.5 (j single trial, at 20.8047193(8) Hz 
(epoch MID 50759) for the 1.9-29.3 keV band. Taken into account 
the number of scan steps of ~ 830 independent Fourier steps the 
detected signal remains still very significant. 


The 1.9-29.3 keV PCA pulse profile with superposed the best 
fit double Gaussian model is shown in Fig. The phase separation 
between the peaks (peak 2 is the most dominant pulse at phase ~ 
0.82 in Fig. ll2b A<l>2-i is 0.509±0.004, consistent with theXMM- 
Newton value of 0.5043 ± 0.0009 (see lMarelli et al.l[2014hl) . The 
X-ray pulses, however, are much narrower, ai = 0.019 ± 0.005 
and (72 = 0.013 ± 0.002, than observed by XMM EPIC-pn with 
(71 = 0.0330 ± 0.0008 and (72 = 0.0305 ± 0.0005, because of the 
much better time resolution of the PCA observation with repect to 
the, XMM-Newton observation, i.e. ~ 122.07 /is versus ~ 5.7 ms. 

For the differential PCA energy bands, 3.9-7.9, 7.9-15.1 and 
15.1-33.5 keV, we obtained the following .Z|-test significances, 
5.3(7, 4.3(7 and 3.3(7, respectively, proving PSR 11813-1246 to be 
a soft 7-ray pulsar. Unfortunately, the statistics are too poor to de¬ 
rive reliable pulsed flux estimates for these PCA energy bands. 

The spectrum of the pulsed emission of PSR 11813-1246 is 
shown in the right panel of Fig. (red-red orange). The X-ray 
data points are derived in this work from the EPIC-pn data (tem¬ 
plate fitting per energy band, followed by a spectral fit adopting a 
power-law photon flux model), wh ile the high-energy 7-ray spec¬ 
tral points and fit are adopted from lAbdo et alj l l201 3ll. This figure 
also includes the 35-80 keV total flux measurement of INTEGRAL 
ISGRI. The latter flux is a good estimate of the pulsed flux, because 
the pulsed fraction is consistent with 100%. 

No TeV counterpart has been detected so far for PSR 11813- 
1246. 


5.12 PSR J1813-1749: The pulsar associated with HESS 
J1813-178/SNR G12.82-0.02/1GR J18135-1751 

lAharonian et al.l 1 I2OO5I) announced the detection of eight unknown 
VHE sources during scans of the inner galactic plane with the 
H.E.S.S. telescope from May to luly 2004. One of these new 
sources is the compact TeV source HESS 11813-178 embedded 
in a highly obscured region, which lacked at the time of discov¬ 
ery plausible counterparts at other wavelength and was classified 
as a “dark” particle accelerator. Soon after its det ection as TeV 
sourc e a soft 7-ray counterpart, IGR 118135-1751 dUbertini et al.l 
l2005t) . of HESS 11813-178 was discovered. Moreover, radio obser¬ 
vations revealed the presence of non-thermal radio emission from 
a young shell-type SNR G12.82-0.02 with a diameter of 2(5 
dBrogan et al.ll200^ . These properties do conjecture that we are 
dealing with a PWN powered by a so far undetected energetic pul¬ 
sar. _ 

_ ^Follow-up X-ray o bservations with Chandra dHelfand et al.l 

l2007h and XMM-Newton dPunk et alj 20071) indeed revealed the pu¬ 
tative pulsar and its nebula, although pulsations still had to be de¬ 
tected, because the operation modes of the involved X-ray instru¬ 
ments had insufficient time resolution to detect pulsations at time 
scales of 10-100 milli-second. 

A 30 ks Chandra ACIS-I image of the near environment of 
G12.82-0.02 is presented in Fig.|13l and shows the extraction re¬ 
gions used in iHelfand et alj d2007ll to derive the spectra of the (pu¬ 
tative) pulsar (circular aperture of 2" in radius centered on the point 
source), the inner pulsar wind nebula (6” x 8” elliptical region ex¬ 
cluding the point source region) and the outer pulsar wind nebula 
encompassing the bulk of the nebula extent (a 80" radius circular 
region offset from the point source excluding the point source and 
the inner nebula) . Fitting absorbed p ower-law models to the vari¬ 
ous components dHelfand et alj|200'% yielded for the outer PWN a 
Hydrogen column density Nh of (9.8 ± 1.2) x 10^^ c m~^ (90% 
confidence), consistent with the XMM-Newton findings dFunk et al.l 





































































































The soft 'j-ray pulsar population 19 



Figure 13. Chandra ACIS-I image (30 ks) of G12.82-0.02 showing clearly 
the point-source, the (putative) pulsar PSR J1813-1749 (hlack 2" circle), 
the inner PWN confined within a 6" X 8 " (yellow) elliptical region and the 
outer PWN confined to a 80” circular region offset from the point source 
(large white circle). 

l2007h . indicating highly absorbed emission. Both the pulsar and the 
outer PWN show hard power-law indices of about —1.3, while the 
inner nebula has even a harder index of —0.4 (see e.g. the spectral 
compilation in Fig. [H for the fits of the inner PWN, total pulsar 
emission and its combination as dotted, solid and dashed orange 
lines, respectively), a clear manifestion of Synchrotron cooling ef¬ 
fects. The total flux in the nebula is about 4.3 times larger than that 
measured for the point-source. 

Targeting at _ &e central point source in G12.82-0.02 

iGotthelf et all ( l2009l) finally discovered 44 .7 ms pulsations 
analysing a 98 ks XMM-Newton observation performed on March 
27, 2009, with the EPIC pn operating in small-window mode pro¬ 
viding a time resolution sufficient to detect timing signals up to ~ 
125 Hz. Follow-up observations with Chandra (ACIS-S3 Contin¬ 
uous Clocking; Feb. 12, 2012) and XMM-Newton EPIC pn (small- 
window mode; March 13, 2011) made it possible to derive a reliable 
period derivative estimate of PSR J1813-1749, yielding a charac- 
teristic age of 5.6 ky r and a spin-down power of 5.6 x 10^^ erg/s 
jHalnem et alj2oT% . 

We re-analyzed the EPIC pn data of the 98 ks XMM-Newton 
observation (obs. id. 0552790101), performed in small-window 
mode, with as main goal to derive the pulsed spectrum (not pub¬ 
lished so far) and from this the pulsed fracti on. The pulsed signa l 
was clearly visible at the predicted frequency jCotthelf et al.l20()^ 
and the pulse profile for the 2-10 keV energy band is shown in 
the upper panel of Fig. M Below ~ 2 keV no pulsed signal can 
be detected. The 21.la 2-10 pulse profile was fitted with a trun¬ 
cated Fourier-series using the fundamental and one harmonic, and 
the model has subsequently been used in the extraction process of 
pulsed excess counts. These pulsed excess counts have been con¬ 
verted to pulsed flux values adopting an absorbed power-law model 
in a forward folding fit procedure using response information tak¬ 
ing into account the 15” extraction radius. 

For comparison purposes with the spectral results determined 
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Figure 14. PSRJ1813-1749; the upper panel shows the pulse profile for 
the 2-10 keV band as measured by XMM-Newton EPIC pn during a 98 
ks obseivation performed in small-window mode. In the lower panel the 
RXTE PCA lightcurve of PSR 11813-1749 is shown for the 8.1-27.5 keV 
band (Z^ = 4.8 it) obtained from dedicated HESS 11813-178 observations 
performed about 500 days before the discovery of its pulsed nature. 

bv lHelfand et al.l ( l2007h for the total point source emission we fixed 
the Hydrogen column density Nh to 9.8 x 10^^ cm~^. We derived 
a photon index of —1.30 ± 0.03 and an unabsorbed/absorbed 2-10 
keV pulsed flux of (9.2 ± 0.45) x 10"^®/(5.8 ± 0.82) x 10"^^ 
erg/cm^s. This translates, using the absorbe d total 2-10 keV flu x 
value of 1.3 x lO”*^^ erg/cm^s as derived bv lHelfand et al.l(l2007h . 
to a (2-10 keV) pulsed fraction of 0 .4 5 ± 0.06, consistent with 
that estimated bv iGotthelf et al.l ( l2009t) ; [HalDern et alJ ( l2012h . The 
EPIC pn pulsed flux measurements and its best fit are shown in the 
spectral compilation depicted in Fig. M as dark-orange/red data- 
points and solid line, respectively. 

We also performed a spatial analysis of the EPIC pn small- 
window mode data by fitting (adopting Poissonian statistics) a 
model composed of the EPIC pn PSF, centered on the pulsar’s 
X-ray counterpart, and a (locally) flat background model to the 
measured 2d-event distribution (extraction radius 30”) for every 
required energy band. Because the EPIC pn PSF is worse than that 
of Chandra ACIS-I the extracted total counts are composed of the 
total emission component of the pulsar and a contribution from the 
inner pulsar wind nebula, which can not be resolved by XMM. 

Initially, we fitted the measured count spectrum with an ab¬ 
sorbed power-law model with Nh fixed to 9.8 x 10^^ cm~^. How¬ 
ever, this rendered a rather poor fit (Xr,i 7-2 = 28.90/15 = 1.93; 
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Figure 15. PSRJ1813-1749; periodogram showing versus trial fre- 
quency in a 50IFS window around the predicted frequency i/g using events 
from PCU-2 with PHA in the range 5-50 and employing all detection lay¬ 
ers. A 4.6 (t single trial peak maximum is found 6.8 IFS shifted from the 
prediction, resulting in an overall detection significance of 3.7 (t. The dot¬ 
ted line indicates the 3cr confidence level for a single trial. 


1.5% probability that a good model gives such a x^)- The flux 
measurements and its best fit, adopting Nh = 9.8 x 10^^ cm“^, 
are shown in Fig. [T^ as red datapoints and solid line, respectively. 
We also left Nh free in the absorbed power-law fit and obtained 
a good fit with Xr,i7-3 = 18.34/14 = 1.31 yielding a Nh 
of (11.7 ± 0.35) X 10^^ cm~ ^, slightly hig her, but still consis¬ 
tent with the value derived by iHelfand et alj l l2007tl using Chan¬ 
dra ACIS-I data. The unabsorbed / absorbed 2-10 keV flux was 
(3.30 ± 0.05) X 10"^V(l-94 ± 0.11) x 10"^^ erg/cm^s and the 
photon-index —1.31 ± 0.01. The absorbed 2-10 keV flux value is 
comparable with the sum of the point s ource - and inner nebu la flux 
of ~ 1.7 X 10“^^ erg/cm^s as given in iHelfand et al.l j2007h . 

We also revisited PCA/HEXTE data from a dedicated RXTE 
observation of HES SJ1813-178 (observation identifier 93022) per¬ 
formed during Nov. 16-20, 2007 (about 500 days before the de¬ 
tection of the pulsed nature) with only one operational PCU and 
at a 28!3 offset angle from PSR J1813-1749, thereby considerably 
reducing the PCA overall sensitivity. T he availability of an (inco¬ 
herent) ephemeris dHalnern et al .l2012h with a sufficiently accurate 
period derivative made a restricted period search around the pre¬ 
dicted (backwards extrapolated) pulse frequency value feasible. 

Selecting PCA events with PHA values in the range 5-50 (~ 
2 — 20 keV) using all detection layers from PCU-2, the only operat¬ 
ing unit during the 4-days covering observation (MJD 54420.101- 
54424.332) of about 125 ks, this search yielded a 4.6cr single trial 
peak shifted 6.8 Independent Fourier Steps (IFS= l/r, in which r 
is the full time stretch of the observation) from the prediction. Tak¬ 
ing into account the number of search steps (50 IFS) this results in 
an overall detection significance of 3.7a. The periodogram cover¬ 
ing 50 IFS and showing Zf versus trial frequency near the predicted 
frequency, is depicted in Fig. [B] We oversampled, for display pur¬ 
poses, one IFS by a factor of 20 to follow the peak structure in 
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Figure 16. PSRJ1813-1749; unabsorbed high-energy spectrum of the 
pulsed emission along with that of its inner PWN as obtained by CXO 
ACIS-I (2-10 keV; inner PWN, total pulsar and its combined emission; 
iHelfand et alJbOOTl) . XMM-Newton EPIC pn (2-10 keV; pulsed, total pul¬ 
sar plus inner PWN; this work), RXTE PCA (2-30 keV; pulsed; this work) 
and INTEGRAL ISGRI (20-300 keV; total=total pulsar plus inner PWN 
and whole PWN; this work). An absorbing Hydrogen column density of 
9.8 X 10^^ cm“^ has been assumed for all X-ray measurements to facil¬ 
itate the comparison of CXO ACIS-I and XMM-Newton EPIC pn spectral 
results. 


detail. This detection of pulsed emission by the PCA at a frequency 
shifted marginall}0 from the prediction is very plausible, because 
(W3)PIMM4i| predicts for the PCA with one operational PCU at 
an offset of about 30', reducing the source count rate by ~ 50%, 
an observation time of 111 ks to reach a 5a signal in the 0-49 PHA 
range given the 2-10 keV pulsed flux measured by XMM EPIC pn. 

We derived pulsed excess counts for three broad PCA energy 
bands (see e.g. panel b of Fig[T4]for the 8.1-27.5 keV lightcurve) by 
(XMM EPIC pn) template fitting and converted these to flux values 
(see Fig. [H aqua data points) using proper response information 
and adopting an absorbed power-law model with Nh fixed to 9.8 x 
10^^ cm“^. In HEXTE (> 15 keV) data the pulsed emission could 
not be detected. 

In the soft 7-ray band we obtained spectral information on 
the total emission of PSR J1813-1749 from the INTEGRAL IS¬ 
GRI skymaps produced in seven differential energy bands. These 
skymaps were centered on PSR J1811-1925 (see Section [S.lOl l. 
which is only 1?7 away from PSR J1813-1749. The flux measure¬ 
ments are shown as purple data points in Fig. M Fitting a power- 
law model we derived a photon index of —1.48±0.08 and a 20-100 


® The incoherent ephemeris given bv iHalpem et all i2Q12l) does not contain 
a i> value, and young pulsars often show glitches making precise predictions 
impossible given the 500 days back extrapolation. 

^ Portable, Interactive Multi-Mission Simulator; a HEASARC tool to 
convert count rates/fluxes observed for one high-energy instrument 
into that for another one, adopting a (photon) spectral model and 
absorbing Hydrogen column. Online documentation can be found at 
http://heasarc.gsfc.nasa.gov/docs/software/tools/pimms.html 
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keV flux of (2. 30 zb 0.13) x 10“^^ erg/cm^s, consistent with the 
value given bv lUbertini et alj ( l2Q05h . Note, that the ISGRI fluxes, 
due to the lack of resolving power at arcsecond scales, represent 
the combined emission from the pulsar (pulsed and unpulsed), the 
inner PWN and outer PWN. 

In the hard 7-ray band no pulsed emission from PSR J1813- 
1749 has been detected by the Fermi LAT instrument jAbdo et^ 

I 2 OI 3 I) . 


Finally, in the radio-b and no pulsed emiss ion has been de¬ 
tected so far at 1374 MHz jHelfand e t alj 200 7h using the ATNF 
Parkes telescope and at 2 GHz iHalpern et ^ 12012) using the 
NRAO Green Bank Telescope (GBT) telescope, rendering one of 
the deepest flux upper limits obtained so far for a young radio pul¬ 
sar (< 0.006 mjy at 2 GHz, equivalent to < 0.01 mjy at 1.4 GHz). 


5.13 PSR J1838-0655 / AX J1838.0-0655 

The detection of soft 7-ray emission up to ~ 300 keV from the 
ASCA source AX J1838.0-0655 has been reported bv iMalizia et al.l 
( I 2 OO 5 I) using INTEGRAL ISGRI data. Its location made an asso¬ 
ciation with TeV source HESS 11837-069 jAharonian et alj[2005ll 
plausible, an d thus suggests a puls ar/PWN origin. This was indeed 
confirmed bv iGotthelf et al.l 1 1200 8) . who detected a young (23 kyr) 
70.5 ms pulsar in the INTEGRAL ISGRI error circle using RXTE 
P CA data. Its sp i n-dow n rate, energetics and spectra were reported 
by iKuiper et~^ (l2008h . So far, no radio emission associated with 
th e pulsar in AX 1183 8.0-0655 has been reported (see e.g. Sect. 5 
of iMalizia et al.ll200l . who report a flux limit of 1-2 mly at 1.4 
GHz). 

Since its identification as a pulsar early 2008 with RXTE 
AX 11838.0-0655 has been monitored up to December 6, 2010. 
We have made accurate phase-coherent timing models (see Table 
[T]( for the period February 17, 2008 - December 6, 2010, during 
which a large timing glitch has been detected, occuring somewhere 
between MID 55002 and MID 55018 (lune 20 - luly 6, 2009) with 
a fractional frequency (Ai z/i/) jump size of (1.55 ± 0.07) x 10~® 
( iKuiper & Herrns^l2010h . 

Applying these ephemerides in a timing study of RXTE PCA 
and HEXTE data collected during observation programs P93429 
and P94303 (February 17, 2008 - December 5, 2009; MID 54513 
- 55170) resulted in the pulse phase diagrams shown in Fig. 117b- 
f. The screened (and PCU-detector break-down cleaned) exposure 
times for units 0-4 are: 60.34, 69.24, 256.60, 45.45 and 72.21 
ks, respectively. For HEXTE the deadtime-corrected on-exposure 
times for the same period are 163.13 and 85.51 ks, for cluster- 
0 (staring on-source) and cluster-1 (rocking), respectively. Below 
~ 16 keV a structured broad pulse is apparent. 

Application of the same timing models in a timing analysis 
of ISGRI data, spanning INTEGRAL revolutions 592 - 910 (Au¬ 
gust 19, 2007 - March 29, 2010; MID 54331 - 552843; 3.3 Ms of 
screened vignetting corrected effective on-axis exposure), revealed 
a 12.1(7 pulsed signal in the 15-150 keV band. The significance of 
the pulsed signal above 50 keV is still 7.2(7 (see Fig.ll7h). 


® For the period August 19, 2007-February 17, 2008 (MJD 54331-54513) 
the first RXTE PCA ephemeris of AX J1838.0-0655 (see Table [T) is used 
for ISGRI data taken before the monitoring of AX J1838.0-0655 com¬ 
menced with RXTE (thus in backwards extrapolation). The pulsed signal 
(with proper alignment) could easily be recognized in the ISGRI data dur¬ 
ing this period not covered by the validity interval of the first ephemeris. 
This proves that the ephemeris is also valid across a larger (backwards ex¬ 
tended) time period. 
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Figure 17. AX J1838.0-0655; pulse profiles in various energy bands for 
RXTE PCA (panels a-d), RXTE HEXTE (panels e-f) and INTEGRAL IS¬ 
GRI (panels g-h). Pulsed emission has been detected up to ~ 150 keV. 


From the RXTE PCA/HEXTE and INTEGRAL ISGRI pulse- 
phase histograms, we have extracted the pulsed excess counts 
through template fitting, and converted these to photon-flux values 
adopting (for the PCA range) an absorbing Hydrogen column Nn 
of 5.4 X 10^^ cm“^ (see e.g. lAnada et al.fl2oSl ; iKargaltsev et al.l 
l2012h . These pulsed flux measurements are shown in Fig.[T^ The 
combined RXTE PCA/HEXTE and INTEGRAL ISGRI pulsed flux 
measurements could be accurately fitted across the 3-150 keV band 
by a “curved” power-law model (a parabola in a (log(izFV) vs. 
log(iz)) spectral representation) with a photon-index of 1.36(2) at 
14.2 keV (this model is indicated by a black solid line in EigllSb. 
Pulsed flux estimates (unabsorbed) are (6.00 ± 0.33) x 10“^^ 
erg/cm^s for the 2-10 keV, and (2.53 ± 0.18) x 10“^^ erg/cm^s 
for the 20-100 keV band, respectively. 

Our derived value for the 2-10 keV unabsorbed flux for the 
pulsed emission is signi ficantly less than the value of 9.00 x 10“^^ 
erg/cm^ determined bv iGotthelf et al.l (l2008l) . who analysed much 
less RXTE PCA data. However, the spectral shape the latter authors 
derived for the PCA band is compatible with our “curved” spectral 
model in the overlapping region. 
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Figure 18. AX J1838.0-0655; high-energy (0.7-300 keV) total and pulsed 
spectra as derived from measurements by Suzaku XIS (total 0.7-10 keV; 
dashed brown), RXTE PCA (pulsed 2-30 keV; aqua), RXTE HEXTE (pulsed 
15-150 keV; blue) and INTEGRAL ISGRI (20-300 keV; total (purple) / 
pulsed (green)). The best fit model to the pulsed flux measurements is su¬ 
perposed as solid black line. The red dashed line and data points are flux 
estimates for the underlying PWN, which suggest a break/bend near 50 keV. 


We also derived total flux measurements of AX J1838.0-0655 
for energies above 20 keV from an INTEGRAL ISGRI imaging 
analysis of the Scutum region (AX JI 838.0-0655 was in the field 
of view of P SR JI846-0258 (see Sect. 15. 14t for which we reported 
the details in lKuiper & Hermseij ll2009l) ). combining observations 
performed between INTEGRAL revolutions 49 and 441 (March 10, 
2003 - May 27, 2006; c overing a considerably longer time span 
than iMalizia et al.l ( l2005h . who used data from Revs 46 to 186). 
Note that for this period no valid ephemerides exist and thus tim¬ 
ing results are lacking. The ISGRI total flux measurements (see 
FigUHfor data points and best fit) can be described by a power-law 
with a photon-index of — 1.72 ± 0.07, consiste nt with the value of 
—1.66 ± 0.23 obtained bv iMalizia et al.l ( l2005h . and a 20-100 keV 
(unabsorbed) energy flux of (3.69 ± 0.21) x 10“^^^ erg/cm^s. The 
20-300 keV flux is (7.35 ± 0.39) x 10"“ erg/cm^s, slightly les s 
than the value of 9 x 10“^^ erg/cm^s from IMalizia et al.l ( l2005l) . 
The pulsed fraction in the 20-100 keV band is 69 ± 6%. 

In the soft X-ray band (< 10 keV) spectral information 
for AX J1838.0-0655 and its PWN has been derived from Chan¬ 
dra ACIS jGotthelf et allbOOSh . Suzaku XIS jAnada et al.ll2009l) 
and XMM-Newton EPIC MOS-2 jKargaltsev et al.ll2012l) data. The 
Suzaku XIS spectral analysis of the pulsar plus its surrounding 
PWN used an extraction radius of 3' encompassing ~ 90% of 
the source photons, while the XMM-Newton EPIC MOS-2 spec¬ 
tral analysis utilized a 40” extraction radius centered on the source. 
Both (unabsorbed) flux measurements in the 2-10 keV band, con¬ 
verted from the 0.7-10 keV flux value of (13.2^0 6) 10“'^^ 

erg/cm^s for the XIS and from the 1-11 keV flux of (12.4 ± 2) x 
10“^^ erg/cm^s for the EPIC MOS-2, are consistent within uncer¬ 
tainties, (10.7tQ 5) x lO”*^^ erg/cm^s and (9.7 ± 1.6) x lO”*^^ 
erg/cm^s, respectively. Moreover, the derived photon indices for the 
XIS and MOS-2 measurements of -1.27 ± 0.11 and -1.25lo;?4, 


respectively, and absorbing Hydrogen column densities of (5.4 ± 
0.5) X 10^^ cm“^ and (5.22 q g) x 10^^ cm“^, respectively, are 
consistent. 

In a 20 ks Chandra ACIS-I observation of HESS J1837- 
069, as reported by ICotthelf etaP (l2008l) . AXJ1838.0-0655 and 
its PWN were detected 5(25 off-axis (from 13 aim-point) degrad¬ 
ing considerably the image quality. These authors performed spa¬ 
tially resolved spectral analyses for the compact central source - 
the pulsar - selecting events from a 5” x 7” elliptical aperture cen¬ 
tered on the source peak, and the PWN selecting events from a 1' 
radius aperture centered at a location slightly offset from the com¬ 
pact source, while excluding a 7" x 9” elliptical aperture around 
the compact source. They obtained an extremely hard spectrum for 
the compact source with a photon index of —0.5 ± 0.2 and unab¬ 
sorbed 2-10 keV flux of 8.8 x 10“^^ erg/cm^s, and for the PWN 
a photon index of —1.6 ± 0.4 and unabsorbed 2-10 keV flux of 
1.0 X 10“^^ erg/cm^s, both absorbed through a Hydrogen col¬ 
umn of (4.5 ± 0.8) X 10^^ cm~^. While the combined compact 
source and PWN unabsorbed flux of 9.8 x 10~^^ erg/cm^s is con¬ 
sistent with that derived by both Suzaku XIS and XMM-Newton 
EPIC MOS-2, the spectral shape does not! The cause of this dis¬ 
crepancy is not yet clear. 

The spectral model of the total (pulsed plus D C) emission 
from AXH838.0-0655 and its PWN, as derived by lAnada et al.l 
l l2009l) . using Suzaku XIS data for energies below 10 keV (=con- 
sistent with XMM-Newton result) is shown in Pig[T8] as a dashed 
brown line. The pulsed fraction in the 2-10 keV band is 0.56±0.04, 
including the PWN contribution in the total emission. Thus the gen¬ 
uine pulsed fraction of AX J1838.0-0655 is larger than 56%. 

We can estimate the spectrum of the underlying DC- 
component (0.7-150 keV), which (mainly) originates from the 
PWN, by subtracting the flux of the pulsed component from that 
of the total emission. These DC-component flux estimates are also 
shown in Pig. [T^ as red line (< 10 keV) and data points (> 20 
keV), and suggest that the PWN spectrum breaks/bends near 50 
keV, providing clues on the magnetic field in the PWN (e.g. its av¬ 
eraged strength). 

At GeV energies (> 100 MeV) no pulsed 7-ray emission has 
been detected so far from AX Jl 838.0-0655 in Eermi LAT data 
dAbdo et al.ll2013h . lLande et al.l ( l2012l) . however, reported recently 
the detection of spatially extended emission in the 10-100 GeV 
band from a location positionally consistent with HESS Jl837-069. 


5.14 PSRJ1846-0258 

ICotthelf et al.l (I 2 OOOI) reported the discovery of PSR J1846-0258 in 
ASCA and RXTE data. It is a “slow” radio-quiet pulsar, P ~ 324 
ms, but it has the smallest characteristic age, r ~ 723 y, of all 
known pulsars. Its surface magnetic field strength of 4.9 x 10^^ 
G is above the quantum critical field strength of 4.413 x 10^^ 
G and this classifies the pulsar as a high-B- field pulsar. It is lo- 
cated at the centre of SN-remnant Kes 75 dHelfand et all I 2 OO 3 I) 
and shows up as a bright hard X-ray source surrounded by a dif¬ 
fuse PWN, also emitting hard X-rays. Prom RXTE monitoring data 
it w as found that PSR Jl84 6-0258 behaves as a very stable rota¬ 
tor dLivingstone et al.ir2006l) . With INTEGR AL point-source emis¬ 
sion has been detected up to ~ 200 keV dMcBride et all |200^ : 
[Kuiper & Hermsen||2009l). Pulsed emi ssion has been detected up to 
~ 150 keV I Kuiper & Hermsenl2009t) . Por details about the pulsed 
and tot al high-energy spectrum a cross the ~ 2-300 keV band we 
refe r tolKuiEgr^Hennsenl d2009[ ). 

iKumar & Safi-HarbI 1 2008l) reported the detection of a dra- 
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Phase 

Figure 19. PSR J1849-0001; RXTE PCA 30 bins pulse profile for the ~ 
2 — 28 keV band (PHA 5-65). Two cycli are shown for clarity. A truncated 
Fourier-series fit adopting 5 harmonics is superposed. 


malic brightening of the pulsar in Chandra observations of Kes 
75 during June, 7-12, 2006. The pulsar’s spectrum softened con- 
sider ahly from a power- law spectrum with index T ~ 1.32 to 
1.97. iGavriil et al.[ ( l2008h showed that during this radiative event, 
lasting for about 55 days, phase coherence in pulsar timing was 
lost due to an unprecedented increase of the timing noise. They 
also discovered five short magnetar-like hursts during the outburst. 
iKuiner & HermsenI j2009l) discovered that the onset of the radiative 
event was accompanied hy a strong glitch in the rotation behaviour 
of the pulsar probably triggering the sudden release of magnetic 
energy. Thus, this over many years very stably behaving source, 
both temporally and spectrally, exhibited suddenly magnetar-like 
behaviour during the outburst, after which it continued again as a 
young energetic rotation-powered pulsar. Therefore, this high-B- 
field pulsar PSR J1846-025 8 might play a crucial role in revealing 
the connection between rotation-powered and magnetically pow¬ 
ered pulsars (magnetars). 

So far, the coherent signal from PSR J1846-0258 remains un¬ 
detected at radio-frequencies with a 4 (t upper-limit at 1.95 GHz 
ranging between 4.9 and 43 )rJy, depending on the assumed duty 
cycle, and at 1400 MHz 27 a round the time of the X-ray 
bursts ( I Archibald et alj[2003) . Also, IParent et aP ( 1201 ih have not 
detected the pulsed signal of PSR J1846-025 8 for energies above 
100 MeV, analyzing about 20 months of Fermi LAT data. How¬ 
ever, the source HESS J1846-029 has been detected at TeV energies 
position ally consistent with Kes 75 with a (post-trial) significance 
of 8.3cr toiannati-Atai et ^|2008|) . Its TeV flux is at the level of 
~ 2% of that of the Crab nebula, and its intrinsic extension is com¬ 
patible with a point-like source. 


5.15 PSR J1849-0001 / IGR J18490-0000 

IGRJ18490-0000 was discovered bv iMolkov et al.l j2004h in a sur¬ 
vey of the Scutum region, followed somewhat later by the detection 
of a faint TeV source, HESS J184 9-000, coincident with the hard 
X-ray source dTerrier et al.l |2Q0^ . suggesting the source to be a 


PWN powered by an energetic puls ar. Swift-XKT ( Rodriguez et al.l 
l2008l 12.3 ks) and XMM-Newton dTerrier et al.ll2008L 10 ks)~^ 
servations, both performed early 2006, revealed the X-ray counter¬ 
part of IGR J18490-0000 at soft X-rays and spectral analysis in¬ 
dicated a(yery)hard_^ectrum from the somewhat blurred coun¬ 
terpart. iRatti et aTI d2010h nailed down the location of IGR J18490- 
0000 to sub-arcsecond level by analyzing a 1.2 ks Chandra HRC 
observation performed on February 16, 2006. No obvious coun- 
terparts have been found at lower-energies in the i ’ and Ks bands 
dRatti et al.ll2010h . Also, no radio counterpart could be detected 
in the 610 MHz band of the GMRT in a systematic follow-up 
on INTEGRAL sou rces, resulting in an upper-limit of 3.5 mJy 
dPandev et Ii]|2003) . 

The true nature of IGR J18490-0000 was found by 
iGotthelf et al.l d201 ll) . who detected, pointing at the Chandra loca¬ 
tion, a young and energetic 38.5 ms pulsar using 112.3 ks of RXTE 
PCA exposure time, collected during dedicated observations span¬ 
ning Nov. 25 - Dec. 15, 2010. Analysing also XMM-Newton EPIC 
pn and MOS data from an 11 ks observation performed on April 
3, 2006 these au thors estimated a puls ed fraction of ~ 25% for the 
2-10 keV band. IGotthelf et alj d201lh found also evidence for ex¬ 
tended emission around IGR 118490-0000 between 20" and 150" 
radius. The spectrum of this diffuse component was much softer, 
with photon index of —2.1 ±0.3, than that of the point-source emis¬ 
sion, with photon index of —1.1 ± 0.2, while it represents about 
25% of the point source emission. Given the large Hydrogen col- 
umn density Nh of (4.3 ± 0.6) x 10^^ cm~^ dCotthelf et alj20ll|) 
the diffuse emission could also be interpreted as a dust-scattered 
halo around a bright point source. 

We (re)analysed the RXTE PCA/HEXTE data, collected dur¬ 
ing the dedicated observations, to obtain a phase coherent tim¬ 
ing model and to derive the pulsed emission spectrum over the 
3 — 150 keV band. Through template cross-correlation we de¬ 
rived pulse arrival times (ToA’s) and from these we determined a 
phase coherent timing solution, which is shown in Table[T] We pro¬ 
duced pulse profiles of IGR J18490-0000 for all 255 PHA channels 
of the PCA, allowing contrary to IGotthelf et all d201lh all detec¬ 
tor layers, by phase folding (using the coherent timing model) the 
barycentered event arrival times that passed our selection criteria, 
and subsequently sort these on PHA. The pulse profile for the PHA 
band 5-65 in 30 bins (~ 2 - 27 keV; Zl = 2791.9 = 35.2cr) 
is shown in Fig. EH and consists of one broad (somewhat struc¬ 
tured) single pulse. A smoothed version of this profile (a truncated 
Fourier series adopting 5 harmonics; shown in Fig. EH as solid 
line) has been used in the extraction procedure of the pulsed ex¬ 
cess counts, because we do not see pulse morphology changes as 
a function of energy. The pulsed excess counts have subsequently 
been converted to pulsed fluxes in a forward folding procedure as¬ 
suming an underlying power-law model absorbed through a col- 


yielding as best fit parameters an unabsorbed 2-10 keV pulsed 
flux of (4.21 ± 0.08) X 10“^^ erg/cm“^s and a photon index of 
— 1.37 ± 0.01. While the latter value is consistent with that derived 
by IGotthelf eFZI ( 1201 ll) . our derived flux value is about 4 times 
larger than the value quoted bv iGotthelf et alj (1201 ih . 

The HEXTE data from the (only active) staring cluster-A de¬ 
tectors have been screened adopting default screening conditions. 
The barycentered events have subsequently been folded using the 
coherent timing model to yield pulse profiles for all HEXTE PHA 
channels. Pulsed emission has been detected significantly up to 
- 60 keV (HEXTE band 31.0-60.1 keV; 4.8cr). Pulsed excess 
counts in various energy bands have been determined through tern- 
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Figure 20. PSR J1849-0001; the left panel shows the (15 bins) HRC-S pulse profile (0.06-10 keV), selecting events within a 2" aperture around its centroid. 
The background level at I.l counts/bin and its Itr uncertainty levels are indicated b y the horizontal solid an d dashed lines, respectively. The unpulsed (DC) 
level and its Itr uncertainty levels, as derived from the bootstrap method outlined bv ISwanenoel et al] J 19961) . are represented by the long dashed and dotted 
lines, respectively. The solid horizontal line segment given at the top indicates the unpulsed “region” as estimated by this bootstrap method. The intrinsic 
pulsed fraction is 0.77 ± 0.04. The right panel shows the radial distribution of HRC-S events centered on IGRJ18490-0000 up to 30” along with a model 
(dotted line) composed of a flat background (dashed horizontal line) and a point source with 406 counts. No indication for structured extended emission is 
found in the near vicinity of IGRJ 18490-0000. 


plate fitting. These counts have been converted to flux values using 
proper HEXTE cluster-A response information and the deadtime 
corrected exposure for cluster-A, adopting the photon power-law 
model derived from PCA data. The PCA and HEXTE pulsed flux 
values are shown in the spectral compilation depicted in Fig.l22las 
aqua (PCA) and blue (HEXTE) data points. 

Because our own analysis of the pulsed X-ray spectrum of 
IGRJ18490-0000 using RXTE PCA data yie l ded a ~ 4x larger 
pulsed flux than that given by iGotthelf et al.l )|201 ih . we obtained 
also a higher pulsed fraction of 88 ± 8 %. To shed light on this 
discrepancy and on the extended emission around IGR J18490- 
0000 we proposed 2 Chandra observations of 25 ks each, one with 
the HRC-SI (timing) and one with the ACIS-S (spatially resolved 
spectral analysis). Our aim was to determine the genuine intrinsic 
pulsed fraction in the X-ray band and to perform accurate spatially 
resolved spectral analysis at arcsec scales. 

The HRC-S observation in timing mode (6' x 30'field of view; 
time resolution 16/rs) was performed on Nov. 20, 2011 for an ef¬ 
fective exposure time of 25.1 ks (CXO observation id. 13292). An 
imaging analysis using Maximum Likelihood (ML) techniques for 
Poissonian distributed pixel cont ent (see e.g. the Supplementary 
Material of iHermsen et al.l 1201 3l . for a similar approach) of the 
HRC-S data showed a clear point-source at ( 02000 , <520oo) =(18 
49 L632,-00 01 17'.'45) with a formal systematic uncertainty of 
(/.'6 (in radius; 90%) and a negligible statistical uncertainty of O'.'02 
(Icr) in each coordinate. Selecting (barycentered) events within 2'' 
from the centroid location of IGR 118490-0000 yielded a number 
of 406 events of which 16.4 ± 1.2 events are expected to be back¬ 
ground events as estimated in an annulus of 4'' inner radius and 8'' 
outer radius centered on IGR 118490-0000. 

A restricted period search around the predicted value using 


the Z|-test statistic yielded a very significant signal at t^max = 
25.9609614(25) Hz, exactly at the extrapolated value from the 
(RXTE PCA based) coherent ephemeris (this work) as given in Ta¬ 
ble [T] Subsequent phase folding yielded the CXO HRC-S pulse 
profile (energy integrated; 15 bins) of IGR 118490-0000 as shown 
in the left panel of Fig. The pulsed fraction, not corrected for 
(the small) background, is 0.74± 0.04, which translates in a (back¬ 
ground corrected) intrinsic pulsed fractio n of 0.77rfc0.04, inco nsis- 
tent with the value of ~ 0.25 derived by IGotthelf et all i201ll) . but 
consistent with the value of 0.88 ± 0.08 we derived earlier from 
our own RXTE PCA spectral analysis of the pulsed emission and 
the XMM-Newton total s ource flux of IGR 118490-0000 as reported 
bv iGotthelf et al.l ) l201ll) . 

To study the extended emission in the vicinity of IGR 118490- 
0000 we applied a (deep) Gaussian smoothing to the central part 
of the raw image centered on IGR 118490-0000 and displayed the 
image on a log scale to emphasize possibly extended emission. 
The resulting map (PHA integrated) did not show evidence for the 
presence of diffuse emission near IGR 118490-0000. It should be 
noted, however, that given the response of the HRC-S with max¬ 
imum sensitivity between 1-2 keV (strongly reduced beyond 2 
keV) this image mainly shows the spatial information at soft X-rays 
(0.1-2 keV). Moreover, the genuine underlying emission around 
IGR 118490-0000 is strongly suppressed by the highly absorbing 
intervening Hydrogen column. 

Analyzing the HRC-S data, the more quantitative ML-imaging 
approach of the near environment of IGR 118490-0000 (in this 
case a region with a radial extent of 10” centered on the centroid 
scanned at O'.'2 scale) also yielded no evidence for structured dif¬ 
fuse emission. This method assigns 406 ±21 counts to the point 
source at a background rate of 1.224 ± 0.066 counts/arcsec^. The 
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Figure 21. PSR J1849-0001; ACIS-S3 2-10 keV image of the near en¬ 
vironment of IGRJ18490-0000 (black circle). The image is Gaussian 
smoothed and scaled logarithmically to emphasize diffuse emission, which 
is clearly visible in a “rectangular” region of size 75" X 37'.'5 (small yel¬ 
low box) around IGR J18490-0000. Just outside the larger yellow box of 
size 150” X 65”, used for the determination of the X-ray spectrum of 
the extended emission, a cluster of three point sources are popping-up at 
about 60" from IGR J18490-0000 (reddish spot). The two circular aper¬ 
tures (white dashed circles) each with a 60” radius indicate sky regions 
used to estimate the background. The purple annulus with radii at 75” and 
150” represents the region selected in a study of the diffuse PWN emission 
around IGR J18490-0000 using XMM-Newton MOS 1 &2 data. 


projection of the 2d-imaging information into a radial event distri¬ 
bution (O'.'25 binsize) is shown in the right panel of Fig. 1201 with 
superposed the point source model for IGR J18490-0000 with 406 
source counts on top of the local background, estimated in an an¬ 
nulus with 10” inner- and 25” outer radius. The latter background 
value is fully consistent with that derived from the 2d-ML method 
fitting the background and source simultaneously. No significant 
deviations between measured distribution and model can be found, 
indicating that there is no evidence for structured extended emis¬ 
sion in the near environment of IGR J18490-0000. 

Our second Chandra observation (CXO observation id. 
13291) with IGR J18490-0000 at the aim point of the back- 
illuminated ACIS-S3 chip was performed on Nov. 16, 2012 for 25.6 
ks yielding an effective exposure of 22.7 ks. A Gaussian smoothed 
and logarithmically scaled (to emphasize diffuse emission) ACIS- 
S3 map for energies in the range 2-10 keV is shown in Fig. HU 
Diffuse emission is now clearly visible confined more or less in a 
“rectangular” region of size 75” x 37'.'5 (small yellow solid rectan¬ 
gle/box) around IGR J18490-0000. A similar map for events with 
energies in the range 0.5-2 keV does not show evidence for dif¬ 
fuse emission, consistent with the HRC-S findings, while the point 
source emission from IGR J18490-0000 is still detected. 

More quantitatively, we produced a Maximum Likelihood ra¬ 
tio map selecting events with energies within 0.5-10 keV for a 
2' X 2' region centered on IGR J18490-0000 to search for point 
sources in the near environment of IGR J18490-0000, polluting the 
determination of the genuine diffuse emission. This method as¬ 
signs 2359 ± 49 counts to the point source IGR J18490-0000 in 
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Figure 22. PSRJ1849-0001; unabsorbed (Nh = 4.5 x 10'^'^ cm”-^) to¬ 
tal and pulsed spectta from 0.3-150 keV combining data from RXTE PCA 
and HEXTE (pulsed), INTEGRAL ISGRI (total), Chandra ACIS-S (total), 
XMM-Newton EPIC pn (pulsed & total) and MOS (total). Total flux mea¬ 
surements have filled square symbols, while pulsed fluxes have filled cir¬ 
cles. A fit (“curved” power-law model) through all pulsed flux measure¬ 
ments is shown as black solid line. The emission of IGR J18490-0000 is 
highly pulsed with a pulsed fraction of about 80%. The spectrum of the 
diffuse emission is also shown below that of the pulsar (a) lime triangles 
ACIS-S; 150” X 65” extraction box minus 5” circular source region; 2) 
green triangles XMM-Newton EPIC MOS 1-h 2; annulus of 75" inner- and 
150” outer radius). 

this band, yielding a source count rate (7^5-iokeV 0.104 ± 0.002 
c/s. In addition, five sources are popping-up with detection signif¬ 
icances > 5 ct of which three (with relatively high significances in 
the range 6.7 — 15.2cr) cluster in a small region at about ~ 60” 
from IGR J18490-0000 (see FigHU just outside the large dashed 
yellow box). In the estimation of the spectrum of the diffuse emis¬ 
sion we chose an extraction region for the diffuse emission exclud¬ 
ing these sources. In fact, these sources clearly pollu ted the diffuse 
spectr um derived from XMM-Newton EPIC data by ICotthelf et al.l 
ll201 ll) . who used an annulus of 30" inner- and 150" outer radius 
ce ntered on IGR J184 90-0000 as extraction region (see also Fig. 2 
of iTerrier et al.ll2008h . 

We used two rectangular extraction regions encompassing 
IGR J18490-0000, excluding a circular region of 5” radial extend 
centered on IGR J18490-0000, to evaluate the spectrum of the dif¬ 
fuse emission: 1) a 75" x 37'.'5 box (small yellow box in Fig|21t 
and 2) a 150" x 65" box (large yellow box in Fig|21t. We also care¬ 
fully chose two 60" circular background regions, free of sources, 
in the vicinity of IGR J18490-0000 (see also FigHU the two white 
dashed circles). The background subtracted count spectrum (0.6- 
7 keV band), properly taking into account the area difference of 
the diffuse- and background extraction regions, has been fitted in a 
forward folding procedure adopting an absorbed power-law model 
with a fixed Hydrogen column density of Nh of 4.5 x 10^^ cm~'^ 
(see further this section) and proper response information. We 
found for the large extraction box a photon index of —1.18 ± 0.05 
and an unabsorbed 2-10 keV flux of (7.1 ± 0.5) x lO”*^® erg/cm'^s 
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(see Fig. for the ACIS-S PWN fluxes and fit; lime coloured), 
about 5-6 times smaller than that of the IGRJ18490-0000 point 
source. For the smaller extraction box we found a similar photon 
index of —1.13 ± 0.06 and a 1.6 times lower flux. Thus, our 
CXO ACIS-S derived spectr um of the diffuse em ission is much 
harder than that derived by ICotthelf et alj ( 1201 ih using XMM- 
Newton EPIC data, while our flux estimates are compatible given 
the different extraction regions and methods used. 

Next, we derived the total ACIS-S (point source) count spec¬ 
trum of ICR 118490-0000 by fitting simultaneously a point-source 
model and a flat background model to the spatial event distri¬ 
bution (using the ML-method adopting Poissonian statistics) for 
each chosen energy band between 0.3 and 7 keV. These spatial 
event distributions have been sampled in a 5” circular region cen¬ 
tered on ICR 118490-0000 in 0”5 x d!h bins. The resulting source 
counts per energy slice have subsequently been converted to pho¬ 
ton flux values using proper energy response information and an 
effective exposure value of 22.67 ks adopting an absorbed power- 
law model. A model fit with Nh free yielded a column density of 
(4.30 ± 0.16) X 10^^ cm~^ and a photon index of —1.08 ± 0.02. 
Fixing Nh to 4.5 x 10^^ cm“^ resulted in a slightly softer pho¬ 
ton index of —1.15 ± 0.02 and an unabsorbed 2-10 keV flux of 
(4.11 ± 0.09) X 10“^^ erg/cm^s. These flux measurements are 
shown in Fig.[^as magenta data points. 

We studied the effects of pile-up on the spectral param¬ 
eter estimation for this CXO ACIS-S observation, because the 
measured count rate (0.5-10 keV) per 3.2 s frame time is 
~ 0.33 cnts/frame which is sufficiently large for a moder¬ 
ate pile-up impact (see e.g. the middle panel of Fig. 3 of 
“The Chandra ABC Cuide to pileup” which is referenced from 
http://cxc. harvard. edu/ciao/ahelp/acis_pileup. html I. 

The pile-up fraction // (i.e. the fraction of frames that have 
detected events containing two or more events) is about 15%, while 
the total fraction ft of events lost, either through grade or energy 
migration, with respect to the expect rate is about 28%. To study the 
impact of pile-up on our spectral results in more detail we simulated 
spectra using webspec adopting different pile-up model parameters 
(no pile-up, pile-up with grade migration parameter a = (0, 0.5,1) 
for an input photon flux model with photon-index -1.08 and nor¬ 
malization at 1 keV of 4 X 10“"^ ph/cm^s keV absorbed by a Hy¬ 
drogen column of 4.5 x 10^^ cm~^. We found for all the three 
different pile-up models a photon flux reduction of about 22 % with 
respect to the unpiled case. The pile-up impact on the reconstructed 
photon index was marginal and its size is about 0 . 1 . 

XMM-Newton has observed ICR 118490-0000 for a second 
time on March 23-24, 2011 for about 54 ks. EPIC MOS operated 
in Prime Full Window mode using a medium filter for an effec¬ 
tive exposure of 52.635 ks and 52.664 ks for MOS 1 and MOS 2, 
respectively. EPIC pn operated in Small Window mode (medium 
filter) for an effective exposure of 37.76 ks allowing studies of the 
pulsed signal given the 5.7 ms time resolution in this mode. 
These XMM-Newton observations do not suffer from pile-up ef¬ 
fects, and allow us to reconstruct the unpiled total (MOS 1 & 2,pn) 
and pulsed (pn, only) spectrum of the source. 

To extract the total emission spectrum of ICR 118490-0000 
(MOS 1 & 2; pn) we applied a similar procedure as used for the 
ACIS-S, now sorting the events with energy between 0.3-10 keV 
and within 60" from the ICR 118490-0000 centroid in 2" x 2" 
spatial bins for each energy band. Eor MOS 1 & 2 we derived, 
adopting an absorbed power-law model, an Hydrogen column den¬ 
sity Nh of (4.51 ± 0 . 10 ) x 10 ^^ cm“^ and a photon index T 
of —1.13 ± 0.01 (Xr,35-3 = 1.20). For the pn we found for 


Nh a value of (4.54 ± 0.10) x 10^^ cm ^ and for the pho¬ 
ton index —1.23 ± 0.01, somewhat softer than determined from 
MOS 1 & 2 data {xl.n-s. = 0.86 ; 1.42-10 keV). The derived 
Hydrogen col u mn de nsity is consistent with the value found by 
ICotthelf et al.l (1201 ih . and in all further spectral analysis of X- 
ray data we fixed Nh to 4.5 x 10^^ cm“^. The unabsorbed 2-10 
keV total flux of ICR 118490-0000 is (4.72 ± 0.05) x 10”^^ and 
(4.92 ± 0.05) X 10“^^ erg/cm^ s, for MOS 1 & 2 and pn, respec¬ 
tively, about 17% higher than the ACIS-S estimate, which suffered 
from pile-up effects. The flux measurements are shown in Fig. 1221 
as red (EPIC-pn) and orange (MOS 1 & 2) datapoints. 

We also studied the diffuse emission around ICR J18490-0000 
using MOS 1 & 2 data. We chose a source region consisting of an 
annulus with an inner radius of 75" (just outside the regions stud¬ 
ied using ACIS-S data) and an outer radius of 150" (purple annulus 
in Fig. l21t . well outside the “source cluster” at an angular distance 
of about 60" from ICR J18490-0000. The background region con¬ 
sisted of an annular region with 150" inner and 225" outer radius. 
The difference count spectrum, taking into account the difference 
in area of the source- and background extraction regions, has been 
fitted with an absorbed power-law model (Nh fixed to a value of 
4.5 X 10^^ cm“^) utilizing response functions produced for the 
annular source region. The fit yielded an unabsorbed 2-10 keV en¬ 
ergy flux of (6.23 ± 0.42) x lO”*^® erg/cm^s and a photon index of 
—1.75±0.05, considerably softer than found by the ACIS-S for the 
inner PWN region. This can best be explained as a manifestation of 
Synchrotron cooling, which makes the outskirts of a PWN softer 
than the inner parts. The MOS 1 & 2 PWN flux measurements are 
shown in Eig.j^as green datapoints. 

We also performed a timing analysis of the EPIC-pn data, 
which, sampled in Small Window Mode, has a time resolution of 
~ 5.7 ms, amply sufficient to study the 38.5 ms pulsed signal. Eor 
this purpose we selected barycentered events within 15" from the 
ICR 118490-0000 centroid. A restricted period search around the 
predicted pulsed period extrapolating the timing parameters of the 
ICR 118490-0000 entry listed in Table[T]yielded a very significant 
peak in the periodogram exactly atop the prediction. Subsequent 
phase folding adopting again the timing model listed in Table [T] 
for ICR 118490-0000 yielded significant pulsed emission down to 
about 1 keV. 

To derive the pulsed spectrum we first applied for each se¬ 
lected energy band, adopting Poissonian statistics, a Maximum 
Likelihood template fit procedure consisting of optimizing the 
number of pulsed excess counts and unpulsed (flat) level simul¬ 
taneously assuming that the pulse shape is described by the pulse 
profile (=template) as shown in Fig.ll9l 

These excess counts are converted to flux values in a forward 
folding procedure adopting an absorbed power-law model (Nh = 
4.5 X 10^^ cm~^) using instrument response functions determined 
for the 15" extraction radius. We derived an unabsorbed 2-10 keV 
pulsed flux of (3.65 ± 0.06) x erg/cm~^s and a photon 

index of —1.14 ± 0.01. The pulsed flux measurements by EPIC- 
pn are shown as dark-red datapoints in Fig.|^ Civen the total and 
pulsed 2-10 keV (XMM) fluxes we determined a pulsed fraction 
of 0.76 ± 0.02 for the 2-10 keV band, in full agreement with the 
HRC-S estimate for its integral energy band and the PCA based 
value. 

The sky region near ICR 118490-0000 was also frequently in 
the field of view of the soft 7 -ray imager ISCRI (15-300 keV) 
aboard the INTEGRAL satellite. A deep image centered on rotation- 
powered pulsar PSR11846-0258 (see Sectionl5.14l( and enclosin g 
ICR 118490-0000 has been shown by I Kuiper & Herm^ ( l2009h . 
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Figure 23. PSRJ1930+1852; pulse profiles in various energy bands for 
RXTE PCA (panels a-c) and RXTE HEXTE (panel d). Pulsed emission has 
been detected up to ~ 50 keV (4.5cr; panel d). 


The total spectrum of IGRJ18490-0000 for the 20-300 keV band 
has been derived from maps produced in 10 logarithmically binned 
energy bands covering the 20-300 keV window. The selected IN¬ 
TEGRAL observations covered INTEGRAL revolutions 49 - 603 
(March 10, 2003 - Sept. 23, 2007) and represent a total GTl ex- 
posure of 7.35 Ms . For m ore details we refer to Section 3.1 of 
iKuiner & HermsenI ( l2009h . Fitting a power-law model we derived 
a photon index of —1.82 ± 0.14 and a 20-100 keV energy flux of 
(2.25 ± 0.16) X 10~^ ^ erg/c m^s, consistent with the values deter¬ 
mined bv lTerrier et alJ l l2008h . who used less INTEGRAL exposure 
time. These ISGRl flux measurements (sum of pulsed/unpulsed 
emission from the pulsar and the diffuse emission) and best fit 
are shown in Fig. [22] as purple data points and solid line, respec¬ 
tively. Finally, in the GeV 7 -ray band no pulsed em ission from 
IGRJ18490-0000 has been detected by Eermi LAT jAbdo et al.l 

I 2 OI 3 I) . 

In summary, we found hard pulsed emission for IGR J18490- 
0000 from ~ 1 keV to ~ 60 keV, which represents about 80% 
of the total emission. The diffuse emission, 5-6 times weaker than 
the point source emission of IGR J18490-0000, is also hard and 
softens towards the outer parts of the PWN. This is characteristic 
for Synchrotron cooling. 


central bright hard point-source - the putative p ulsar - and several 

diffuse structures surrounding it i ILu et al.ll2002h. _ 

Targeting at the central source ICamilo et ^ (l2002h discov¬ 
ered the energetic pulsar PSRJ1930-1-1852 at radio frequencies. 
Its period of ~ 136 ms was also detect ed at soft X-rays in the 
archival ASCA GIS data. ILu et al.l (|2007|) showed, analysing ad¬ 
ditional Chandra AClS-1 data taken on June 30, 2003 for a total 
exposure time of 58.4 ks in CC-mode, that the pulsed emission is 
hard with a photon-index F of 1 . 2 ( 2 ) and that the source has a high 
pulsed fraction of 71 ± 5%. The absorbed pulsed flux in the 2-10 
keV band was 1.2 x 10“^^ erg/cm^s, which translates to an unab¬ 
sorbed flux of ~ 1.33 X 10~^^ erg/cm^s in the same band taking 
into account the 1.6 x 10^^ cm~^ absorbing Hydrogen column. 
These authors also analyzed data from a ~ 80 ks RXTE observa¬ 
tion taken during September 12-14 and December 23-25, 2002, 

but only for timing purposes. _ 

Guided by the available radio ephemeris dCamilo et al.ll2003) 
we generated a phase coherent ephemeris based on these X-ray data 
for the 2002 September 12 - December 25 period, listed in Table 
[T] We re-analysed the RXTE PCA and HEXTE data using template 
fitting to obtain the pulsed emission spectrum over a much broader 
bandpass than Chandra. Pulse profiles for various energy bands 
are shown in Eigj^and pulsed emission has been detected up to 
~ 50 keV (see Fig. 12311). The reconstructed photon spectrum of the 
pulsed emission across the 2.8-32.1 keV range, ad opting a Hydro - 
gen column density Nh of 1.6 x 10^^ cm“^ (see ILu et alj|2002h . 
has a power-law index of —1.30 ± 0.02 and an unabsorbed 2-10 
keV flux of (1.57 ± 0. 07) x erg/cm^s, slightly higher than 

the value obtained by ILu et alj (l2007l) who used Chandra ACIS 
data. Using Chandra ACIS-S data of observations performed over 
the period 8- 15 July 2008, totaling 290.77 ks of exposure time, 
iTemim et al. I ( I 2 OIOI) ' derived for the total pulsar spectrum a pho¬ 
ton index of —1.44 ± 0.04 and an unabsorbed 0.3-10 keV flux 
of 3.26 X 10“^^ erg/cm^s, which translates to a 2-10 keV flux of 
2.25 X 10“^^ erg/cm^s. From this total 2-10 keV flux value and 
our pulsed 2-10 keV flux we derive a pulsed fraction of ~ 70%, 
consistent with the value derived by using only Chan dra data. The 
pulse d {RXTE PCA and HEXTE) and total (C/iflnrfra: iTemim et al.l 
iTOldl emission spectra of PSRJ 1930-1-1852 are shown in Pig.[24l 
Soft 7 -ray emission in the 20-60 keV band from a location 
consistent with PSR J193 0-tl852 has also bee n detected in an imag¬ 
ing study of ISGRl data llSegreto et al.ll20l5l . 

PSR J1930-I-1852 is not listed a s high-energy 7 -ra y pulsar in 
the Fermi second pulsar catalogue (lAbdo et al.ll2013h . however, 
at TeV energies a pointlike source, VER 119304-188, has been 
detected by the VERITAS ground-based gamma-ray observatory 
from the direction of SNR G54.1-I-0.3 with an integral flux above 
1 TeV of 2.5% of the C rab nebula, very lik ely the PWN associated 
with PSR J1930-H1852 dAcciari et al.ll2010h . 


5.16 PSRJ1930+1852 


5.17 PSRJ2022-F3842 


ISeward et alj ( Il989h reported the detection of the Crab-like SNR 
G54.1-I-0.3 at X-rays using Einstein IPC (0.5-4 keV) data from an 
observation performed on May 7, 1980 collecting about 50 pho¬ 
tons, too faint to determine accurate spectral infor mation. More de- 
tailed X-ray spectral information was obtained by ILu et alJ l l200lh 
using data from the ROSAT PSPC (0.1-2.5 keV; - 20 ks) and the 
ASCA GIS (16.5 ks) and SIS (~ 20 ks) detectors, both operating 
in the 0.5-10 keV range. Spatially resolved spectral analysis in the 
X-ray band became possible using the superb imaging quality of 
the Chandra ACIS instrument, revealing in a 30.9 ks exposure a 


PSRJ2022-I-3842 was discovered at 1.95 GHz with the GET 
SPIGOT spectrometer targeting at the point source detected 
in the X-ray band at the center of SNR G76.9-t-1.0 (see e.g. 

1993) in a 54 ks Chandra 


IWendker et al.lll99l: Landecker et alj 


ACIS observation Jat: 


zoumanian 1 


et al.ll201lh . The radio pulsations 


of PSR J2022-I-3842 are severely affected by dispersion (DM = 
429.1 ±0.5) and scattering in the interstellar medium, and these 
effects explain the non-detections of pulsed emission in targeted 
searches at lower rad i o freq uencies, 618 and 1170 MHz, with the 
GMRT dMarthi et al.l 1201 ll) . Subsequent radio-observations with 
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E [MeV] 

Figure 24. PSR J1930+1852; high-energy (0.5-250 keV) total and pulsed 
spectra as derive d from measurement s by Chandra ACIS (total pulsar 0.5- 
10 keV; solid redi lTRinim i2O10l) l. RXTE PCA (pulsed 2.5-30 keV and 
best fit; aqua; this work) and RXTE HEXTE (pulsed 15-250 keV; blue; this 
work). 


the GBT of this weak pulsar (60 )rjy at 2 GHz with a very fiat 
spectrum, index a = —0.33) settled the timing properties of 
PSR J2022-I-3842 llArzoumanian et aT]l201 ll) . A spin-up glitch oc- 
curing between MJD 54400 and 54950 was also evident in these 
radio observations. 

Dedicated RXTE monitoring observations performed between 
Jan. 27, 2010 and Febr. 4, 2010 revealed the pulsed signal (a sin¬ 
gle sh arp pulse at pulse period 24.3 ms) also in the 2-20 keV PCA 
band l lArzoumanian et al]|201 ih . These authors found a hard spec¬ 
tral index of —1.1 ± 0.2 and an absorbed 2-10 keV pulsed flux of 
5.4 X 10“^® erg/cm^s, which is comparable to the total measured 
absorbed 2-10 keV flux by Chandra, indicating a pulsed fraction 
of ~ 100%. 

An XMM-Newton observation of PSR J2022-I-3842 of 116 ks 
performed on April 14, 2011 revealed that the pulse period is ac¬ 
tually twice (48.6 ms) the initially found period due to the pres¬ 
ence of an interpulse at about ^ 0.5 phase separation of the main 
pulse llArumugasamv et al.l2014l) . The adapted spin-down luminos¬ 
ity and characteristic age are now: Esd ~ 2.97 x 10®^ and 8.9 kyr, 
respectively. 

We re-analyzed the RXTE PCA data (88.73 ks of screened 
PCU-2 exposure) in order to deal with this “change” from single 
peaked to double peaked nature and to better characterize the hard 
X-ray spectrum by including all PCA detection layers, being more 
sensitive to photons with energies in excess of ~ 10 keV. We also 
analyzed the HEXTE data. A high-statistics template pulse profile 
(double peaked) was used in the correlation analysis to derive the 
pulse time-of-arrivals from which we constructed a phase coherent 
timing model (see TableO. This ephemeris has subsequently been 
used in a folding procedure creating the event distributions as a 
function of pulse phase and energy for both the PCA and HEXTE. 
Pulsed emission has been detected up to ~ 30 keV (see Fig.l25|pan- 
els a & b for the PCA 1.7-11.2 (28.2cr) and 11.2-27.3 keV (lO.Scr) 



Pulse Phase 

Figure 25. PSR J2022-I-3842; pulse profiles for RXTE PCA (panels a & b) 
and RXTE HEXTE (panel c). Pulsed emission has been detected up to ~ 30 
keV. Note the sharpness of the pulses of about 1.2 ms (EWHM). 

and panel c for the HEXTE 14.8-29 keV (3.Ter) lightcurves, respec¬ 
tively). 

The broad-band PCA pulse profile can satisfactorily be fitted 
by a model composed of two Gaussians plus background yield¬ 
ing a pulse phase separation of 0.484 ± 0.001 between main- and 
interpulse and pulse widths (FWHM) of 0.0254 ± 0.0014 (main 
pulse) and 0.0273 ± 0.0024 (interpulse). The X-ray pulse profile 
of PSR J2022-1-3842 bears st rong similarities wi th the pulse profile 
of PSR J0205-I-6449 (see e.g. lKuiner et'aDl2010l) . 

The pulsed excess counts for the main pulse (PI), interpulse 
(P2) and its sum (total pulsed; TP=Pl-l-P2) have been determined 
across the 2.5-32 keV PCA band (for which a proper energy cali¬ 
bration exists) through fitting the measured pulse profiles for vari¬ 
ous differential energy bands in terms of the previously derived best 
Gaussian models (fixed position and shape) and background. Sub¬ 
sequently, these excess counts have been converted to pulsed fluxes 
in a forward folding spectral fitting procedure assuming a power- 
law model and adopting an absorbing Hydrogen column density 
Nh of 1.7 X 10^^ cm"2 l lArzoumanian et ^1201 ll) . 

We found for the 2.5-32 keV band a photon index of —1.31 ± 
0.02, —1.04 ± 0.03 and —1.20 ± 0.014 for the main pulse, inter 
pulse and total pulse, respectively. The corresponding unabsorbed 
2-10 keV band fluxes are (3.60 ±0.12) x 10■^^ (1.91 ±0.11) x 
10““ and (5.50 ± 0.17) x 10““ erg/cm^s, for PI, P2 and TP, 
respectively. The main pulse thus comprises ~ 65% of the total 
pulsed flux. The absorbed/unabsorbed 2-10 keV total pulsed flux of 
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Figure 26. PSRJ2229+6114; pulse profiles in various energy bands for 
RXTE PCA (panels a-b), RXTE HEXTE (panels c). Pulsed emission has 
been detected up to 22 keV. 


(4.94±0.37) X 10 ~^^/(5.50±0.17) X 1 0~^^ erg/cm^s is somewhat 


less than found bvlAi 


rzoumanian 


et al.l ll 201 ih . but is still compara¬ 


ble (93 %) to the total absorbed 2-10 keV flux of the point source 
derived from Chandra (imaging) data. 

PSR J2022-I-3842 is not listed as high-energy 7 -ray pulsar in 
the Fermi second pulsar catalogue l|Abdo_etalj20]3^jhowever, the 
Fermi LAT second source catalogue ( Nolan et alj|2012h reports on 
a 10a 7 -ray source, 2FGL J2022.8-l-3843c, positionally consistent 
with G76.9-H1.0. At TeV energies (100 GeV - 50 TeV) the VERI¬ 
TAS survey of the Cygnus region, reaching a VHE sensitivity of 
~ 4% of the Crab nebula flux for an energy threshold of 200 
GeV, did not result in a detectio n of a TeV sourc e at the location 
of PSRJ2022-I-3842/ G76.6-tl.0 jWard et alj2010t) . 


5.18 PSRJ2229-t6114 

PSR J2229-t6114 was discovered at radio frequencies in the error 
region of the unidentified EGRET source 3EG J2227-I-612 after the 
end of the CGRO mission jHalpem et al.l20oT9) . The pulsed signal 
(P ~ 51.6 ms) was also detected in the X-ray band using archival 
ASCA GIS (0.7-10 keV) data, taken 1.6 year before the detection of 
the radio pulsations, showing 2 pulses separated about 0.5 in phase. 

A non-thermal incomplete radio shell, G106.6-I-2.9, surrounds 
PSR J2229-I-6114. This compact radio nebula of size ~ 3', called 
the “Boomerang”, has a flat spectrum and is polarized at a level of 
~ 25%, and represents the PWN associated with PSR J2229-I-6114. 
It is located at the northeastern tip of the radio continuum 
source G106.3-I-2.7 of size 0?5 x 1° which is class ified as the 
remnant of a supernova (see e.g. iKothes et alj 1200^ . that gave 
birth to PSRJ2229-I-6114. Diffuse emission from the near field 
of PSR J2229-I-6114 has also been clearly detected at X-rays by 


iNg & Romanil (|2004|) in a 95 ks Chandra ACIS-S observation 
(taken March 15, 2002), who fitted the observed morphology with 
their pulsar wind torus model. 

In recent years the pulsed fingerprint of PSR J2229-I-6114 has 
also securely bee n found in AGILE and Fermi LAT data at e nergies 
above 100 MeV jPellizzoni et al.ll2009tlAbdo et al.ll2009bl) . show¬ 
ing a single asymmetric peak at a pha se separation of 0.49 from the 
(single) radio pulse. At TeV energies [Acciari et al.l ( l2009ll reported 
the detection of extended emission from G106.3-I-2.7, as measured 
by VERITAS. 

We analyzed the - so far unpublished - RXTE PCA and 
HEXTE data from RXTE observation runs 60130 and 80092 (per¬ 
formed on Dec. 7-8, 2001 and Nov. 25-29, 2003 for 90.4 ks and 
129.8 ks, respectively) in order to characterize the pulsed emission 
spectrum of PSR J2229-t6114 over the 3-250 keV band. In both 
observation runs the pulsed signal could easily be detected in the 
PCA ~ 2-10 keV band (PHA range 4-27) taking into account the 
phase smearing due to the frequency derivative by including the 0 
value of —2.92006(5) x 10“^^ Hz/s, found at radio frequencies 
for a later epoch. The event cubes (2d distribution of events versus 
pulse phase (60 bins) and energy (PHA; 256 bins)) for both runs 
have been made and combined after applying an appropriate phase 
shift to align both distributions in phase. The combined (cleaned) 
exposure for PCU-2 amounts 220.144 ks. Similar folding and align¬ 
ment (with known phase shift) procedures have been followed for 
HEXTE data. The pulse profiles of PSR J2229-I-6114 obtained for 
the PCA bands 1.9-9.4 keV (PHA 5-22) and 9.4-22.4 keV (PHA 
23-53) are shown in Pig. 126b (12.3a) and b (8.3a), respectively. 
The double peaked morphol ogy, already visible i n the ASCA GIS 
data for the 0.8-10 keV band jHalnern et al.l2001bh . is evident, now 
also for energies above 10 keV. Pulsed emission has been detected 
up to ~ 22 keV. In HEXTE data we found only indications for the 
pulsed signal for energies above ~ 15 keV (see Pig. 126b; 2.4a). 

We also analyzed archival ASCA GIS dOhashi et al.ll99^ data 
obtained during an observation run (57038000) on August 4-7, 
1999 with an e ffective exposure of ab out 113.8 ks. We confirmed 
the findings of iHalpem et alj l l2001bh . but derived a more signifi¬ 
cant pulse profile for the GIS 0.8-10 keV range with a more pro¬ 
nounced double peak morphology by using 2 ' as optimum extrac¬ 
tion radius instead of 4'. We fitted this lightcurve with 2 symmetric 
Lorentzians on top of a flat background and derived a phase sep¬ 
arations of 0.49 ± 0.02. The individual Lorentzian profile shapes 
are subsequently used as templates in the derivation of the pulsed 
excess counts for both ASCA GIS and RXTE PCA pulse profiles 
for user defined energy bands. Applying proper response informa¬ 
tion, (effective) exposure times, correction factors for the missing 
fraction of the PSF (ASCA GIS; the extraction region was a circular 
aperture of 2', which covers about 57% of the PSF) these pulsed ex¬ 
cess counts have been converted to pulsed fluxes in a forward fold¬ 
ing spectral analysis procedure, adopting an absorbed power-law 
model assuming an absorbin g Hydrogen column Nh of 6.3 x 10^^ 
cm“^ dHalpern et al.l2001^ . 

For the PCA 2.8-32.2 keV band the pulsed emission (sum of 
the two pulses) is properly described by a power-law with (hard) 
photon index of —1.11 ± 0.03 and an unabsorbed 2-10 keV flux of 
(5.2±0.3) x 10~^^ erg/cm^s, consistent with the value derived by 
iHalnem et al.l (l2002h using ASCA GIS data. We found no evidence 
for different spectral behaviour of the pulses in the X-ray band. 
However, between ~ 20 keV and 100 MeV the pulse morphology 
drastically cha nges from a double peaked to an asymmetric single 
peaked profile dAbdo et al.l2009bl) . indicating different spectral be¬ 
haviour as a function of pulse phase. Unfortunately, there are no 
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sensitive measurements in this energy band to further constrain the 
pulsed emission properties. The total pulsed emission spectrum of 
PSR J2229+6114 from soft X-rays to high-energy 7 -rays is shown 
in Fig. |28] including ASCA GIS (0.8-10 keV), RXTE PCA (2.7-32 
keV) and Fermi LAT measurements as orange data points/symbols 
and hatched region. 


6 CANDIDATE SOFT 7 -RAY PULSARS 

The sample size of detected soft 7 -ray pulsars is still small, despite 
efforts to increase this sample by searching for timing signatures 
for energies above 20 keV from very promising candidate RPPs. 
Such strong candidates for detection at hard X-rays are Fermi LAT 
detected very young (< 15 kyr) energetic (Lad ^ 10^® erg/cm^ 
s) pulsars, often coincident with a TeV source, but also similarly 
young and energetic pulsars which are not detected by the LAT. In 
this section we will discuss the status (from literature and our work) 
for the candidate soft 7 -ray pulsars as defined above. 


6.1 PSRJ1833-1034inSNRG21.5-0.9 


An interesting example is the weak (51400 = 71/iJt/ ) radio-loud 
and Fermi LAT pulsar PSR J1833-1034 (r ~ 4.8 kyr; L^a ~ 3.4 x 
lO®’^ erg/s) at the centre of SNR G21.5-0.9, associated with TeV 
source HESS J1833-105 dPiannati-Atai et ^ 2008ll. whic h has a 
bright INTEGRAL soft 'v-rav counterpart dDe Rosa et al.l l2009h . 
Therefore, it evidently is a soft 7 -ray pulsar candidate. However, 
deep timing searches in NuSTA R data (3-79 keV) did not reveal a 
pulsed signal from this pulsar dNvnka et al.ll2014ll . This suggests 
that the soft 7 -ray beams are for this pulsar not pointed in the ob¬ 
server’s direction, contrary to the high-energy 7 -ray beams. 


6.2 PSR J1119-6127, PSR J1124-5916 and PSR J1357-6429 

Other puzzling examples are the three, very young, energetic pul¬ 
sars PSR HI 19-6127 (r ~ 1.6 kyr; L,a 2.3 x 10®® erg/s), 
PSR 11124-5916 (r ~ 2.8 kyr; L,a 1.2 x 10®^ erg/s) and 
PSR 11357-6429 (r ~ 7.3 kyr; L^a ^ 3.1 x 10®® erg/s), all ra¬ 
dio and Fermi LAT pulsars associated with SNRs. From these pul¬ 
sars, pulsed thermal X-ray emission has been detected (well below 
2 keV and with a broad single pulse), but not the expected strong 
non-thermal (> 2 keV) pulsed emission (s ee e.g. iGonzalez et al.l 
1 20051; iHughes et al.l2003l; IChang et alj2^. respectively). For two 
of the three pulsars also TeV counterparts have been identified, 
HESS 11119-614 for PSR 11119-6127 and HESS 11356-645 for 
PSR 11357-6429. 


6.3 PSR J1420-6048 and PSR J1418-6058 in the Kookaburra 

Two other potential s oft 7 -rav pulsars are located in the “Kook¬ 
aburra” radio complex dRoberts et al.ll99^ . coindicent with CGRO 
EGRET source 3EG 1420-6038 (and GeV 11417-6100). Wi t hin 
this complex, encompassing an area of about 1 ° x 1 °, two young 
and energetic pulsars have been detected at an angular separa¬ 
tion of only 16(5. One, PSR 11420-6048 (r ~ 13.5 kyr; L,a ^ 
1.0 X 10®^ erg/s) is locate d in the North-Eastern wing as radio pul¬ 
sar dp ’ Amico etjilj [200l[^ and is late r also detected as high-energy 
-Y-rav pulsar dWeltevrede et al.ll2010h . The other, PSR 11418-6058 
(r ~ 10.2 kyr; Lsa — 5.0 x 10®® erg/s) is in the South-Western 
wing, coincident with a bright radio nebula called the “rabbit”. 


and is detected in a blind search at high-energy 7 -rays dAbdo et^ 
l2009dl) . _ _ 

Moreover. lAharonian et alJd2006bll detected TeV counterparts 
for both pulsars HES S 11420-607 and HESS 11418-609, while 
iHoffmann et al.l d2007h reported on a soft 7 -rav source positio nally 
consistent with PSR 11420-6048 (see also lpiocchi et alfcoiol) . 

A (very) hard X-ray point source, embedded in a weak dif¬ 
fuse nebula, has been identified for PS R 11420-6048 wi thin the (ra¬ 
dio) K3 nebula dRoberts et al.ll999h bv iNg et ^ d2005ll who found 
for a 10 ks Chandra ACIS-S observation only 26 counts within 
an aperture of 1" around PSR 11420-6048. Further spectral anal¬ 
ysis yielded a photon index F of —1 with a very large uncer¬ 
tainty due to the poor source statistics. We analysed a 90.7 ks on- 
axis Chandra ACIS-I observation performed on December 8 , 2010 
(obs.id. 12545; unpublished so far) to determine the (total) spec¬ 
trum of the pulsar. Adopting an absorbed power-law model and 
using a Maximum Likelihood (energy resolved) PSF fit method 
to extract the source counts (now 233 counts for the 1.5-10 keV 
range), we derived a photon index F of —0.46 ± 0.07, an unab¬ 
sorbed/absorbed 2-10 keV (total) flux of (1.33 ± 0.10) ■ 10“®® 
/ (1.14 ± 0.30) • 10“®® erg/cm^s and a Hydrogen column den¬ 
sity Nh of ( 3 . 35 lo' 5 i) X 10^^ cm“^, which is slightly larger 
than the Galactic column of (1.58 — 2.13) x 10^^ cm“®. This 
confirms the hardness and weakness of the source with greatly 
improved accuracy. Our derived column density is also consis¬ 
tent with earlier X-ray derived values, using different instruments 
aboard XMM-Newton and Suzaku and adopting much larger ex¬ 
traction radii, an d thus including also the softer diffuse nebula 
emission (see e.g. IVan Etten & Romanill2010l ; iMarelli et al.ll201 J ; 

iKishishita et 

However, so far in the canonical X-ray band (0.3-10 keV) 
pulsed emis sion has not b een securely detected yet from PSR 
J1420-6048 e^^^J^ZOO^^^^n spile of some weak indications m 

ASCA GIS data teoberts et all 1200 ih . Given the hardness of the 
source spectrum and encouraging results from W3PIMMS con¬ 
cerning its detectability in RXTE PCA dat43 we searched for the 
non-thermal pulsed component of PSR J1420-6048 above 2 keV 
analysing a 66.4 ks on-axis RXTE PCA observation (obs.id. 80088) 
performed in the period August 4—6, 2003. The pulsar was observed 
in a post-glitch episode with v alid (radio) ephemerides before and 
after the glitch (I Yu et alj2013h . thereby considerably narrowing the 
period search window. However, also now for the 2-30 keV PCA 
band the non-thermal pulsed emission component could not be de¬ 
tected within the predicted period search window. 

The region around the other pulsar, PSR 11418-6058, within 
the “Rabbit” nebula has been resolved by Ch andra (obs.id. 2794; 

10.1 ks; September 22, 2002; 1^ et aljEoO^ into two unresolved 
point sources , denoted R1 and R2, and an extended source North 
of these sources. In a later Chandra observation (obs.id. 7640; 

71.1 ks; June 14, 2007) source R2 was undetectable and is likely 
associated with an AGN, while the other source R1 and the ex¬ 
tended source behaved as steady sources. Updated positional in¬ 
formation for PSR 11418-6058 from Fermi LAT timing analysis 
(see the TEMPO paramete r file in the supplementary material pro¬ 
vided in lAbdo et alj|2013h yielded a position consisten t with R1 
(see also Section on PSR 11418-6058 in iRav et alj201 Ibh . Searches 


^ Assuming a pulsed fraction of 100% for PSR J1420-6048 the measured 
Chandra spectrum predicts for the RXTE PCA with 3 PCU’s operating a 
5(T (3 (t) detection for a total observing time of 125 (45) ks in the PCA 2-22 
keV band. 
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for pulsed X-ray emission from PSRJ1418-6058 so far focussed on 
source R2 (ACIS in Continu ous Clocking m ode and XMM-EPIC 
pn in Small Window Mode: Ixg et al.l[20Q5h . while the e xtended 
source North to R1 and R2 was favoured as counterpart in iMarellil 
( I 2 OI 2 I 1 (see Chapter 4 page 154-157), who performed a spectral 
analysis on this proposed counterpart subsequently. 

We performed a period search using EPIC-pn Small Win¬ 
dow Mode data from a long 125.2 ks XMM-Newlon observa¬ 
tion performed on Februari 21-23, 2009 (obs.id. 0555700101) se¬ 
lecting events from a 15” aperture centered on Rl. We could 
not detect pulsed emission in the expected (a Fermi LAT based 
ephemeris overlapping the XMM observation period is available) 
period range, neither in the 0.3-2 keV nor in the 2-10 keV band. 


6.4 The Fermi blind-search pulsars PSR J1023-5746, PSR 
J1838-0537 and PSR J1826-1256 

PSR J1023-5746, detected in a blind search in (Fermi) LAT data, 
is located near star forming region RCW 49 and its OB asso¬ 
ciation Westerlund 2. It has a characteristic age of only 4.6 kyr 
and s pin-down luminosity of 1.1 x 10^^ erg/s dSaz Parkinson et al.l 
bold) . Archival off-axis Chandra ACIS observations revealed a 
hard X-ray source, CXOU J102302.84-574606.9, coincident with 
the Fermi LAT timing position with a typical (hard) pulsar-like 
power-law spectrum with a photon index L of — 1.2 ± 0.2 and un- 
absorbed flux of (1.2±0.3) x 10“^® erg/cm^s iSaz Parkinson et aP 
l20ld) . The detection was later confirmed in an on-axis 9.9 ks 
Chandra ACIS o bservation within the C handra Pulsar Survey 
(ChaPS) project jKargaltsev et alj l2012bh . The pulsar location 
is a lso positionally consistent with T eV source HESS J1023- 
575 ta.E.S.S. Collaboration et al.l201lh . Unfortunately, the current 
(hard) X-ray archive does not contain suitable observations to per¬ 
form period searches in order to reveal the possible/expected pulsed 
fingerprint of the pulsar. 

PSR J 1838-0537 (r ^ 5 ky r; L,i ~ 5.9 x 10®® erg/s) was dis¬ 
covered bv IPletsch et al.l 1 I 2 OI 2 I) as Fermi LAT blind-search pulsar 
applying a novel data-analysis search technique. Its sky position 
coincides with the extended TeV source HESS JI84I-055, which 
is likely composed of multiple sources. Recently, a 43.3 ks on-axis 
XMM-Newton observation (obs.id. 0720750201; Oct. 14, 2013) has 
been performed. The EPIC-pn exposure was split in two chunks of 
24.7 and 12.5 ks, while it was operating in Full Frame mode (time 
resolution 73.4 ms, which is just too low to catch the 145.77 ms 
timing signal). The second chunk was heavily contaminated by soft 
proton flares (only 3.7 ks of exposure left after screening) and sub¬ 
sequently ignored in our further analysis. The first exposure part, 
however, was unpolluted, and has an effective exposure of 22.1 ks. 
In this exposure we detected at a 6 cr level in the 2-10 keV band a 
(very) weak source positionally consistent with PSR J1838-0537. 
We found 61±13 counts for the source in the 2-10 keV band, while 
below 2 keV no detection could be claimed. The count rate, how¬ 
ever, is too low to determine and constrain the X-ray spectral char¬ 
acteristics of the pulsar. Unfortunately, there are no further suitable 
X-ray observations in the (HEASARC) archive providing sufficient 
time resolution to search for a possible timing signal of this pulsar. 

Another Fermi LAT blind-search pulsar is PSR J1826-1256 
(r ~ 14.4 kyr; Lsd — 3.6 x 10®® erg/s; lAbdo et al.ll2009 j) . not 
to be confused with the nearby energetic pulsar PSR J1826-1334 
(PSR B1823-13; r ~ 21 kyr; ~ 2.8 x 10®® erg/s) about 38' 
away, which has not be en detected (yet) as h igh-energy 7 -ray pul¬ 
sar contrary to its PWN jCrondin et ^J20ll]) . PSR 11826-1256 co¬ 
incides with the Northern extension of the extended TeV source 


HESS 11825-137 jAharonian etaDl2006ch . This TeV structure is 
positionally consistent with CGRO EGRE T source (> 1 GeV) GeV 
11825-1310 (also 3EG 11826-1302; see [Roberts etZll2001bh . A 
plausible soft X-ray counterpart is ASCA GIS source AX 11826.1- 
1257 located within the error contours of GeV 11825-1310. The X- 
ray source region was later resolved by Chandra (obs.id. 3851; 15.2 
ks; ACIS-I; Feb. 17, 2003) in a point-like component, positionally 
consistent with the pulsar, with a faint, ~ 4' l ong trail of hard (> 2 
keV) X-ray emission, resembling an “Eel” dRoberts et aHUOOTh . 
Using data from a much deeper more recent Chandra ACIS-I ob¬ 
servation (obs.id. 7641; 74.8 ks; luly 26, 2007) the spectral char¬ 
acteristics of the counterpart could be determined yielding a (very 
hard) photon index L of —0.79 ± 0.39 and an unabsorbed 0.3-10 
keV flux of (1.12 ± 0.25) x 10“®® erg/cm^s absorbed through a 
Hydrogen column Nh of (1.26j ln ®®) x 10^^ cm“^ (see Chapter 4, 
page 202-204 ofiM arellill201^ . In the X-ray archive we found a 
92.9 ks RXTE observation (obs.id. 90069; Nov. 24-29, 2004) with 
PSR 11826-1256 5(7 off-axis. A subsequent search for a possible 
timing signal in the 2-30 keV band yielded a negative result, not 
surprising given the source flux and its prediction for the PCA in¬ 
strument requiring observation times in the range 450-650 ks for a 
5 (t detection. A decisive answer about the X-ray timing properties 
will be given soon from the data analysis of a 140 ks XMM-Newton 
observation performed on Oct. 11, 2014. 


6.5 Young and energetic pulsars not detected by Fermi LAT: 
PSR J1301-6305, PSR J1341-6220, PSR J1614-5048, PSR 
J1803-2137 and PSR J1856+0245 


Table 13 of lAbdo et alj ) l2013h lists 28 pulsars not detected by the 
LAT with spin-down powers exceeding 10®® erg/s. Apart from PSR 
11747-2809 (in GO.9-1-0.1) all top 12 pulsars listed in this table, 
those with the highest spin-down power, are detected as soft 7 -ray 
pulsars (note; PSRB0540-6 9 has now been detected also as Fermi 
pulsar; iMartin et al.ll2014l) . Of the additional 15 pulsars we con¬ 
sider those pulsars which are young and/or associated with TeV 
counterparts the best candidates to detect as soft 7 -ray pulsars. 

In particular, PSR 11301-6305 (r ~ 11 ky r; ~ 1.7 x 
10®® erg/s), associated with HESS 11303 -631 jAharonian et al.1 
l2005d;lH.E.S.S. Collaboration et al.ll20i2^ . and PSR 11341-6220 
(PSR B1338-62; r ~ 12.1 kyr; ~ 1. 4 x 10®® erg/s), as¬ 
sociated with SN R G308.7-H0.0/G308.8-0.1 dCaswell et aLlJll992l; 
iKaspi et^l 19921) are good candidates to harbour a soft 7 -ray pul¬ 
sar. For the latter pulsar, a frequent glitcher, now also a weak X-ray 
counterpart is visible (this work) in archival XMM Newton data 
(obs.id. 0301740101; luly 26, 2005, 37.8 ks). 

Other interesting candidates are the 7.4 kyr old pulsar PSR 
11614-5048 (PSR B1610-50; ~ 1.6 x 10®® erg/s) show- 

ing weak > 2 keV X-ray emission in AS CA GIS images 

I Kawai & Tamural 19961 ; iBrinkmann et alJll999l) . but doubted by 
IPivoyaroff et al.l )l200(]|) . and PSR 11803-2137 (PSR B1800-21; r ~ 
15.8 kyr; Lsd — 2.2 x 10®® erg/s) possibly as sociated with HESS 
11804-21 6 jKargaltsev et aO l2007bl; llln et al J |20 1 3h . For the lat¬ 
ter pulsar [iCargaltsev et al.l l l2007ah detected weak, but hard pulsar 
emission in a 30.2 ks Chandra ACIS observation. 

Noteworthy is also the 81 ms pulsar PSR 11856-1-0245 (r ~ 21 
kyr; Lsd — 4.6 x 10®® erg /s) discovered in the Arecibo PALFA 
survey dHessels et al.l 120081) coincident with TeV source HESS 
11857-1-026. At X-rays the pulsar was clearly detected in a 39 
ks Chandra ACIS-I observation (obs.id. 12557; Feb. 28, 2011 
[Rousseau et alj[^12t) . but the low number of source counts pre¬ 
vented an accurate spectral characterisation of the pulsar at X-rays. 
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Table 2. Rotation powered (non-recycled) pulsars with (securely detected) pulsed emission in the hard X-ray band (> 20 keV). 


name 

period 

(ms) 

age 

(kyr) 

10 log(Lsd) 

Sl400 

(mJy) 

pulse 

shape 

photon 

index 

Fermi 

LAT/Pulsed 

TeV 

PWN/Pulsed 

PSRJ0205-t6449 (3C58) 

65.7 

5.4 

37.43 

O.O 45 I 

two sharp pulses 

-1.1(1) 

yes 

yes 

PSRJ0534-H2200/B0531-I-21 (Crab) 

33.5 

1.23 

38.66 

14 

two pulses 

curved 

yes 

yes/yes 

PSR J0537-6910 (N157B in LMC) 

16.1 

4.9 

38.69 

< 0.062 

single, sharp 

-1.57(1) 

no 

yes 

PSR J0540-6919/B0540-69 

50.5 

1.7 

38.18 

0.106® 

structured broad 

curved 

yes 


(N158A in LMC) 

PSRJ0835-4510/B0833-45 (Vela) 

89 

11 

36.84 

1100 

multiple sharp 

- 1.1 

yes 

yes/yes 

PSRJ1400-6325 IGRJ14003-6326 

31.2 

12.7 

37.71 

0.242 

single, broad 

-1.95(4) 

no 

no 

PSR J1513-5908/B1509-58 

150 

1.6 

37.26 

0.94 

single broad 

curved 

yes 

yes 

(MSH 15-52) 

PSRJ1617-5055 

69.0 

8.0 

37.20 

0.54 

single, broad 

-1.42(2) 

no 

? 

PSRJ1640-4631 (G338.3-0.0) 

206 

3.4 

36.64 

< i.qs 

single 

-1.3 

no 

yes 

PSRJ1811-1925 (Gll.2-0.3) 

65.0 

24.0 

36.81 

< 0.072 

single, broad 

-1.11(1) 

no 

no 

PSRJ1813-1246 

48.1 

43.0 

36.80 

< 0.017® 

two pulses 

-0.85(3) 

yes 


PSRJ1813-1749 (G12.82-0.02) 

44.7 

5.6 

37.75 

< 0 . 01 '^ 

single, broad 

-1.30(3) 

no 

yes 

PSRJ1838-0655 AXJ1838.0-0655 

70.5 

23.0 

36.75 

< 2 

structured broad 

-1.12(1) 

no 

yes 

PSRJ1846-0258 (Kes 75) 

324 

0.72 

36.91 

< 0.027® 

single, broad 

-1.20(1) 

no 

yes 

PSR J1849-0001 IGR J18490-0000 

38.5 

42.8 

36.99 

< 0.853® 

single, broad 

-1.37(1) 

no 

yes 

PSR J1930-H1852 (G54.1-I-0.3) 

136 

2.9 

37.08 

O.O 45 I 

single, broad 

-1.21(1) 

no 

yes 

PSRJ2022-H3842 (G76.9-I-1.0) 

48.6 

8.9 

37.47 

0.067“ 

two sharp pulses 

-1.20(2) 

no 

no 

PSRJ2229-H6114 (G106.6-H2.9) 

51.6 

10.5 

37.34 

0.25 

two pulses 

-1.11(3) 

yes 

yes 


Notes. 

Column 1 gives source name(s) and associated SNR or PWN (between brackets), when applicable. 

Column 3 gives the characteristic age (t = —O.bv/v) 

Column 4 the spin-down power (Lsj = Att^Iuu), in erg/s. 

Column 5 gives the radio flux density (or upper limit) at 1400 MHz (S 1400 ; see Section 6), taken from the ATNF database except for the noted entries, where: 
(1) Camilo et al. (2002) . f2) Crawford et al. (1998) . (3) Manchester et al. (1993) .f4) Kaspi et al. 11998) . 15) Castelletti et al. 12011) . (6) Ray et al. 1201 lb) . 
(7) Ifalgeme^^QO^, (8) Archibald^^^QOO^, (9) Pmide^et^QOO^, (lO)^rzoumaniane^^ (20^. 

Columns 6 and 7 give for the 2-150 keV band a description of the pulse shape and the energy spectrum (power-law photon index or curved spectrum). 

Column 8 indicates whether the Fermi LAT detected a pulsed signal. 

Finally, column 9 indicates whether a PWN and pulsed emission have been detected at TeV energies. 


Data from an archival 55.3 ks XMM-Newton observation (obs. id. 
0505920101; March 27, 2008) have been analysed in this work ap¬ 
plying our Maximum Likelihood source fitting method per energy 
band to disentangle the source contributions from the pulsar and 
an unrelated source at an angular distance of about 20 ”, rendering 
about 300 pulsar counts above 2 keV for further spectral analysis. 
Adopting an absorbed power-law model we obtained a very hard 
photon index L of —0.34 ± 0.08 absorbed through a Hydrogen 
column Nh of (I.SIq 7 ) x 10^^ cm“^, consistent with the Galactic 
value of (1.5-1.8)x 10^^ cm“^. The unabsorbed 2-10 keV flux was 
(9.52ln'an) X 10~^'^ erg/cm^s, consistent with the value given in 
[Rousseau et al.l ^20121) . Our results are , howe ver, somewhat differ¬ 
ent from those reported by iNice et al] ( l2013h . who found a softer 
photon index and a somewhat lower unabsorbed 2-10 keV flux, 
mainly due to the higher absorbing Hydrogen column density Nh 
of (4.9t|j) X 10^^ cm“^ that they derived. 


6.6 Another soft 7 -ray pulsar candidate: IGR J11014-6103 

A puzzling so ft 7 -rav pulsar ca ndidate is INTEGRAL source IGR 
J1101 4-6103 l lBirdetal.l[20l^ . studied in detail by IPavan et"^ 
( I 2 OIII) using all available archival X-ray data of the region. These 
authors concluded that the source is probably a PWN surrounding 
an energetic pulsar given the morphology of the source composed 


of two distinct sources separated by ~ 22 ” and a dimmer elongated 
structure (tail) of ~ 4' size. 

A subsequent Chandra observation shed light on the spatial 
structure of the emission components, and showed that the putative 
pulsar moves away from SNR MSH 1 1-61 A, located 11'away from 
IGR H1014-6103 jPavan et al.ll2014h . The modulated main jet and 
counter jet structures are perpendicularly oriented with respect to 
the propagation direction away from the SNR. Spectral analysis 
yielded for the pulsar counterpart, adopting an absorbed power-law 
model, a hard photon index Tof—1.1±0.2 with a 2-10 keV flux 
of (6.1 ± 0.6) X 10“^® erg/cm^s. 

The pu t ative pulsar has been detected recently by 
iHalpern et al.l i2014l) . who found 62.8 ms pulsations in two 
different XMM-Newton observations performed on July 21, 2013 
and June 8 , 2014. From the timing characteristics they deduced a 
spin-down luminosity Lsd of ~ 1.36 x 10®® erg/s and a character¬ 
istic age r of 116 kyr, rather old for such a system. These values 
should be taken with some care, because it had been assumed that 
no glitch or glitches had occured between the two observations. 

In this work we have also determined the pulsed spectrum 
of IGR 111014-6103 (analoguous to the method outlined for 
PSRJ1813-1749 in Sect. 15421 ). employing EPIC-pn data from 
the two XMM-Newton observations. We found for spectrum of the 
pulsed emission a photon index T of —1.13 ± 0.11 and an un¬ 
absorbed 2-10 keV flux of (3.93 ± 0.35) x 10“^® erg/cm^s, ab- 
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sorbed through an Hydrogen column Nh of (O.QSIq'j®) x 10^^ 
cm~^. T he latter value is c onsistent with the column density de¬ 
rived by IPavan et alJ ll2014li . who used Chandra ACIS-I data, for 
the total pulsar emission. From the pulsed and total (unabsorbed) 
2-10 keV fluxes we derived a pulsed fraction of about 64%. This 
demonstrates that IGR J11014-6103 is an excellent candidate to 
show up as soft 7 -ray pulsar in the future, given its spectral hard¬ 
ness and the already established soft 7 -ray nature for its total source 
emission. 


7 SUMMARY & DISCUSSION 

In section 5 we presented in detail the multi-wavelengths charac¬ 
teristics of the currently 18 members of the soft 7 -ray pulsar pop¬ 
ulation, and provided for 14 of them new high-energy results from 
this work. For PSRB0531-H21 (PSRJ0534-I-2200), PSRB0540-69 
(PSRJ0540-6919), PSRJ0205-t6449 and PSRJ1846-0258 we re¬ 
viewed the current status at high energies referring to up-to-date 
literature. 

For all of them the pulse profiles have been detected signifi¬ 
cantly above 20 keV and the spectral shapes of the pulsed-emission 
spectra have been determined (curved or power-law shapes). How¬ 
ever, these pulsars are weak emitters in the hard X-ray window, 
with PSRJ0534-H2200 (the Crab pulsar) and PSR J1513-5908 be¬ 
ing the exceptions. By exploiting the extensive RXTE and INTE¬ 
GRAL archives for accumulating sufficient exposures on these pul¬ 
sars at hard X-rays, we could achieve this major progress. Phase- 
coherent ephemerides were required to find the hard X-ray pulsed 
signals for all observations used. An extreme example is the case 
of PSR J0537-6910 for which we derived and used 30 ephemerides 
spread over 7.7 years to arrive at a lOfr detection of the pulse 
profile between 15 and 50 keV, in this case in RXTE HEXTE data. 
In several cases we also obtained new or statistically improved re¬ 
sults at lower X-ray energies using archival XMM-Newton data. 

An interesting example is PSR J0835-4510 (the Vela pulsar), 
for which we present in Fig. |4] pulse profiles from the radio to the 
high-energy 7 -ray band. High statistics INTEGRAL ISGRI, RXTE 
PCA and XMM-Newton EPIC-pn data reveal a plethora of sharp 
pulses, which have counterparts at lower energies in the radio, op¬ 
tical and/or UV bands as well as at the high-energy 7 -rays. One can 
identify five, possibly more sharp pulses, exhibiting very different 
variations in shape and intensity over the electromagnetic spectrum. 
This complexity has not yet been addressed in theoretical works, 
and it is beyond this paper to address this unique behaviour of the 
Vela pulsar. 

A pulsar, for which we present an extensive study is 
PSR J1849-0001, discovered as INTEGRAL source IGR J18490- 
0000 and associated with the TeV source HESS J1849-000. The 
latter association suggested the pulsar to be located in a PWN. 
Using again archival RXTE PCA and XMM-Newton data, we de¬ 
rived a new high-statistics spectrum for the pulsed emission ( 2-10 
keV) that can be well described with a power-law shape with in¬ 
dex of -1.37 ± 0.01 and pulsed flux of (4.21 ± 0.08) x 10"^^ 
erg/cm^s. The latter flux is a factor about four higher than re¬ 
ported by iGotthelf et al.l 1 I 2 OI ih . This pulsed spectrum could be 
extended up to 60 keV by analysing RXTE HEXTE data. Using 
INTEGRAL ISGRI data, we obtained a new spectrum for the to¬ 
tal emission in the energy band 20-300 keV with a softer spec¬ 
tral index of —1.82 ± 0.14. Taking the latter spectrum into ac¬ 
count, the spectrum of the pulsed emission is better described with 
a curved spectral shape from 2 keV up to ~ 300 keV (see Fig. 


\22\ . We applied for CXO observations (with the HRC-S for timing 
and the ACIS-S for spatially resolved spectral analysis) to resolve 
the contributions from the point source and the PWN. These HRC- 
S and the RXTE PCA observations both give a pul sed fraction of 
80 %, about three times the value published by IGotthelf et al.l 
( I 2 OIII 1 . The ACIS-S observation reveals for energies 2-10 keV 
structured diffuse emission in a rectangular region (75” x 37('5 
box in Fig|21b around the point source with a photon spectral index 
of —1.13 ± 0.06, to be compared with an index of —1.75 ± 0.05 
that we measured with XMM for a more extended region. This dif¬ 
ference we explain as Synchrotron cooling towards the outskirts of 
the PWN. 

7.1 Characteristics of the soft 7 -ray pulsar population 

In the following we will summarize and discuss the characteris¬ 
tics of the population of soft 7 -ray pulsars in comparison with 
those of the sample of non-recycled pulsars in the Fermi LAT sec¬ 
ond pulsar catalogue. Table summarizes the main characteristics 
of the 18 soft 7 -ray pulsars. They are all fast rotators with spin- 
periods (P = I/v) between 16.1 ms and 324 ms, most of them 
(14) even spinning well below 100 ms (see column 2 of Table 
H. The characteristic ages (r = —0.5tz/i>) of this pulsar sample 
range from 1.23 kyr (PSRB0531-I-21) up to 43 kyr (PSRJ1813- 
1246), thus representing a very young population, well below 50 
kyr. Moreover, the population is very energetic with spin-down 
powers (Lsd = Ttt^/zzP; I is the moment of inertia of the neutron 
star assumed to be 10'*® g cm^) above 4.4 x 10®® erg/s (the lowest 
value measured for PSR J1640-4631). 

The radio flux densities (or upper limits) at 1400 MHz (5i4oo) 
are shown in column 5 of Table taken from the ATNF Pul¬ 
sar Catalog, except for the noted entries. When there are no pub¬ 
lished values at 1400 MHz, we extrapolated to S'1400 from mea¬ 
surements/upper limits at other frequencies, assuming Sv oc v°‘ 
with a = —1.7 or using a published a value for a particular pulsar. 
For 8 members no radio pulsations have yet been detected, how¬ 
ever, two (PSR J1640-4631 and PSR J1849-0001) have upper limits 
at 1400 MHz as high as ~ 1 mJv. lAbdo et al.U2013h defined for the 
Fermi LAT catalogue a pulsar as “radio-loud” if S'1400 > 30/i Jj/, 
and “radio-quiet” if the measured radio flux density is lower. In 
their sample of 77 young non-recycled pulsars 41 (53%) are radio- 
loud, and the fraction of radio-loud pulsars increases with increas¬ 
ing Lsd (consistent with 100% for Lsd > 10®^ erg/s). Adopting that 
definition, in Tablel^there are three radio-quiet pulsars PSR J1813- 
1246, PSR J1813-1749 and PSR J1846-0258 with upper limits be¬ 
low 30/rJi/. There are also three weak radio-loud pulsars listed 
with flux densities 30/r Jy < S1400 < lOO/rJj/ and two with up¬ 
per limits in that range. Only the Vela and Crab pulsars have flux 
densities above 1 mjy (11%). The sample of LAT radio-loud non- 
recycled pulsars has only th ree entries below I QQiiJy and 18 (44%) 
above 1 mJy (see Fig. 3 of lAbdo et al.l ( l2013h ). This indicates that 
the radio-loud soft 7 -ray pulsars appear to be on average weaker 
radio pulsars than their LAT counterparts. However, the fraction of 
radio-loud pulsars of the energetic (Lsd > 4.4 X 10®® erg/s) soft 
7 -ray sample is similar to that fraction of the LAT pulsars, namely 
in the range 56 % to 83%, depending on how many of the so far not 
detected pulsars are radio-quiet, and seems to agree with the trent 
reported above for the Fermi LAT pulsars. 

Another noteworthy finding is that only 7 of the 18 members 
of the soft 7 -ray pulsar population are detected as 7 -ray pulsars 
in the hard 7 -ray band above 100 MeV (see column 8 , with as 
latest addition the Fermi detection of PSR J0540-6919, the twin 
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Pulse Phase 


Figure 27. The pulse shapes of the soft 7 -ray pulsar population members: a) PSR J0205-H6449, PCA 8.2-27.8 keV, b) PSR J0534-I-2200/B0531-H21, ISGRI 
20-300 keV, c) PSRJ0537-6910, PCA > 10 keV, d) PSR J0540-6919/B0540-69, HEXTE 12^8 keV, e) PSR J0835-4510/B0833-45, ISGRI 20-300 keV, 
f) IGR114003-6326, PCA 7.6-32.2 keV, g) PSR I1513-5908/B1509-58, ISGRI 20-300 keV, h) PSR 11617-5055, PCA 8.2-32.0 keV, i) PSRI1811-1925, 
ISGRI 15-135 keV,y) PSR J18I3-1749, PCA 8.1-27.5 keV, k) AXI1838.0-0655, ISGRI 20-150 keV, 1) PSRJ1846-0258, ISGRI 20-150 keV, m) IGRI18490- 
0000/PSR 11849-0001, HEXTE 14-150 keV, n) PSRI1930-H1852, PCA 8.2-32.1 keV, o) PSR 12022-1-3842, PCA 11.2-27.3 keV, p) PSR12229-H6114, PCA 
9.4-22.4 keV, and the recently reported hard X-ray pulsars q) PSR 11640-4631, AtiSTAR 3-25 keV, r) PSR11813-1246, XMM EPIC-pn 8-12 keV (significantly 
smeared due to the rather coarse time resolution of ~ 5.7 ms given its period of 48.1 ms; see also Fig. ED. 
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E [MeV] E [MeV] 

Figure 28. The high-energy pulsed emission spectra for 17 of the 18 soft 7 -ray pulsar population members from 0.1 keV - 10 GeV spread across two panels. 
The spectra of PSRB0531-1-21 (Crab), PSRB0833-45 (Vela) and PSRB1509-58 are shown in both panels for reference purposes. The names of the soft 7 -ray 
pulsar members, ordered according to their right ascension increasing from bottom to top, are indicated in the plots along with color coding. Symbols/data 
points are for flux measurements and solid lines (filled regions for the Fermi passband) for the model fits. The pulsed spectrum of PSRJ1640-4631 is not 
shown in either of the panels, because only its total spectrum has been reported on. 


of the Crab pulsa r in the LMC: [Martin et al.|l2014t) . Furthermore, 
lAhdo et ^ 1 I 2 OI 3 I) list in their Table 13 28 rotation-powered pul¬ 
sars which are not (yet) detected hy Fermi despite of their high 
spin-down powers (above 10®® erg/s). Of the top 12 in that table, 
those with the highest spin-down power, 11 are listed in Table of 
this paper. Interestingly, those most energetic pulsars missed by the 
LAT, are reported by us to be soft 7 -ray pulsars. 

In the last column of Table one can see that at least 
12 members do have, often bright, TeV counterparts, associated 
with the TeV emission from their PWNe. For three of these pul¬ 
sars (PSRI1813-1749 , PSRI1838-0655 and PSR11849-0001) we 
present in this paper the total (PWN-l-pulsar) source spectra as mea¬ 
sured with the imaging instrument INTEGRAL ISGRI from 20 
keV up to 150 keV, together with the pulsed-emission spec¬ 
tra. The difference between these two spectra gives an estimate of 
the PWN spectrum at hard X-rays. This PWN spectrum is shown 
for PSRI1838-0655 in Fig[T 8 ]for energies between 3 and ~ 60 
keV, and suggests a break/bend around 50 keV. For PSRJ1617- 
5055 and PSR J1811-1925 the latter located in SNR G11.2-0.3, we 
also present the pulsed and total-emission spectra, but these pul¬ 
sars are not detected at TeV energies, so far. In only two cases, 
PSRJ0534-H2200 (Crab) and PSR10835-4510 (Vela), pulsed TeV 
emission has been detected. 

The variation in morphology of the pulse profiles (described 
in column 6 of Table in the hard X-ray band can be seen in 
Fig.|27] This figure shows that the large majority (11 members) 
has broad/structured single pulses, and one (PSR J0537-6910) has 
a sharper single pulse, making the fraction with single pulses 67%. 
These are all, except the canonical soft 7 -ray pulsar PSRJ1513- 
5908 and PSR10540-6919, not detected at high-energy 7 -rays (> 
100 MeV) by Fermi LAT. The remaining 6 (33%) have double or 
even multiple (PSR 10835-4510; Vela pulsar) pulses, of which 5 are 


detected at high-energy 7 -rays (except so far PSR12022-I-3842). It 
is interesting to note that of the 77 non-recycled pulsars in the sec¬ 
ond Fermi pulsar catalogue 58 (or 75%) show two strong, caustic 
peaks significantly separated, thus a very much larger fraction than 
for the soft 7 -sample. In fact, 91% (10/11) of the soft 7 -ray pul¬ 
sars which are not detected by Fermi LAT show a single, mostly 
broad/structured pulse. 

Finally, column 7 of Table|^gives information about the spec¬ 
tral shape at hard X-rays. All members, except PSR 10534-1-2200, 
PSRJ0540-6919 and PSRJ1400-6325, show very hard spectra in 
the soft 7 -ray band with photon indices in the range -0.85 - -1.57. 
Fig. 1^ shows in two panels the spectral energy distributions over 
the high-energy band from ~ 2 keV up to 10 GeV of the pulsed 
emissions of the soft 7 -ra y pulsars. The pre liminary Fermi spec¬ 
trum of PSR 10540-6919 dMartin et alj|2()14h is not included, but 
the pulsed spectrum over the total 1 keV - 10 GeV band mimics 
the shape of that of PSR10534-H2200. Also PSR11640-4631 is not 
included, because spectral information on the pulsed emission has 
not yet been reported for energies above 10 keV For reference, 
the high-energy pulsed spectra of PSR 10534-1-2200, PSR 10835- 
4510 and PSR11513-5908 are shown in both panels. PSR10835- 
4510 exhibits a pulsed spectrum that is ’typical’ for the Fermi 
LAT-detected pulsars with maximal luminosity at GeV energies. 
PSR11513-5908 has also been (weakly) detected by Fermi above 
100 MeV, but reaches its maximum luminosity at MeV energies. 
It is clear from this figure that the pulsed emission spectra of the 
soft 7 -ray pulsar population, given the hardness of the spectra at 
hard X-rays/soft 7 -rays and the non-detections above 100 MeV (11 
members) or the detections as soft Fermi LAT pulsars (5 members), 
also reach maximum luminosities in the MeV range. This is thus 
remarkably different from the typical spectral characteristics of the 
Fermi LAT-detected pulsars. Note, that PSR10835-4510 has been 
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Figure 29. The disteibutions of the soft 7 -ray (dark grey; 18 members) and hard 7 -ray (light grey, Fermi LAT) pulsar populations as a function of spin-down 
luminosity (left panel; including only non-recycled Fermi LAT pulsars; 77 members), characteristic age (middle panel) and surface (dipole) magnetic field 
strength (right panel). The means (vertical) and standard deviations (horizontal) for both populations are indicated by the solid lines. The recycled pulsar subset 
(43 members) of the Fermi LAT population has been ignored in the estimation of these numbers. The soft 7 -ray pulsar population is on average ~ 50 times 
more luminous (left panel) and ~ 10 times younger than the Fermi LAT pulsar population, while the Bs distributions are comparable. 


detected in the soft 7 -ray band, only, because its distance is ~ 290 
pc. 

We also compared the distributions of the spin-down power, 
characteristic age and surface magnetic field strength at the pole 
(Bs ci 3.2 X 10^®a/— iz/iz® G) of the soft- and hard 7 -ray pulsar 
populations. The results are shown in the three panels of Fig. 
for spin-down power (left), characteristic age (middle) and surface 
magnetic field strength (right). The averages (vertical lines) and 
standard deviations (horizontal lines) of the log of the spin-down 
powers for the Fermi LAT non-recycled (77 members; light grey 
histogram left panel Fig.l29t and soft 7 -ray (18 members; dark grey 
histogram left panel Fig. I29t pulsar population are 35.73 ± 1.04 
and 37.36 ± 0.63, respectively. Thus the latter is on average ~ 
43 X more energetic than the Fermi LAT population. The averages 
and standard deviations of the log of the characteristic ages are 
4.78±0.75 and 3.81±0.52 for the non-recycled Fermi LAT pulsar- 
and soft 7 -ray pulsar population, respectively. This indicates that 
the soft 7 -ray pulsar population is on average ~ 9.3x younger 
than the non-recycled Fermi LAT pulsar population. Finally, the 
right panel of Fig. shows that the distributions of the surface 
magnetic field strengths are comparable with averages log Bs of 
12.36 ± 0.37 and 12.52 ± 0.49 for the non-recycled Fermi LAT 
pulsar- and soft 7 -ray pulsar population, respectively. 

7.2 Too many young luminous 7 -ray pulsars? 

The advent of the Fermi launch and the milestone later achieved 
with the release of the second Fermi pulsar catalogue triggered 
studies of the collective properties of the pulsar population for 
comparisons with predictions from competing 7 -ray pulsar mod- 


els, constraining the 

birth properties, beaming and evolution 

jGonthier et al.ll2007l; 

Watters & Romanill201 ll; iTakata et al.ll201 ll: 

IPierbattista et alj |2012 

). One of the conclusions from these studies 


is that (one-pole) outer-gap type models and slot-gap type models 
can reproduce many of the average radio and high-energy 7 -ray 
characteristics. However, significant differences in properties be¬ 
tween the simulated and the observed Fermi LAT populations are 
also noted. The sample of soft 7 -ray pulsars is too small for a sim¬ 
ilar study, and the peculiar characteristics of the young and ener¬ 
getic pulsars not detected at soft 7 -ray energies does not make it a 


fully coherent picture. Nevertheless, it is interesting to shortly ad¬ 
dress how this sample fits in the general picture emerging from the 
theoretical studies of the Fermi population of young non-recycled 
pulsars. 

In Table 2, we are dealing with very young energetic pul¬ 
sars, which are mostly not detected by Fermi LAT because they 
reach their maximum luminosities in the MeV band and become 
too weak for detection above 100 MeV. All high-energy emission 
models propose contributions to the non-thermal pulsar spectrum 
of synchrotron radiation produced deeper in the magnetosphere, 
reaching its maximum luminosity at MeV energies, and of higher- 
energy curvature radiation and inverse Compton radiation produced 
higher-up in the magnetosphere, closer to the light cylinder. Appar¬ 
ently, the contribution from synchrotron emission is in our sample 
dominating in the shape of the measured average spectrum over the 
components produced closer to the light cylinder. These curvature 
and inverse Compton components appear to be too weak for detec¬ 
tion by Fermi LAT. Some of these pulsars will likely be detected 
by Fermi after the accumulated exposures have been sufficiently 
increased over the coming years (a first example is the recently re¬ 
ported detection of PSRB0540- 69). This is particularly interesting 
because [Pierbattista et alj l l2012h conclude that their simulations for 
all models (Polar Cap, Slot Gap, Outer Gap and One Pole Caustic) 
under predict the detected number of LAT pulsars with high Lsd- 
Our sample of soft 7 -ray pulsars increases even further the number 
of pulsars with the highest Bsd that emit non-thermal gamma-ray 
radiation, thus increases the found discrepancy. 

In section]^ we discussed a sample of young luminous pul¬ 
sars which are not yet detected in the soft 7 -ray band, but are ex¬ 
cellent candidates for detection in dedicated observations with X- 
ray instruments with large sensitive areas and operating in modes 
with sufficient time resolution. NuSTAR (in the 3-80 keV band) 
and XMM-Newton (2-12 keV band) are currently the only observa¬ 
tories, which can provide the required capabilities. An increase of 
the size of the detected soft 7 -ray pulsar population and, in partic¬ 
ular, a more detailed measurement of the timing and spectral char¬ 
acteristics of individual pulsars and the sample in the soft 7 -band 
up to MeV energies will give important constraints for the above 
mentioned comparisons between the different models aiming at ex¬ 
plaining the production of high-energy radiation in the pulsar mag- 
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netospheres. It is therefore unfortunate that no future space mis¬ 
sions are currently planned to bridge the observational gap between 
100 keV and 100 MeV. 
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